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Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 
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1  Introduction 


Purpose 

The  presence  of  contaminants  in  the  near  surface  ground  water  has  resulted 
in  numerous  geotechnical  and  ground  water  related  studies.  Most  of  these 
studies  have  focused  on  some  particular  area  of  interest,  such  as  the  delinea¬ 
tion  of  a  particular  plume  or  the  geotechnical  characterization  of  a  specific  site 
for  waste  disposal. 

The  amount  of  data  required  for  characterizing  any  such  site  is  driven  by  a 
general  unwillingness  to  accept  any  degree  of  risk  and  a  lack  of  understanding 
of  the  subsurface.  Both  result  from  being  unable  to  directly  "see"  below  the 
surface. 

Obtaining  the  required  data,  such  as  from  borings  or  monitor  wells,  at  a 
contaminated  site  is  hazardous  and  costly.  Other  factors  which  must  be 
considered  beyond  die  normal  costs  associated  with  site  investigations  are,  die 
cost  of  protective  clothing  and  protective  equipment  needed  for  workers,  the 
potential  exposure  to  die  workers  of  hazardous  material  during  the  investiga¬ 
tion,  and  die  risk  of  further  environmental  contamination. 

To  obtain  the  necessary  subsurface  data,  the  site  is  usually  sampled  on  a 
grid  pattern.  This  occurs  because  most  ground  water  studies  are  conducted  on 
die  premise  that  the  aquifers  are  homogeneous  and  isotropic.  For  regional 
ground  water  studies  to  determine  yield  from  an  aquifer,  these  assumptions 
can  usually  be  used  satisfactorily.  However,  at  small  sites,  variability  in  the 
aquifer  is  critical  for  contaminant  movement.  Here  die  geometry  of  the  more 
permeable  materials  is  a  major  factor  in  the  flow  of  contaminants  in  ground 
water. 

For  example,  if  a  discontinuous  sand  body  is  discovered  at  a  site,  addi¬ 
tional  data  collection  is  initiated  to  define  the  extent  of  that  sand  body.  This 
is  commonly  done  with  more  grid  style  sampling,  on  a  closer  spacing  and  this 
"regridding"  may  go  through  several  iterations,  until  data  points  which  suffi¬ 
ciently  define  die  discontinuous  sand  body,  are  established. 

A  grid  pattern  does  not  consider  the  geology  of  a  site,  and  results  in  exces¬ 
sive  da**  points,  many  of  which  do  not  add  pertinent  information.  This  exces¬ 
sive  data  collection  increases  the  risk  of  exposure,  as  well  as  the  expense  of 
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the  operation.  A  27  square  mile'  area  at  Rocky  Mountain  Arsenal,  Colorado 
has  about  3,400  soil  borings  and  1,600  monitoring  wells  for  foe  purpose  of 
site  characterization  for  hazardous  waste  cleanup  (Duplancic  and  Buckle 
1989).  Had  an  information  based,  rather  than  grid  based,  boring  plan  been 
followed  some  of  these  might  not  have  been  necessary. 

Many  urban  areas  throughout  foe  world  are  built  on  sites  adjacent  to 
streams  and  rivers.  Thus,  foe  industrial  developments  associated  with  these 
areas  are  located  on  fluvial  deposits,  which  contain  discontinuous  sand  bodies. 
Therefore,  a  need  exists  to  bring  geology  into  foe  site  investigation  process  in 
order  to  limit  foe  number  of  data  points  needed  to  properly  characterize  dis¬ 
continuous  fluvial  sand  bodies. 

A  model  predicting  a  minimum  number  of  data  points  necessary  to  charac¬ 
terize  discontinuous  sands,  such  as  those  commonly  found  in  fluvial  settings, 
could  reduce  foe  number  of  data  locations  needed  and  thus  foe  risk  of  foe 
hazards  and  costs  of  drilling  at  contaminated  sites. 

Paleogeomorphic  features,  such  as  discontinuous  fluvial  sand  bodies,  are 
important  in  controlling  foe  movement  of  ground  water.  Detailed  geomorphic 
and  statistical  analysis  of  discontinuous  sand  bodies  can  predict  foe  type  of 
sand  body.  If  foe  sand  body  is  penetrated  by  one  or  more  borings,  and  foe 
width  of  foe  sand  body  can  be  estimated,  foe  correct  data  point  spacing  can  be 
determined.  Once  foe  shape  and  orientation  of  a  particular  sand  body  has  been 
determined,  foe  appropriate  hydrologic  or  physical  control  can  be 
implemented. 

From  foe  above  ideas,  a  predictive  model  was  developed.  The  model  was 
designed  primarily  for  characterization  of  sites  with  geotechnical  and  ground 
water  applications,  however,  it  could  also  be  adapted  to  other  uses.  For 
example,  in  oil  and  gas  exploration,  it  could  be  used  for  following  a  discon¬ 
tinuous  sand  reservoir.  It  could  also  be  used  for  following  a  discontinuous 
sand  body  for  pursuit  of  reduction-oxidation  fronts  where  mobile  metals  such 
as  uranium  may  have  been  precipitated.  Other  uses  for  which  foe  model 
could  be  adapted  are  certainly  possible. 


Scope  of  Research 

This  study  was  conducted  to  develop  a  predictive  model  for  locating  sam¬ 
ple  points  needed  to  characterize  discontinuous  fluvial  sand  bodies,  and  thus, 
minimize  data  needed  to  define  foe  sand  body,  minimize  foe  exposure  to 
hazards  and  to  reduce  foe  expense  and  time  spent  obtaining  such  data.  This 
was  accomplished  by: 


1  A  table  of  factors  for  converting  SI  (metric)  units  of  measurements  to  non  SI  units  is 
presented  page  ix. 
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a.  developing  a  rationale  for  exploration, 

b.  selecting  sites  for  testing  the  exploration  rationale, 

c.  selecting  and  applying  statistical  techniques  to  determine  locations  of 
future  sample  points, 

d.  establishing  confidence  of  data  point  locations  from  established 
geology,  and 

e.  developing  a  predictive  model. 

The  following  methodology  was  used  to  develop  the  model  to  determine 
die  optimum  number  and  location  of  data  points  needed  to  adequately  charac¬ 
terize  discontinuous  sand  bodies. 

A  rationale  for  exploration  was  developed  using  minimum  requirements  for 
die  number  of  data  points  needed,  taking  into  consideration  the  probability  of 
encountering  a  discontinuous  sand  body  and  the  site  boundary. 

A  systematic  evaluation  was  conducted  to  determine  how  the  hazard  and 
expense  of  obtaining  data  for  characterizing  a  discontinuous  sand  body  at  a 
contaminated  site  could  be  reduced.  Figure  1  shows  the  results  of  a  compar¬ 
ison  of  die  typical  method  of  gridding  Opart  A)  with  a  systematic  approach 
(part  B).  The  systematic  approach  places  data  points  around  the  perimeter  of 
a  site  and  defines  a  discontinuous  sand  body  by  predicting  its  location  from 
stratigraphic  data,  rather  than  just  gridding  and  regridding  the  entire  site.  It 
was  evident  that  bringing  geology  into  the  data  location  selection  portion  of 
site  investigations  as  soon  as  possible  would  reduce  the  number  of  data  points 
needed. 

To  conduct  the  systematic  evaluation,  how  die  optimal  data  should  be 
collected  to  reduce  die  hazard  must  be  determined.  To  determine  how  die 
optimal  data  is  collected,  the  control  of  die  scale  on  the  optimal  data  collection 
was  determined.  This  is  shown  in  Figure  1,  part  B,  as  follows: 

a.  The  control  of  the  scale  was  determined  by  the  variability  of  the  site. 

To  establish  the  variability  of  the  site,  the  environment  of  deposition 
was  determined. 

b.  The  environment  of  deposition  was  determined  by  interpretation  of 
stratigraphic  information,  such  as  sedimentary  structure,  mean  grain 
size,  and  relative  percentages  of  quartz  and  matrix  material  (Berg 
1986). 

c.  Once  die  environment  of  deposition  was  determined,  the  morphology 
was  predicted,  which  dictated  the  scale  of  the  data  collection. 
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Figure  1 .  Diagrammatic  representation  of  (A)  the  typical  method  of  defin¬ 
ing  e  discontinuous  sand  body  and  (B)  a  systematic  approach  of 
defining  a  discontinuous  sand  body 
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d.  With  the  scale  of  the  data  collection  established,  die  distribution  of  the 
sand  was  determined,  using  Lorenz  and  others  (1985)  method  for  a 
meander  belt  sand.  The  distribution  of  the  sand  controls  die  probability 
of  establishing  a  meaningful  data  point. 

e.  The  probability  of  obtaining  a  meaningful  data  point  controls  the  expo¬ 
sure  to  the  hazards  and  the  expense  of  obtaining  the  data  points. 

This  rationale  was  tested  in  study  areas,  selected  from  a  literature  survey 
and  established  criteria  met  by  the  sites.  The  sites  contained  a  specific 
geologic  condition  (a  discontinuous  sand  body)  and  a  data  base  which  had  an 
abundance  of  data  already  available.  The  data  base  itself  was  defined  to 
determine  the  amount  and  type  of  data  available  and  its  density  (which  con¬ 
trolled  the  resolution  of  interpretations). 

Such  a  site  allowed  for: 

a.  The  extent  of  die  site’s  discontinuous  sand  body  to  be  established,  such 
as  from  isopact  and  sand  distribution  maps. 

b.  Comparison  of  predicted  parameters  with  that  which  is  actually  present. 
This,  in  turn,  allowed  for  a  specific  degree  of  confidence  for  various 
data  point  spacing  to  be  determined. 

Statistical  techniques  were  selected  which  could  be  applied  to  determine  the 
probability  of  the  location  of  a  data  point  being  meaningful  and  the  degree  of 
confidence  of  that  location. 

Kriging  was  the  geostatistical  technique  chosen,  based  on  numerous 
authors’  descriptions  of  the  applicability  of  kriging  to  geologic  data  sets  (Clark 
1979,  Davis  1986,  Di  1989,  and  many  others).  Kriging  interpolates  irregu¬ 
larly  spaced  data  to  a  regular  grid,  which  was  used  for  contour  plotting. 

Kriging  compared  pairs  of  known  data  points  to  generate  a  curve,  called  a 
variogram.  The  variogram  shows  the  variation  of  die  numeric  variable  versus 
distance  from  control  points.  Kriging  estimates  die  value  of  a  variable  away 
from  control  points.  These  estimates  are  produced  in  a  grid,  which  can  be 
contoured  to  provide  a  map  of  the  estimated  variable  values.  Additional 
information  on  Kriging  is  contained  in  Appendix  A:  Kriging  Tutorial. 

Once  die  sand  body  morphology  was  estimated  from  stratigraphic  methods, 
sand  body  width  and  probability  of  dual  penetration  were  used  as  the  pairs  of 
known  data  points  to  generate  different  data  point  spacing.  This  allowed  for  a 
specific  data  point  spacing  to  be  established,  based  on  the  probability  desired, 
for  reencountering  die  sand  body. 

Also,  sand  thickness  was  contoured  and  confidence  intervals  determined  for 
die  sand  thickness.  With  abundant  data,  the  number  of  data  points  was 
varied.  Sand  thickness  and  confidence  was  determined  for  each  variation. 
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resulting  in  a  minimum  amount  of  data  sufficient  to  define  the  sand  body. 
This  minimum  data  being  based  on  die  confidence  desired. 

The  rationale  for  exploration,  stratigraphic  methods  for  predicting  extent, 
and  statistical  methods  for  predicting  data  point  locations  and  confidence  for 
those  locations  were  combined  into  a  predictive  model. 


Related  Studies 

Previous  studies  have  been  concerned  with  either,  (1)  the  methods  for 
interpretation  and  prediction  of  environments  of  deposition  (Berg  1986,  1970; 
Ethridge  et  al.  197S  and  LeBlanc  1972)  or,  (2)  statistical  predictive  methods 
(David  1977,  Davis  1986  and  Mousset-Jones  1980). 

May  (1985)  dealt  specifically  with  die  application  of  kriging  techniques  in 
conjunction  with  paleogeomorphic  predictive  techniques  delineating  the  overall 
trend  of  a  fluvial  sand  sequence  and  delineating  areas  where  additional  data 
were  needed. 

Di  and  others  (1989)  used  geostatistics  in  designing  sampling  strategies  for 
soil  surveys.  Their  study  determined  the  utility  of  geostatistics  in  assisting 
design  of  a  sampling  scheme  for  soil  morphological  properties  in  an  alluvial 
system.  Kriging  was  also  used  in  the  study  for  predicting  soil  property  values 
at  unsampled  locations. 

McBratney  and  Webster  (1983)  also  applied  kriging  to  regional  soil  sam¬ 
pling  and  suggested  that  the  actual  efficiency  achieved  in  their  studies  was 
three  -  to  nine  -  times  greater  than  that  estimated  by  die  classical  statistical 
methods.  However,  few  other  studies  of  this  kind  have  been  made  to  substan¬ 
tiate  die  claim  (Di  et  al.  1989). 
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2  Rationale  for  Exploration 


Sampling  Strategy 

The  first  step,  as  in  any  geological  investigation  is  a  literature  survey.  A 
literature  survey  should  give  at  least  an  idea  of  die  geology  at  a  particular  site 
and  thus,  whether  or  not  a  fluvial  sand  deposit  is  expected.  Literature  may  or 
may  not  give  information  as  to  any  expected  trend  (i.e.,  dip)  of  "bedrock”  on 
which  fluvial  deposits  may  have  developed  and  thus  the  expected  trend  of  any 
potential  fluvial  deposits.  In  some  areas  the  published  information  may  be 
detailed  enough  to  show  that  there  are  indeed  fluvial  sand  bodies  expected  in 
the  specific  horizon  (formation)  of  interest.  There  are,  however,  areas  for 
which  no  information  at  all  is  published,  so  that  the  first  information  would  be 
from  die  site  of  interest  itself.  This  could  be  from  surface  morphology  as  in 
the  case  of  surface  or  very  near  surface  interest,  or  from  actual  sampling  as  in 
the  case  of  below  surface  interest. 

This  model  will  be  based  on  die  premise  that  little  detail  is  known  at  the 
site  of  interest  itself.  The  only  information  that  may  be  known  is  that  there  is 
die  possibility  of  a  fluvial  sand  body  at  die  site. 

In  nearly  any  she  investigation  there  will  be  boundaries  within  which  the 
investigation  will  be  restricted,  at  least  initially.  In  contaminated  or  character¬ 
ization  sites,  for  example,  the  boundary  will  usually  be  die  property  boundary 
of  die  project  itself,  at  least  initially.  For  oil  and  gas  or  mineral  exploration, 
die  boundary  will  be  the  lease  block. 

For  any  site  investigation  there  will  be  lines  bounding  the  area  of  interest. 
To  encounter  a  fluvial  sand  body  that  may  enter  or  exit  the  site,  die  boundary 
should  be  die  first  region  in  which  to  obtain  data.  The  spacing  of  die  data 
collection  locations  will  be  based  on  die  legal  requirements  and/or  die  smallest 
size  sand  body  that  is  important  in  die  project,  or  which  can  be  predicted. 

The  spacing  requirements  can  range  from  about  one  hundred  feet  to  thousands 
of  feet.  More  discussion  about  die  predicted  size  limitations  will  come  in  later 
sections. 

In  die  case  of  contamination  and  characterization  sites,  die  order  in  which 
die  boundary  is  drilled  can  be  dictated  by  logistics,  since  a  minimum  number 
of  data  points  may  be  required,  regardless  of  die  geology.  If,  however,  there 
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it  tome  variation  possible,  due  to  the  geology  encountered  or  if  literature 
indicates  an  expected  trend,  a  specific  approach  is  more  desirable. 

There  are  two  possible  cases:  (1)  die  sand  body  is  at  or  near  die  surface, 
and  (2)  the  sand  body  is  below  the  surface.  In  either  case,  surface  geophysics 
is  a  tool  which  may  be  brought  into  use.  In  the  case  of  near  surface  contami¬ 
nation  and  characterization  sites,  resistivity,  electromagnetic,  or  ground  pene¬ 
trating  radar  may  be  used,  running  lines  along  die  boundaries  itself.  For  sites 
suspected  of  having  deeper  sand  bodies  seismic  lines  may  be  run  along  the 
boundaries. 

Each  of  these  surface  geophysical  methods  have  drawbacks  which  may 
cause  them  to  be  of  little  or  no  use.  These  drawbacks  include:  cultural  inter¬ 
ference,  which  inhibit  die  use  of  resistivity  and  electromagnetic  surveys;  and 
layers  of  different  characteristics,  which  mask  anything  below  than,  prevent¬ 
ing  ground  penetrating  radar  or  seismic  from  "seeing"  the  horizon  of  interest. 

Since  boundaries  of  many  sites  correspond  with  roadways,  fences,  utilities, 
etc.,  die  utilization  of  surface  geophysics  for  these  sites  is  often  severely 
limited.  Masking  and  resolution  needed  for  deeper  zones  limits  use  of  surface 
geophysics  in  the  case  of  many  fluvial  sand  bodies. 

If  little  or  no  geophysical  information  is  available  or  obtainable,  then  the 
locations  for  data,  such  as  from  drilling,  are  needed  in  an  order  to  optimize 
encountering  a  fluvial  sand  body. 

The  dip  of  the  bedrock  in  die  zone  of  interest  generally  will  dictate  the 
orientation.  Such  information  may  have  been  gathered  from  die  literature. 

For  a  site  with  no  information  available,  a  minimum  of  three  drill  holes  to 
define  die  general  stratigraphy  under  a  site  would  be  placed,  and  a  three  point 
problem  solved  to  obtain  the  dip  of  die  bedrock  for  die  horizon  in  question. 
These  three  stratigraphic  locations  should  be  placed  along  boundaries,  so  that 
they  can  provide  other  data  of  interest  for  the  site. 

Once  the  dip  of  die  bedrock  in  question  is  obtained,  die  placement  of  data 
locations  can  be  optimized  for  encountering  a  fluvial  sand  body.  The  applica¬ 
tion  for  which  die  model  is  used  will  dictate  die  actual  locations.  For  example 
if  a  north  -  south  trend  is  expected  in  the  fluvial  sand  body  (where  die  drain¬ 
age  control  is  only  due  to  a  southward  dipping  bedrock)  in  a  contaminated  or 
characterization  site,  the  southern  most  boundaries  would  be  drilled  first.  The 
southern  boundary  would  be  the  expected  direction  for  a  sand  body  to  exit  die 
site,  due  to  stream  flow,  assuming  the  dip  of  the  bedrock  is  or  was  the  control 
for  any  stream  developed  upon  the  bedrock.  This  would  also  be  die  direction 
for  any  contaminant  or  potential  contaminant  migration  off  site,  based  on 
expected  ground  water  flow  in  such  a  sand  body. 

Once  the  first  (south  in  the  example)  boundary  is  drilled,  die  drilling  pro¬ 
gram  should  work  from  this  first  boundary  on  both  sides  (to  the  north  along 
die  east  and  west  boundaries  in  die  example).  These  would  be  die  next  most 
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likely  boundaries  for  an  exiting  fluvial  sand  deposit,  based  on  the  bedrock 
trend. 

Finally,  foe  remaining  boundaries  (north  in  foe  example)  would  be  drilled 
in  foe  event  that  a  fluvial  deposit  may  extend  through  it.  A  sand  body  may 
not  occur  at  an  expected  boundary  for  various  reasons.  One  reason  is  that  foe 
sand  body  may  not  be  continuous  through  foe  site.  Another  reason  is  that  a 
stream  may  have  developed  under  different  controls  than  bedrock  dip,  or 
perhaps  foe  bedrock  dip  was  different  at  foe  time  the  stream  developed  than  it 
is  today. 

The  boundary  drilling  should  be  carried  to  completion,  regardless  of 
whether  or  not  a  fluvial  sand  body  has  been  encountered  early  in  foe  boundary 
drilling,  because  more  than  one  may  exist  at  a  site. 

If  no  fluvial  sand  is  found  in  foe  boundary  exploration,  drilling  should 
progress  inward  in  foe  same  manna*  as  foe  boundary  was  drilled,  because  a 
fluvial  sand  body  may  exist  within  foe  site,  but  not  extend  out  of  it.  This 
could  be  caused  by  a  facies  change,  or  a  fault  whae  a  sand  body  would  have 
been  displaced.  This  is  important  in  contamination  or  characterization  investi¬ 
gation  as  a  collection  or  potential  collection  area  within  the  site  for 
contamination. 


Environmental  Interpretations 

Once  a  sand  has  been  encountered  at  a  site,  foe  stratigraphic  information 
must  be  interpreted  to  establish  foe  environment  of  deposition,  whether  or  not 
it  is  a  fluvial  sand,  and  whether  or  not  it  is  from  a  meandering  stream.  This 
is  done  by  analyzing  sedimentary  structure,  texture,  and  composition  (Visher 
1965). 

Of  these  sources,  foe  sedimentary  structures,  are  of  most  importance, 
because  they  reflect  foe  processes  that  caused  foe  sediment’s  distribution. 

This  is  followed  in  order  of  importance  by  textural  change  and  composition. 
Otha  knowledge  which  can  be  helpful  in  interpretation  include  regional  strati¬ 
graphic  setting,  nature  of  adjacent  sediments,  types  of  associated  fossils,  and 
lateral  variations  in  foe  sandstone  (Berg  1986).  Information  can  also  be 
obtained  secondarily  from  geophysical  log  responses,  and  porosity  and  perme¬ 
ability.  From  a  vertical  sequence  through  foe  sediments  in  foe  horizon  of 
interest,  foe  general  morphology  can  be  predicted.  Figure  2  shows  common 
sequences  of  sedimentary  structure,  texture,  and  composition  from  sandstones 
of  different  depositional  environments.  A  sandstone  could  be  determined  to 
be  of  fluvial  origin  if:  (1)  foe  sequence  of  sedimentary  structures  range  from 
dune  (representing  high  flow  regime)  at  foe  base  to  ripples  (representing  low 
flow  regime)  at  foe  top;  (2)  grain  size  decreases  upward;  and  (3)  composition 
is  mostly  quartz  with  some  matrix,  foe  matrix  increasing  near  foe  top  of  foe 
sequence. 
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Figure  2.  Common  sequences  of  sedimentary  structure,  texture,  and  com¬ 
position  of  sandstones  of  different  depos'itional  environments 
(after  Berg  1986) 
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There  are,  however,  two  types  of  fluvial  deposits:  braided  and  meander¬ 
ing.  Sand  bodies  deposited  from  braided  streams  have  a  lower  sinuosity  (the 
ratio  of  channel  length  to  valley  length)  than  the  meandering  sand  body  depos¬ 
its.  Braided  channels  have  a  sinuosity  of  less  than  l.S  and  meandering  chan¬ 
nel’s  sinuosities  are  equal  to  or  greater  than  l.S  (Leopold  and  Miller  1964). 
Unfortunately,  the  depositional  sequence  remains  the  same  regardless  of  the 
curvature,  since  the  sequence  is  based  on  the  flow  regimes.  But,  braided 
sandstone  bodies  average  on  the  order  of  2000  ft  in  width  while  the  meander¬ 
ing  sand  bodies  have  much  wider  distribution,  on  the  order  of  10,000  ft  or 
more  (Berg  1986).  An  exception  to  this  is  a  meandering  stream  carrying  a 
large  suspended  load,  resulting  in  small  point  bars  because  of  sand-poor  sedi¬ 
ment.  From  one  vertical  sequence,  the  type  of  stream  (i.e.  braided  vs.  mean¬ 
dering)  may  not  be  obtained.  If  the  vertical  sequence  is  thick  (thickness  vs. 
width  relationships  will  be  discussed  in  the  next  section)  the  sand  body  could 
be  assumed  to  be  meandering. 

The  morphologies  of  sand  bodies  that  are  not  modern  are  based  on  those 
observed  in  modern-day  environments  of  deposition.  Each  environment  of 
deposition  has  a  characteristic  distribution  of  sediment  (LeBlanc  1972). 


Morphology  Estimation 

Once  determined  (or  assumed)  to  be  from  a  meandering  stream,  the  width 
of  the  sand  body  can  be  estimated  by  using  meander  belt  amplitude  to  approx¬ 
imate  the  sand  body  width.  Meander  belt  amplitude  is  calculated  from  chan¬ 
nel  widths.  The  channel  widths  are  calculated  from  channel  depths.  The 
channel  depths  being  equivalent  to  sand  thickness  (Lorenz  et  al.  1985).  This 
can  be  done  from  one  data  location  by  making  certain  assumptions.  However, 
in  estimating  the  width  of  the  sand  body,  some  assumptions  must  also  be  made 
initially.  A  primary  assumption  is  that  meander  belt  amplitude  approximates 
sand  body  width  (Lorenz  et  al.  1985).  Another  necessary  assumption  is  that 
sinuosities  are  greater  than  1 .7  (ratio  of  length  of  channel  to  down  valley 
distance).  This  is  necessary  because  with  less  than  a  1.7  sinuosity  there  is 
little  relationship  between  width  and  depth  (Leeder  1973).  Successive  sweeps 
of  a  migrating  channel  will  erode  tops  of  previously  deposited  sands,  resulting 
in  incomplete  thicknesses,  however  it  must  be  assumed  that  a  mature  meander 
belt  is  isolated  within  the  floodplain  deposits  and  that  the  thickness  observed 
represents  that  of  the  maximum  thickness. 

The  sand  body  is  assumed  to  be  deposited  as  point  bars,  the  thickness  of 
which  is  an  approximation  of  the  depth  of  the  channel  (Allen  1965).  This  is 
shown  in  Figure  3,  Part  A.  If  the  sand  body  is  sandstone,  rather  than  uncon¬ 
solidated  sand,  the  thickness  measured  must  be  corrected  for  compaction  from 
sand  to  sandstone.  Ethridge  and  Schumm  (1978)  suggest  a  10  percent  factor 
as  reasonable  for  this  compaction  factor. 
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Point  bar 


Figure  3.  Relationships  between  point  bars,  channel  width,  and  meander 
belt  amplitude  (after  Lorenz  et  al.  1985).  Wc  is  stream  channel 
width;  Wm  is  meander  belt  width. 
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Once  die  depth  is  obtained.  Lender’s  (1973)  formula 


Wt  -  6MlM 


where 


We  *  channel  width,  m 
h  =  channel  depth,  m 

yields  channel  widths.  This  formula  is  in  meters,  so  any  measurements  in 
english  units  must  be  converted  to  metric  units  before  using  the  formula. 

Once  die  width  is  calculated,  meander  belt  amplitude  is  obtained  from 
Lorenz  and  others’  (1985)  equation 

Wm  -  7.44 Wf101 


where 


Wm  *  meander  belt  amplitude,  ft 
We  *  channel  width,  ft 

The  meander  belt  amplitude  being  the  approximation  of  the  sand  body  width. 
This  is  shown  in  Figure  3,  Part  B.  This  equation  is  for  english  units,  so  the 
metric  width  obtained  from  Leeder’s  equation  must  be  converted  to  english 
units  before  being  used. 

Lorenz  and  others  (1985)  used  both  Leopold  and  Wolman’s  (1960)  and 
Carlston’s  (1965)  data  to  establish  the  relationship  between  meander  belt  width 
and  channel  width.  These  data  contained  a  range  of  meander  belt  widths  of 
approximately  10  ft,  from  physical  models,  to  near  15,000  ft,  from  natural 
streams. 

The  estimated  sand  body  width  will  be  an  average  for  sand  body  popula¬ 
tions  over  the  sedimentary  section,  and  are  subject  to  normal  geologic  vari¬ 
ability,  such  as  differences  in  thicknesses  due  to  location  within  die  stream 
channel  or  facies  changes  due  to  different  periods  of  deposition.  The  original 
channel  depths  are  assumed  in  this  model  and  the  original  data  set  relating  that 
depth  to  channel  and  meander  belt  width  contain  variabilities  also.  The 
model,  however,  produces  satisfactory  sand  body  width  estimates. 
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Prediction  of  Additional  Data  Point  Locations 


Once  die  width  of  the  sand  body  has  been  estimated  and  the  location  of  one 
point  in  the  sand  body  is  known,  the  spacing  of  another  data  point  location  is 
established  based  on  die  probability  of  encountering  die  sand  body  again. 

Data  point  locations  in  the  sand  body  and  outside  die  sand  body  are  necessary 
to  define  the  boundary  of  die  sand  body. 

Two  theorems  of  elementary  statistics  are  necessary  for  derivation  of  die 
probability  that  a  second  data  location  encounters  a  sand  body  if  one  has 
already  encountered  it.  The  two  theorems  are:  the  law  of  compound  probabil¬ 
ity  and  the  law  of  total  probability. 

The  law  of  compound  probability  is  that  die  probability  of  two  events  (say 
A  and  B)  occurring,  P(AB),  is  equal  to  the  product  of  the  probability  of  one 
event  occurring  (say  P(B)),  and  the  conditional  probability  of  the  other  event 
(A)  occurring  given  that  the  first  has  already  occurred,  or  P(AB)  = 

P(B)P(A  |  B). 

The  law  of  total  probability  states  that  if  there  are  two  mutually  exclusive 
events,  (say  C  and  D)  then  the  probability  of  either  occurring  (P(CUD» 
equals  the  sum  of  die  probability  of  each  occurring,  or  P(CUD)  =  P(C)  + 
P(D). 

Using  these  probability  statements,  Lorenz  and  others  (1985)  developed  die 
probability  function  for  dual  penetration  of  a  sand  body,  assuming  an  ideal 
case,  as  depicted  in  Figure  4.  In  this  case,  die  meander  belt  sand  body  is 
infinitely  long,  has  a  width  of  Wm,  and  die  data  location  spacing  is  w.  The 
probability  of  both  data  points  intersecting  die  sand  body  if  one  data  location 
point  is  a  distance  x  from  die  center  of  die  sand  body  is  6/x,  where  6,  in 
radians,  is  defined  on  Figure  4  and  is  a  function  of  x.  So,  P(A  |  B)  =  0/t, 
where  A  is  a  dual  intersection  and  B  is  a  data  point  intersection  a  distance  x 
from  the  colter  of  the  sand  body.  The  probability  of  a  data  point  location 
bong  within  the  distance  of  from  x  +  dx/2  to  x  -  dx/2  from  center  of  die 
sand  body  is  2dx/Wm  (dx  is  the  mathematical  notation  for  die  differential, 
where  x  is  die  variable  of  integration).  So,  P(B)  =  2dxAVm.  Symmetry 
about  die  centef  of  die  sand  body  results  in  the  factor  of  2. 

From  the  law  of  compound  probability,  the  probability  of  both  occurring  is 
P(AB)  =  P(B)P(A  |  B)  =  ZGdx/OrWm).  It  does  not  matter  where  the  data 
points  intersect  die  sand  body,  as  long  as  there  is  dual  intersection. 

If  a  data  point  intersects  die  sand  body  at  a  different  distance  than  in  event 
B,  (say  event  C)  then  the  law  of  total  probability  gives  die  probability  as 
P(ABUAC)  =  P(AB)  +  P(AC). 
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Being  extending  to  all  possible  distances 


P(DuaI  Intersection )  =  (2/irWJ)  J^Vm/2  (ddx) 


Since 


6  -  v  -  0,  -  02 


where 


0i 


cos-1  ([WJ2  -  jcJ/w)  i/z  2;  Wm/2  -  w 
0  ifx  <  WJ2  -  w 


02 


cos-1  {[WJ2  *  JC]/W)  ifx  £  w  -  WJ2 
0  ifx  >  w  -  WJ2 


that 


P(Dual  Intersection)  =  (2/xW.n)  f  0Wm/2(x  -  <t>{  -  ^)dx 

Using  these  probability  relationships,  probability  of  dual  penetration  vs. 
sand  body  width  can  be  plotted  for  various  data  point  location  spacing.  From 
die  resulting  curves,  probability  for  dual  intersection  can  be  obtained  for  any 
combination  of  sand  body  width  (as  estimated)  and  data  point  location  spacing. 
Figure  S  shows  probability  of  dual  penetration  vs.  sand  body  width  for  several 
data  point  location  spacings.  A  computer  program,  MathCAD,  was  used  to 
calculate  die  probabilities  for  the  curves  shown  in  Figure  S.  MathCAD  is  a 
computational  software  copyrighted  (1986-1989)  by  Mathsoft,  Inc.,  Cam¬ 
bridge,  Massachusetts.  The  results  of  the  calculations  are  contained  in  Appen¬ 
dix  B.  The  program  was  verified  by  comparison  of  its  computed  values  with 
those,  for  die  same  variables,  computed  by  May  (1985)  and  by  Lorenz  and 
others  (1985).  Both  May’s  and  Lorenz  and  others’  curves  are  contained  in 
Appendix  C. 

Entering  die  plot  of  sand  body  width  vs.  probability  with  a  high  probability 
produces  a  data  point  spacing  likely  to  intersect  the  sand  body  again.  Con¬ 
versely,  using  a  low  probability  value  produces  a  data  point  not  likely  to 
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intersect  the  sand  body.  In  proceeding  radially  from  die  first  intersection 
point,  it  is  necessary  to  use  appropriate  probabilities  and  resulting  radii  to  both 
intersea  and  not  intersect  the  sand  body.  In  this  manner,  the  sand  body  boun¬ 
dary  location  is  delineated. 

The  spacing  for  location  of  the  additional  data  point,  in  any  case,  takes  the 
form  of  die  radius  of  a  circle,  the  center  of  which  is  the  data  point  location 
which  previously  intersected  die  sand  body.  This  is  shown  in  Figure  6.  In 
the  case  of  boundary  drilling,  there  would  not  be  a  foil  circle  for  the  location 
of  the  additional  data  point  due  to  the  portion  of  the  circle  falling  outside  of 
die  boundary  line. 

The  location  of  the  data  point  itself  on  the  circle  (or  portion  of)  can  be 
arbitrary.  However,  it  should  be  placed  in  a  manner  as  to  lend  itself,  as 
easily  as  possible,  to  the  locating  method  used.  Certainly,  accessibility  will 
play  an  important  role  in  the  actual  placement  of  the  data  point  location. 

If,  or  when,  more  than  one  data  point  location  is  to  be  placed  on  the  circle, 
the  distance  between  the  successive  data  points  should  be  the  same  as  the 
radius  of  the  circle,  thus  keeping  the  same  spacing.  As  new  data  are  added, 
the  spacing  can  of  course  be  adjusted  as  necessary,  based  on  new  calculations. 

Once  sufficient  data  are  obtained,  such  as  from  boundary  borings,  kriging, 
as  described  in  foe  Introduction  section,  can  be  used  to  help  select  foe  location 
on  foe  circle  for  data  location  spacr  ,  which  will  give  a  high  probability  of 
being  a  useful  location.  Kriging  is  a  geostatistical  technique  that  can  be  used 
to  estimate  irregularly  spaced  data,  -ach  as  drilling  information,  to  a  regular 
grid  which  can  be  used  for  contour  plotting,  such  as  that  of  a  geologic  sur¬ 
face.  The  kriging  variance  shows  foe  areas  where  additional  data  are  needed. 
An  area  shown  with  most  error  and  overlapping  foe  circle  would  be  foe  area 
to  place  foe  next  data  point. 
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Circle  on  which 
additional  data 
point(s)  are  to 
be  placed 


Figure  6.  Additional  data  point  location  spacing  based  on  probability  of 
dual  penetration 
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To  test  die  rationale  for  exploration  described  in  the  previous  section,  a  site 
was  chosen  which  was  known  to  have  a  discontinuous  sand  body  and  which 
had  a  large  amount  of  data  already  collected  and  available.  Thus,  the  site 
would  laid  itself  conveniently  to  the  application  of  this  rationale. 

The  site  chosen  was  a  part  of  die  Rocky  Mountain  Arsenal  near  Denver, 
Colorado.  Figure  7  shows  the  general  location.  The  Arsenal  is  located  north¬ 
east  of  Denver,  Colorado,  adjoining  die  northern  portion  of  the  Denver 
Stapleton  International  Airport.  Figure  8  shows  die  Arsenal  and  its  vicinity. 
The  application  of  the  rationale  for  exploration  was  limited  to  a  four  square 
mile  portion  of  die  Arsenal  because  of  the  availability  of  usable  data  in  that 
particular  area  and  due  to  the  location  of  die  particular  sand  body  chosen  for 
this  application.  This  four  square  mile  area  is  shown  in  Figure  9.  It  is  com¬ 
posed  of  Sections  1  and  2,  Township  3  South,  Range  7  West  and  Sections  35 
and  36,  Township  2  South,  Range  7  West. 

Due  to  contamination  being  detected  in  die  groundwater  below  die  Arsenal, 
numerous  studies  have  been  conducted  to  delineate  the  route(s)  of  migration 
(May  1985).  This  has  resulted  in  a  large  quantity  of  geological  data. 

Geological  data  includes  geotechnical  reports,  geophysical  logs,  boring 
logs,  and  core  samples.  This  data  was  in  sufficient  quantity  and  of  a  quality 
that  was  useable  in  a  hypothetical  drilling  program. 

Several  discontinuous  sands  have  been  delineated  from  die  data.  One  of 
these  sands  was  chosen  as  die  sand  to  "target"  during  die  application. 


Geologic  Setting 

Rocky  Mountain  Arsenal  is  located  in  the  Denver  basin,  a  structural 
depression  occupying  a  large  area  of  Northeastern  Colorado,  Southeastern 
Wyoming  and  Western  Nebraska.  Within  the  basin  die  Lower  Cretaceous  and 
Tertiary  rocks  occupy  an  area  of  670  square  miles,  from  Greeley  in  die  north 
to  Colorado  Springs  in  die  south  and  from  the  Rocky  Mountain  front  range  in 
the  west  to  near  Limon  in  the  east  (Robson  and  Romero  1981).  The  location 
of  die  Arsoial  within  the  basin  is  shown  in  Figure  10.  The  basin  is  filled 
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Figure  8.  Rocky  Mountain  Arsenal  and  vicinity  (from  May  at  al.  1 983) 


CtMptar  3  Com  Study 


CtMptar  I  Cm*  Study 


23 


Figure  9.  Case  study  area 


Figure  10.  Rocky  Mountain  Arsenal  in  relation  to  the  Denver  basin  (after 
Robson  and  Romero  1981) 
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with  approximately  15,000  ft  of  sediments  composed  of  limestones,  sand¬ 
stones,  shales  and  conglomerates  (Robson  and  Romero  1981). 

The  stratigraphic  position  of  the  Arsenal  within  the  Denver  basin  is  shown 
in  Figure  11.  At  the  Arsenal,  the  geologic  formations  of  concern,  with 
respect  to  this  application,  are  die  Late  Cretaceous  and  Early  Tertiary  Denver 
formation  and  foe  Tertiary  and  Quaternary  alluvium. 

May  (1985)  states  that  foe  Denver  formation  is  composed  of  250  to  400  ft 
of  interbedded  clay  shale,  claystone,  siltstone,  sand,  and  sandstone.  Low- 
grade  coal,  lignite,  and  carbonaceous  clay  shale  are  also  present.  The  pre¬ 
dominant  olive-gray,  brown,  and  green-gray  colors  in  foe  formation  are 
caused  by  rock  fragments  derived  from  foe  erosion  of  basaltic  and  andesitic 
lavas.  The  lowest  elevations  for  foe  base  of  foe  Denver  formation  in  foe 
entire  Denver  basin  are  found  along  foe  eastern  side  of  foe  city  of  Denver, 
including  foe  southern  portions  of  foe  Arsenal.  The  sands  in  foe  Denver 
formation  are  generally  weakly  cemented  sandstones  or  compact  fine-  to 
medium-grained  sands.  Many  of  these  sandy  units  represent  deltaic  channel 
deposits  which  grade  laterally  and  vertically  into  silts  and  clay  shales.  The 
individual  sand  or  sandstone  layers  are  commonly  lens-shaped  and  range  in 
thickness  from  several  inches  to  as  much  as  60  ft. 

May  (1985)  describes  foe  surficial  deposits  as  having  been  deposited  pri¬ 
marily  by  ancient  streams.  In  some  areas,  a  veneer  of  aeolian  deposits  occur, 
but  for  discussion  purposes,  these  are  included  as  alluvium.  The  surficial 
alluvial  deposits  are  of  Pleistocene  and  Recent  age.  Pre-Wisconsin  deposits 
contain  alluvial  silts,  sands  and  gravels.  The  Wisconsin-age  alluvium  was 
deposited  as  glacial  outwash  from  foe  Rocky  Mountain  front  range.  The 
aeolian  deposits  were  derived  from  glacial  outwash.  The  deposits  immediately 
overlying  foe  Denver  formation  are  identified  as  foe  Verdos  alluvium  of 
Kansan  age.  The  Verdos  is  composed  of  boulders,  cobbles,  pebbles,  and 
sands  derived  from  granites  and  pegmatites  and  Cretaceous  shale.  The  Verdos 
is  up  to  100  ft  thick  in  foe  Arsenal  area.  Recent  alluvial  sand  deposits  have 
accumulated  from  several  separate  periods  of  deposition.  Figure  12  shows  foe 
typical  geology  below  foe  Arsenal. 

From  foe  several  sand  bodies  that  incised  foe  Denver  formation  in  foe 
subsurface,  foe  highest,  thickest,  and  best  delineated  one  was  chosen  as  foe 
"target"  fot  application  of  foe  rationale  for  exploration. 


Application  of  Rationale  for  Exploration 

The  location  of  foe  sand  body  chosen  as  the  "target  sand"  dictated  foe  area 
within  foe  Rocky  Mountain  Arsenal  for  foe  application  of  foe  rationale  for 
exploration.  Figure  13  shows  foe  location  of  this  sand  body  within  foe  Arse¬ 
nal.  The  quantity  of  usable  data  dictated  foe  actual  boundaries  for  this  appli¬ 
cation.  Figure  9  shows  foe  selected  site  boundaries. 
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Figure  1 1 .  Stratigraphic  position  of  the  Rocky  Mountain  Arsenal  within  the  Denver  basin  (after  Robson  and 
Romero  1981) 


Figure  12.  Typical  geology  below  the  Rocky  Mountain  Arsenal  (after 
Robson  and  Romero  198') 
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With  die  site  selected,  the  rationale  for  exploration  described  in  the  previ¬ 
ous  chapter  was  used  for  a  hypo  theoretical  drilling  exploration  program. 

Since  a  substantial  amount  is  known  about  the  geology  and  discontinuous 
sand  bodies  at  the  Arsenal,  it  was  presumed  that  only  the  regional  geology 
was  known.  From  the  regional  geology  a  possibility  existed  for  a  sand  body 
to  be  incised  into  die  Denver  formation. 

An  assumption  was  also  made  that  there  is  too  much  cultural  interference 
at  the  Arsenal  to  use  surface  geophysics  to  investigate  the  boundary  of  the 
site.  Indurated  sands  and  gravels  in  the  area  mask  geophysical  signatures  at 
depth. 

Since  the  regional  geology  shows  the  top  of  the  Denver  formation  dipping 
to  the  south,  the  southern  boundary  of  the  site  was  hypothetically  drilled  first. 
A  spacing  of  one  mile  was  arbitrarily  chosen  for  the  boundary  drilling  spac¬ 
ing,  although  this  could  be  initial  spacing  for  a  site  investigation  based  on  an 
expected  sand  body  width  or  logistically  based  on  Sections  of  the  Public  Land 
Survey.  A  Section  is  a  square  mile,  and  this  land  location  system  is  present  at 
the  Rocky  Mountain  Arsenal.  Existing  borings  located  as  close  to  this  mile 
spacing  as  possible,  were  chosen  to  represent  the  hypothetical  drilling  loca¬ 
tions.  As  in  the  case  for  these  locations,  as  well  as  others  that  follow,  the 
borings  chosen  had  to  have  usable  information,  and  thus  were  often  not  in  the 
exact  chosen  location.  This  is  acceptable,  since  logistics  in  actual  drilling 
programs  often  dictate  offsetting  actual  borings  from  desired  locations.  This 
offsetting  can  be  due  to  inaccessibility  which  is  due  to  rugged  or  swampy 
terrain,  to  buildings  and  other  structures,  or  buried  utilities. 

The  target  sand  was  not  encountered  in  any  of  the  southern  boundary 
borings.  Next,  borings  to  complete  the  western  and  eastern  boundaries  were 
chosen.  Again,  the  target  sand  was  not  encountered.  Then  the  boring  to 
complete  die  northern  boundary  was  chosen.  This  boring  did  encounter  die 
target  sand.  Table  1  contains  die  location,  boring  number,  and  "target*  sand 
thickness  for  the  borings  chosen  in  the  hypothetical  boundary  drilling  and 
subsequent  borings.  Figure  14  shows  the  locations  of  the  boundary  drilling 
borings.  These  borings  are  located  with  the  Colorado  state  grid  coordinates. 

Using  the  data  available  from  the  boring  (1228)  which  encountered  the 
target  sand,  the  environment  of  deposition  was  determined  for  that  sand.  The 
data  available  included  grain  size  analysis,  x-ray  diffraction  analysis,  geo¬ 
physical  log,  and  detailed  core  description. 

Sands  generally  occur  in  regular  predictable  sequences,  that  are  character¬ 
ized  by  vertical  changes  in  composition,  texture,  and  sedimentary  structures. 
These  ordered  sequences  contain  information  that  provides  keys  to  their 
method  of  transportation  and  deposition.  By  studying  properties  such  as 
sedimentary  structure,  mean  grain  size,  and  relative  percentages  of  quartz  and 
matrix  material,  a  sandstone  can  be  categorized  as  to  its  environment  of  depo¬ 
sition  (Berg  1986).  This  methodology  was  used  to  determine  the  environment 
of  deposition  of  the  sand  body  encountered  between  thirteen  and  thirty 
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Table  1 


Rocky  Mountain  Arsenal  Hypothetical  Exploration  Boring  Locations 
and  Sand  Thicknesses  in  Sequence  of  Hypothetical  Borings 


lection 

Qrid  Coordinate 
Location,  It 

Boring 

Number 

Other 

Number 

Target  Sand 
Thlcknaae,  ft 

Kaat 

North 

2 

2,178,931 

176,779 

1124 

SP13 

0 

1 

2,183,891 

176,445 

1155 

SP15 

0 

1 

2,187,216 

175,608 

1143 

SP16 

0 

36 

2,188,139 

180,921 

1160 

SP02 

0 

36 

2,188,353 

186,171 

1188 

E01 

0 

2 

2,178,446 

179,361 

1123 

SP08 

0 

35 

2.178,426 

185,575 

1185 

NO  6 

0 

36 

2,183,794 

185,108 

1228 

AP01 

19.0 

35 

2,181,156 

184,639 

653 

— 

0 

35 

2,183,065 

182,320 

757 

— 

10.5 

36 

2,186,014 

184,035 

758 

— 

0 

35 

2,183,045 

180,900 

1251 

AP25 

28.6 

i 

2,184,234 

177,874 

746 

— 

10.4 

2 

2,181,060 

178,469 

1153 

SP09 

0 

36 

2,186,235 

180,686 

756 

— 

0 

2 

2,181,205 

180,372 

1247 

AP21 

0 

2 

2,181,444 

176,282 

1148 

SP12 

0 

1 

2,186,789 

178,332 

1239 

API  2 

o  1 

two  feet  depth  in  boring  1228.  The  various  procedures  that  were  used  are 
summarized  in  the  following  paragraphs.  More  details  on  these  procedures 
are  given  in  the  Stratigraphic  Interpretations  section. 

Detailed  descriptions  of  the  core  from  boring  1228  were  generated  during  a 
previews  study.  These  descriptions  provided  information  on  the  lithology  of 
the  sand  body. 

In  addition  to  sieve  grain-size  analysis,  grain-size  data  was  statistically 
analyzed  to  determine  standard  deviation  and  the  relative  degree  of  sorting  of 
the  grains  by  using  grain  size  distribution  curves.  The  degree  of  sorting  was 
helpful  in  confirming  the  history  of  the  sand. 
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The  geophysical  log  for  boring  1228  was  used  to  aid  in  establishing  die 
abrupt  lower  contact  of  the  sand  body  with  die  underlying  Denver  formation. 
The  clay  shale  of  die  Denver  formation  is  easily  recognizable. 

Selected  samples  were  tested  by  x-ray  diffraction  to  determine  the  percent 
of  quartz  present.  The  percentage  of  quartz  is  used  in  the  indication  of  the 
environment  of  deposition.  Mosdy  quartz  with  a  little  clay  matrix,  which 
increases  upward,  is  typical  for  a  fluvial  sequence. 

From  this  boring  information,  which  is  summarized  in  Figure  15,  the  envi¬ 
ronment  of  deposition  was  interpreted  by  comparison  with  diagnostic  charac¬ 
teristics  for  a  fluvial  sand,  which  are  shown  in  Figure  2. 

Hie  coal  seam  in  boring  1228,  at  a  depth  of  about  50  ft,  was  deposited  in  a 
swampy  environment.  The  coal  seam  was  overlain  by  a  layer  of  olive-brown 
clay  shale  containing  abundant  organic  material  and  sandy  silt  lenses.  The 
clay  shale  which  contains  volcanic  ash  and  worm  borings  was  deposited  in  a 
low  energy  environment  such  as  a  shallow  lake.  A  fine  to  medium  grained 
non-cohesive  to  slighdy  cohesive  alluvial  sand  was  deposited  over  die  clay 
shale.  The  abrupt  lower  contact  of  die  sand  with  die  underlying  clay  shale 
indicates  that  a  stream  deposited  the  sand,  cutting  down  into  the  clay  shale. 
The  upper  finer  portions  of  the  fluvial  sand  sequence  was  removed  by  a  much 
younger  stream  which  deposited  gravel,  sand,  silt  and  clay  on  top  of  the  Den¬ 
ver  alluvial  sand  sequence. 

With  die  environment  of  deposition  determined  to  be  that  of  .a  meandering 
stream,  die  sand  body  width  was  estimated.  Channel  width  was  calculated 
from  die  sand  thickness  of  19.0  ft.  Since  the  sand  is  consolidated,  die  thick¬ 
ness  is  corrected  by  10  percent  to  20.9  ft.  This  was  then  converted  to  metric 
6.37  meters  for  use  in  Leeder’s  equation.  Solving  Leader’s  equation  gave  a 
channel  width  of  1 17.73  meters  or  386.27  ft.  Using  Lorenz’s  equation  and 
die  channel  width  of  386.27  ft,  die  meander  belt  amplitude  was  calculated  to 
be  3050  ft.  This  is  die  estimated  sand  body  width.  This  data,  along  with 
subsequent  calculations,  is  summarized  in  Table  2.  Sample  calculations  are 
contained  in  Appendix  D. 

Using  the  curves  in  Figure  5  which  were  calculated  probabilities  for  vari¬ 
ous  sand  body  widths  and  various  well  spacings,  a  well  spacing  was  chosen. 

A  well  spacing  for  a  50  percent  probability  of  encountering  the  sand  again 
was  used,  since  definition  of  the  sand  body  is  die  object  in  this  application. 
The  well  spacing  for  the  first  set  of  data,  based  on  Figure  5,  was  approxi¬ 
mately  2500  ft.  This  2500  ft.  was  used  as  die  desired  spacing.  Table  3 
contains  die  well  spacing  for  boring  1228  and  subsequent  borings  which 
encountered  the  sand. 

The  2500  ft  spacing  results  in  an  arc  or  semi-  circle  around  boring  1228, 
having  a  radius  which  is  that  of  the  well  spacing. 

Using  die  borings  hypothetically  drilled  up  to  this  point  (numbers  1 124  - 
1228  in  Table  1),  sand  thickness  was  kriged  using  Geo-EAS,  a  public  domain 
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Figure  15.  Summary  of  information  for  boring  1228 


Tallin  2 

Data  from  Estimating  Sand  Body  Width 

Thickness.  It 

Corrected 

Thickness, 

It 

Metric 

Thickness,  m 

Channel 
Width,  m 

Channel 

Width. 

It 

Approximate 

Meander  Balt 

AmpKtude.  It 

19.0 

20.90 

6.37 

117.73 

386.27 

3050 

14.8  (Average) 

16.23 

4.96 

79.64 

261.96 

2060 

19.4  (Average) 

21.34 

6.5 

121.45 

398.48 

3150 

Table  3 

Data  from  Esl 

imating  Well  Spacing 

Boring  Number 

Estimated  Sand 

Body  Width.  It 

Probability  Percent 

Approximate  Data 
Point  Spacing,  It 

1228 

3060 

50 

2500 

767 

2060 

60 

1500 

1261 

3150 

50 

2500 

computer  program  produced  by  the  U.S.  Environmental  Protection  Agency. 
The  Geo-EAS  program  was  verified  by  comparison  of  the  program’s  results 
with  results  of  Clark’s  (1979)  kriging  of  a  simulated  iron  ore  deposit.  Both 
are  contained  in  Appendix  C.  The  estimated  percentage  of  iron  and  die 
kriging  standard  deviations  for  those  estimates  produced  by  the  Geo-EAS 
program  are  similar  to  Clark’s.  The  slight  differences  can  be  attributed  to  die 
differing  interpretations  in  die  drafting  of  the  contours. 

Once  kriged,  the  error  for  the  kriged  thickness  was  obtained,  as  die  kriging 
standard  deviation.  The  kriging  standard  deviation  was  then  contoured. 

By  comparing  die  kriged  standard  deviation  with  die  well  spacing  "arc”, 
some  portions  of  the  arc  fell  in  areas  where  the  sand  thickness  standard  devia¬ 
tions  were  more  than  the  standard  deviation  of  the  whole  data  set,  and  some 
portions  of  die  arc  fell  in  areas  where  the  sand  thickness  standard  deviations 
were  less  than  the  standard  deviation  of  die  whole  data  set.  The  areas  on  the 
arc  that  were  in  the  area  where  the  standard  deviations  of  the  sand  thickness 
was  greater  than  that  of  the  whole  data  set  were  then  hypothetically  drilled. 
The  area  with  die  highest  standard  deviation  was  selected  for  die  Erst  boring, 
dioi  additional  borings  were  selected  along  the  arc  in  die  area  where  die 
standard  deviation  of  the  sand  thickness  was  greater  than  die  standard  devia¬ 
tion  of  the  whole  data  set  These  data  locations  were  separated  by  die  same 
well  spacing  as  the  distance  from  boring  1228.  Figure  16  shows  die  arc  for 
die  well  spacing  from  boring  1228,  die  kriged  standard  deviations  and  die 
locations  for  the  next  borings. 
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Figure  16.  Weil  spacing  arc,  as  heavy  line;  contoured  kriging  standard 
deviations  (in  feet)  of  sand  body  thickness;  and  additional  boring 
locations  away  from  boring  1 228 
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Hie  area  where  die  standard  deviation  is  greater  than  about  6  is  considered 
sufficiently  defined,  that  which  is  less  than  6  is  not  considered  sufficiently 
defined.  This  is  based  on  the  standard  deviation  in  relation  to  the  standard 
deviation  for  the  attire  sample  population.  If  the  specific  location  standard 
deviation  is  greater  than  the  population  standard  deviation,  then  the  thickness 
at  that  location  is  not  sufficiently  defined.  If  the  specific  location  standard 
deviation  is  less  than  die  population  standard  deviation,  the  thickness  is  con¬ 
sidered  reliable  (Gark  1979).  In  this  case  the  standard  deviation  for  die 
population  is  6.28.  This  standard  deviation  was  obtained  by  taking  the  square 
root  of  the  population  variance,  which  is  given  by  die  Geo-EAS  program. 
Table  4  contains  the  variances  given  by  die  Geo-EAS  program  and  their  calcu¬ 
lated  standard  deviation  for  this  first  set  of  borings  and  each  subsequent  set  of 
borings. 


Table  4 

Data  Sat  Sand  Thickness  Population  Standard  Deviations  for  the 
Rocky  Mountain  Arsenal 

Population 

Variance,  ft2 

Standard  Deviation,  It 

Sot  1 

39.48 

6.3 

Sot  2 

46.31 

6.7 

Sot  3 

83.99 

9.2 

Sot  4 

72.30 

8.6 

SotS 

83.15 

8.0 

Full  Data 

194.90 

14.0 

Moy'a  Data 

316.81 

17.8 

The  second  set  of  borings  were  hypothetically  drilled.  Two  did  not 
encounter  the  target  sand,  one  did.  This  information  is  contained  in  Table  1. 
Only  boring  log  information  was  available  for  these  borings,  so  confirmation 
of  die  depositional  environment  could  not  be  made.  The  depositional  environ¬ 
ment  was  assumed  to  be  the  same,  meandering  fluvial.  A  new  sand  thickness 
was  estimated  from  the  average  thickness  of  die  two  borings  which  had 
encountered  die  sand.  A  new  sand  body  width  and  subsequent  spacing  was 
calculated.  This  information  is  shown  in  Table  2. 

Thickness  was  again  kriged,  obtaining  the  standard  deviation  of  die  sand 
body  thickness,  which  was  contoured.  An  arc  from  boring  757,  which 
encountered  the  sand,  only  allowed  for  one  additional  boring  in  an  area  where 
die  standard  deviations  of  die  sand  thickness  was  greater  than  the  standard 
deviation  of  the  whole  data  set,  die  highest  standard  deviation  on  the  arc. 
Figure  17  shows  this.  Table  4  contains  die  population  standard  deviation  as 
"Set  2". 
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The  new  location  was  hypothetically  drilled.  Based  on  stratigraphic 
relationships,  it  was  determined  that  the  target  sand  was  encountered.  This 
boring,  1251,  had  sufficient  information  to  confirm  the  environment  of 
deposition.  This  information  is  summarized  in  Figure  18.  The  stratigraphic 
section  in  the  sand  sequence  is  almost  identical  to  that  described  for  boring 
1228.  A  stream  had  cut  down  into  underlying  dark  gray  silty  clay  shale. 

A  new  average  thickness  was  calculated,  giving  a  new  sand  body  width  and 
subsequently  a  well  spacing  of  approximately  2500  ft,  the  same  spacing  as 
obtained  from  the  first  data  set.  These  are  shown  in  Tables  2  and  3. 

Thickness  was  kriged  and  standard  deviation  contoured.  The  area  of 
highest  standard  deviation  on  the  arc  for  well  spacing  from  boring  1251  could 
not  be  hypothetically  drilled,  because  there  were  no  usable  borings  in  that 
area.  Three  borings  with  approximately  die  same  welt  spacing  along  die  arc, 
as  from  boring  1251,  were  chosen.  These  were  in  areas  where  die  standard 
deviations  were  greater  than  the  standard  deviation  of  die  whole  data  set. 

There  are  fewer  borings  in  die  lower  portion  of  the  site  that  penetrate  die 
target  sand,  because  it  is  becoming  deeper  with  the  direction  of  dip. 

Two  of  die  additional  borings  did  not  encounter  the  target  sand,  one  did,  as 
shown  in  Table  1.  This  is  shown  in  Figure  19. 

This  set  of  borings  was  kriged  showing  die  area  containing  the  sand  body 
to  be  defined,  but  due  to  die  variation  in  the  standard  deviation,  three  more 
boring  sites  were  chosen  in  areas  of  relatively  higher  standard  deviation  in  an 
effort  to  make  the  standard  deviation  more  uniform.  This  is  shown  in 
Figure  20. 

These  areas  are  approximately  die  same  spacing  from  several  wells  as  die 
last  derived  well  spacing. 

These  sites  were  hypothetically  drilled,  none  encountering  the  target  sand. 
These  borings  were  added  to  the  data  set  and  die  thickness  was  kriged.  The 
resulting  standard  deviations  are  fairly  uniform.  This  is  shown  in  Figure  21. 

With  die  data  used  in  die  hypothetical  drilling  program  (Table  1),  the 
kriged  thickness  was  contoured.  This  is  shown  in  Figure  22.  The  purpose  of 
contouring  the  kriged  sand  thickness  was  to  compare  with  that  already  pro¬ 
duced  by  other  studies.  Although  similar  in  some  respects,  die  test  case  was 
noticeably  different  from  that  of  May’s  (1985)  and  ESE’s  (EBASCO  1989). 
These  are  shown  in  Figures  23  and  24  respectively. 

Since  neither  used  die  Geo-EAS  program.  May’s  (1985)  and  ESE’s 
(EBASCO  1989)  sand  thickness  data  were  kriged  with  the  Geo-EAS  program. 
The  data  point  (boring)  locations  are  listed  in  Appendix  E  and  are  shown  in 
Figures  25  and  26.  The  resulting  sand  thickness  contours  are  shown  in  Fig¬ 
ures  27  and  28.  The  kriging  standard  deviations  for  these  two  data  sets  were 
also  contoured  and  are  shown  in  Figures  29  and  30. 
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DEPTH  MEDIAN  GRAIN  SIZE  LITH.  E-LOG 
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Figure  18.  Summary  of  information  for  boring  1251 


Figure  19.  Well  spacing  arc,  as  heavy  line;  contoured  kriging  standard  devi¬ 
ations  (in  feet)  of  sand  body  thickness;  and  additional  boring 
locations  away  from  boring  1251 
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Figure  20.  Well  spacing  arcs,  as  heavy  lines;  contoured  kriging  standard 
deviations  (in  feet)  of  sand  body  thickness;  and  additional  boring 
locations 
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Figure  26.  Location  of  all  data  points  used  for  kriging  at  the  Rocky  Mountain  Arsenal 


When  comparing  the  kriged  standard  deviation  to  the  sand  thickness  stan¬ 
dard  deviations,  in  both  cases.  May’s  and  all  he  data,  the  thickness  standard 
deviations  are  greater  than  die  standard  devkuon  of  the  whole  data  sets  in  the 
southern  portions  of  the  site.  This  is  due  to  the  lack  of  sufficient  borings 
penetrating  die  depth  necessary  to  encounter  the  target  sand  in  that  area  of  the 
site.  The  variances  and  standard  deviations  for  the  sample  populations  are 
given  in  Table  4. 

From  these  results,  it  would  appear  that  die  hypothetical  borings  had 
defined  die  sand  body  better  than  May’s  (1985)  data  or  all  die  available  data 
from  EBASCO  (1989). 

However,  when  making  confidence  statements  of  May’s  (1985),  and  of  all 
the  data,  a  different  conclusion  is  obvious.  The  confidence  statement  is  made 
for  sand  thickness  which  is  derived  from  die  kriging.  This  is  commonly  given 
as  95  percent  confidence  with  a  +  or  -  factor  which  is  double  the  standard 
deviation. 

This  results  in  a  +  or  -  factor  which  is  greater  than  the  sand  thickness  for 
the  test  case,  even  though  die  sand  thickness  standard  deviations  were  less 
than  the  standard  deviation  of  the  whole  data  set.  Table  5  shows  die  confi¬ 
dence  of  locations  shown  in  Figure  31  and  whether  die  sand  thickness  stan¬ 
dard  deviation  was  less  than  die  standard  deviation  of  die  whole  data  set 
(reliable)  or  not,  for  selected  locations  in  the  test  case,  as  well  as  for  May’s 
(1985)  data,  and  all  die  data. 

The  obvious  solution  to  die  confidence  statement  problem,  is  to  add  more 
boring  locations.  However,  since  all  die  data  available  has  already  been  used, 
with  die  sand  thickness  standard  deviations  being  greater  than  the  standard 
deviation  of  the  whole  population  in  some  areas,  there  are  not  enough  boring 
locations  to  be  able  to  do  this. 

Another  approach  had  to  be  taken  to  further  investigate  this  portion  of  die 
rationale  for  exploration.  This  approach  will  be  described  in  die  Supplemen¬ 
tary  Case  Study  section. 

Despite  the  confidence  statement  problem,  the  rationale  for  exploration  did 
follow  the  target  sand  through  the  site.  This  can  be  seen  in  Figure  32  which 
shows  the  subtle  shift  of  die  boring  locations  in  die  direction  of  die  target 

sand. 

Another  problem  was  encountered  during  the  application  of  the  rationale 
for  exploration.  This  was  a  problem  of  justifying  foe  variograms  for  foe 
kriging  with  limited  data.  Variograms  usually  contain  considerably  more  data 
than  that  of  foe  boring  locations  of  foe  hypothetical  drilling. 

A  variogram  was  obtained  with  foe  Geo-EAS  program  for  each  data  set. 
The  same  type  and  model  was  obtained  in  each  case.  These  are  contained  in 
Appendix  F.  To  justify  foe  types  and  models  as  usable  for  these  data  sets. 
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TabU  5 

Confidence  Statements  of  Selected  Locations  at  the  Rocky 
Mountain  Arsenal 


Taat  Cm*  Data 


May's  Data 


Al  Data 


95  Paroant  Confidartca  for 
Thieknaaa.  ft 


Yas 

5  ±  16 

Yaa 

5  ±  16 

Yaa 

5  ±  12 

Yaa 

5  ±  12 

Yaa 

5  ±  12 

Yaa 

5  ±  12 

Yaa 

0  ±  14 

Yaa 

7.5  ±  34 

Yaa 

8  ±  28 

Yaa 

32.5  ±  36 

Yaa 

7.5  ±  34 

Yaa 

20  ±  30 

No 

45  ±  36 

No 

>45  ±  40 

No 

0  i  32 

Yaa 

0  ±  26 

Yaa 

5  ±  22 

Yaa 

0  4  24 

Yaa 

17.5  ±  20 

No 

15  ±  24 

some  modem  streams  were  used  to  obtain  vanograms  from  larger  data  sets 
This  will  be  more  fully  described  in  die  Geostatistics  section. 
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This  section  provides  the  details  of  establishing  and  confirming  the 
environment  of  deposition  for  the  "target"  sand  body  at  the  Rocky  Mountain 
Arsenal  and  the  reason  for  thickness  variation  of  and  termination  of  that  "tar¬ 
get"  sand  body. 


Interpretations 

From  various  data  available  from  borings  at  the  Rocky  Mountain  Arsenal, 
the  environment  of  deposition  of  the  "target"  sand  was  determined.  As  stated 
in  die  previous  section,  this  data  included  grain  size  analysis.  X-ray  diffraction 
analysis,  geophysical  logs,  and  detailed  core  descriptions.  The  information 
available  for  each  boring  of  die  hypothetical  drilling  was  used  to  define  the 
sand  body.  This  information  is  shown  in  Figure  33. 

Core  descriptions  from  borings,  geophysical  logs,  and  grain  size  analyses 
were  available  for  the  four-numeral  boring  locations  with  die  exception  of 
boring  1239,  which  had  no  core  description.  Additionally,  borings  1228  and 
12S1  contained  grain  size  statistical  analyses  and  X-ray  diffraction  analyses. 
Both  of  those  boring  locations  intersected  die  "target"  sand  during  die  hypo¬ 
thetical  exploration  drilling. 

Information  for  die  three-numeral  boring  locations  was  obtained  from 
cross-section  and  tabular  data  contained  in  previously  produced  geological 
reports. 

The  cross-section  plus  tabular  data  and  die  available  core  descriptions  were 
used  to  produce  the  lithologic  characterizations  at  the  boring  site  locations. 
Those  are  shown  for  individual  borings  in  a  portion  of  Figures  15  and  18  and 
in  cross-  sections  contained  in  Figures  34,  35,  and  36.  Core  descriptions  are 
contained  in  Appendix  G.  Cross-sectional  and  tabular  information  are  con¬ 
tained  in  Appendix  H. 

Little  sedimentary  structure  was  described  in  die  cores.  This  may  have 
been  due  to  the  poorly  consolidated  condition  of  die  sand  combined  with  die 
method  of  sampling.  Split  spoons  and  pitcher  samplers  were  pushed  into  the 
sediment  for  sampling.  This  type  of  sampling  of  poorly  consolidated  saturated 
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Grain 

X-Ray 

Boring 
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Number 
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1124 
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X 

X 

X 

1155 

SP15 

X 

X 

X 

1143 
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X 

X 

X 

1160 

SP02 

X 

X 

X 

1188 

E01 

X 

X 

X 

1123 

SP08 

X 

X 

X 

1185 

N06 

X 

X 

X 

1228 

AP01 

X 

X 

X  X 

746 

- 

Tabular  and  Cross-Sectional  Information  Only 

1153 

SP09 

X 

X 

X 

756 

— 

Tabular  and  Cross-Sectional  Information  Only 

1251 

AP25 

X 

X 

X  X 
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— 

Tabular  and  Cross-Sectional  Information  Only 

757 

— 

Tabular  and  Cross-Sectional  Information  Only 

758 

- 

Tabular  and  Cross-Sectional  Information  Only 

1247 

AP21 

X 

X 

X 

1148 

SP12 

X 

X 

X 

1239 

API  2 

X 

X 

Note:  X  -  Available 


Figure  33.  Available  information  for  specific  boring  locations  at  the  Rocky 
Mountain  Arsenal 
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Elevation  in  Feel 


Figure  34.  Cross-sectional  view  down  dip  through  the  "target”  sand  body  at 
the  Rocky  Mountain  Arsenal 
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35.  Northern  cross-sectional  view  across  the  "target"  sand  body  at  the  Rocky  Mountain  Arsenal 


sands  would  destroy  any  faint  structure  that  might  be  present.  The  sand  is 
called  massive  in  some  descriptions,  however.  A  massive  sand  should  be 
expected  in  the  basal  portions  of  the  fluvial  sequence  due  to  higher  transport 
energy  in  the  deeper  portion  of  the  channel.  Although  any  sedimentary  struc¬ 
tures  that  might  have  been  present  in  the  sand  were  not  obtained,  the  lithology 
was  obtained  and  was  instrumental  in  interpreting  the  depositional  environ¬ 
ment.  This  is  exemplified  by  looking  at  the  cross-sections  in  Figures  34, 

33,  and  36.  Figure  34  shows  the  "target”  sand  body  in  the  direction  of  dip. 

It  should  be  noted  that  die  top  of  the  sand  body  is  relatively  smooth,  dropping 
in  elevation  in  the  direction  of  dip.  The  contact  of  this  sand  with  the  underly¬ 
ing  clay  is  sharp,  but  irregular.  The  bottom  of  such  sand  bodies  are  often 
irregular,  but  the  irregularity  in  the  cross-section  may  be  due  to  the  location 
of  the  boring  within  the  width  of  the  sand  body.  Boring  1251  would  be  in  the 
thickest  portion,  while  borings  1228  and  746  are  closer  to  the  edges  of  the 
sand  body,  with  boring  757  being  even  closer  to  the  edge  of  the  sand  body. 
The  cross-section  in  Figure  34  also  shows  why  the  "target"  sand  does  not  con¬ 
tinue  through  the  southern  most  boundary  of  die  site.  Quaternary  and 
Tertiary  alluvial  deposits  have  replaced  die  sand,  indicating  that  it  was  eroded 
and  removed  by  die  Quaternary  and  Tertiary  processes  which  deposited  the 
alluvium. 

The  cross-section  in  Figure  35  shows  that  the  western  side  of  the  sand  may 
also  have  been  terminated  by  the  Quaternary  and  Tertiary  alluvium. 

The  cross-section  in  Figure  36  shows  that  the  Quaternary  and  Tertiary 
alluvium  is  continuous  along  the  southern  boundary  at  the  elevation  where  the 
"target"  sand  would  be  expected. 

The  electrical  resistivity  portion  of  the  geophysical  logs  were  helpful  in 
determining  relative  grain  sizes  and  establishing  contacts  between  different 
geological  materials  as  depicted  in  Figure  37  and  a  portion  of  Figures  15 
and  18.  The  geophysical  logs  that  were  available  on  the  borings  of  the  hypo¬ 
thetical  drilling  program  are  contained  in  Appendix  I. 

Vertical  distribution  of  the  grain  size  within  the  "target"  sand  was  deter¬ 
mined  by  sieve  analyses  from  tbo  core  samples,  taken  at  various  depths  within 
the  sand.  Only  two  locations,  borings  1228  and  1251,  which  encountered  the 
"target"  sand  in  the  hypothetical  exploration,  had  sieve  analyses  in  the  "tar¬ 
get”  sand.  The  gradation  curves  for  the  analyses  for  these  two  borings  are 
contained  in  Appendix  J.  The  median  grain  size  of  each  analysis  was  plotted 
with  depth  for  each  boring  and  is  shown  in  a  portion  of  Figures  15  and  18. 

In  addition  to  the  sieve  analyses,  the  inclusive  graphic  standard  deviation 
was  used  as  a  measure  of  sorting.  Using  Folk’s  (1980)  formula: 

to  include  90°  of  the  cumulative  percent  distri- 

4  6.6 

bution  of  grain  size  in  a  sample  (4  =  -logj  mm).  The  <£84  is  the  <f>  value 
where  the  cumulative  curve  crossed  the  84  percent  line.  Alluvial  deposits  are 
moderately  to  poorly  sorted.  This  statistical  analysis  was  performed  on  the 
two  borings,  1228  and  1251,  which  encountered  the  "target"  sand  and  which 
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had  sieve  analyses.  The  determination  of  sorting  is  contained  in  Table  6.  As 
expected  for  a  fluvial  sand,  the  samples  from  boring  1228,  which  were  taken 
from  the  *****  itself  show  moderately  well  sotting  and  the  sample  from  bor¬ 
ing  1251,  which  was  taken  from  the  upper  portion  of  die  fluvial  sequence 
shows  poor  sorting. 


Table  6 

Sorting  aa  Determined  from  Inclusive  Standard  Deviation 

Boring 

Depth.  It 

Sample 

Number 

Standard 

Deviation 

Sorting 

1228 

14.0-16.3 

4 

0.89 

Moderately  Well 

1228 

18.0-20.0 

5 

0.73 

Poorly 

1261 

40.0-41.4 

15 

1.36 

Poorly 

This  tendency  to  change  toward  die  upper  portion  of  the  sequence  is  sup¬ 
ported  by  die  increase  in  matrix  material  and  decrease  in  quartz  in  the  upper 
portions  of  the  sequence.  X-ray  diffraction  analysis  was  performed  on  sam¬ 
ples  from  different  intervals  of  the  two  borings,  1228  and  1251,  which 
encountered  die  "target”  sand.  The  results  are  summarized  in  Table  7  and  are 
shown  in  a  portion  of  Figures  15  and  18.  Boring  1251  shows  a  decrease  in 
percentage  of  quartz  at  die  upper  portion  of  the  sequence.  Bodi  borings  show 
around  40  percent  quartz  through  die  middle  and  lower  portions  of  the  sand. 


Table  7 

X-Ray  Diffraction  Analysis  Results 

loring  Numbf 

Depth,  ft 

Percent  Quartz 

—  1 

Wm 

1228 

14.0  - 15.3 

38 

15 

78 

1228 

18.0  -  20.0 

38 

15 

78 

1251 

40.9-41.4 

15 

10 

— 

1251 

61.5-52.0 

42 

17 

- 

1251 

63.2  •  63.7 

37 

16 

- 

1261 

74.3  -  74.8 

39 

16 

81 

1251 

78.3  -  78.8 

42 

17 

- 

1251 

80.0  -  80.5 

35 

16 

- 

The  information  as  just  described,  when  available,  was  compiled  for  each 
boring  in  die  hypothetical  drilling  that  encountered  die  "target"  sand  as  shown 
in  Figures  15  and  18  for  borings  1228  and  1251  respect!  Tiis 


Chapter  4  Stratigraphic  Interpretatione 


63 


information  was  then  compared  to  that  for  various  depositions!  environments 
which  can  form  sand  bodies,  shown  in  Figure  2.  From  the  comparison,  the 
diagnostic  characteristics  of  the  fluvial  channel  environment  were  met  by  the 
“target"  sand  with  the  exception  of  the  sedimentary  structures,  which  were  not 
obtained.  Figure  38  shows  the  diagnostic  characteristics  for  a  fluvial  channel 
sand. 

The  sand  penetrated  by  boring  1228  is  interpreted  to  be  a  fluvial  sand.  A 
number  of  factors  confirm  this  interpretation.  The  fact  that  no  marine  or  mar¬ 
ginal  marine  fossils  were  observed  in  the  sand  sequence  or  in  strata  above  or 
below  the  sand  suggest  dud  deposition  occurred  inland,  away  from  direct 
marine  influence.  The  absence  of  glauconite  is  another  factor  pointing  toward 
a  fluvial  origin  for  the  sand.  Finally,  the  stratigraphic  position  of  die  sand 
sequence,  in  reference  to  the  large  volume  of  sediments  within  the  Denver 
basin,  places  it  in  an  area  where  fluvial  deposits  have  prograded  out  and  over 
deltaic  and  marine  sediments. 
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5  Supplementary  Case  Study 


As  discussed  in  the  Case  Study  section  there  were  a  few  problems 
encountered  in  applying  the  rationale  for  exploration  to  the  case  study  at  the 
Rocky  Mountain  Arsenal.  This  section  addresses  one  of  those  problems,  that 
of  insufficient  data  locations  of  die  depth  needed  to  encounter  foe  Arsenal 
"target”  sand. 

To  alleviate  this  problem,  another  site  was  selected.  This  "supplementary" 
she  was  chosen  from  a  modern  day  stream.  A  modern  stream  allowed  for  a 
dose  to  absolute  definition  of  a  sand  body,  for  comparison  of  what  is  visibly 
defined  with  that  which  is  defined  by  foe  rationale  for  exploration. 


Site  Description 

The  site  chosen  for  this  supplemental  application  of  foe  rationale  for 
exploration  is  a  portion  of  foe  Brazos  River  floodplain  in  central  Texas.  The 
general  location  is  shown  in  Figure  39.  Availability  of  satellite  imagery  and  a 
corresponding  topographic  map  dictated  foe  choice  of  foe  Brazos  River. 

Satellite  imagery  and  a  U.S.  Geological  Survey  topographic  map  of  foe 
same  area,  foe  Austin,  Texas  1:250,000  scale  map,  were  used  to  obtain 
approximate  floodplain  boundaries.  These  floodplain  boundaries  subsequently 
were  assumed  to  represent  foe  width  of  foe  stream’s  sand  body.  The  site  is  a 
275,000  ft  east  to  west  and  225,000  ft  north  to  south  area,  bounded  by  foe 
central  Texas  state  grid  coordinate  system  lines  of  3,100,000  ft  and  3,275,- 
000  ft  East,  and  275,000  ft  and  500,000  ft  North.  It  encompasses  foe  Brazos 
River  from  just  west  of  Heame,  Texas  to  near  where  foe  Navasota  River  joins 
foe  Brazos  River  at  Navasota,  Texas.  Figure  40  shows  foe  she  boundaries 
and  foe  Brazos  River  floodplain  boundaries. 

From  measuring  foe  width  of  foe  Brazos  River  floodplain  boundaries,  foe 
sand  body  width  was  mapped  to  be  approximately  24,000  ft.  From  this 
width,  a  sand  body  thickness  of  approximately  80  ft  was  back  calculated  from 
Lorenz’s  and  Leader’s  equations  (described  in  Rationale  for  Exploration  sec¬ 
tion).  This  80  ft  thickness  was  assumed  to  be  uniform  for  any  location  within 
foe  sand  body.  Thus  both  a  thickness  and  boundary  for  foe  hypothetical  sand 
body  were  obtainable  for  any  location  in  foe  she. 
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Figure  40.  Brazos  River  site  and  floodplain  boundaries 
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Application  of  the  Rationale  for  Exploration 


With  die  assumptions  described  above,  the  rationale  for  exploration  was 
applied  to  the  site.  Thus,  any  location  selected  by  this  rationale  provided  the 
needed  information,  and  the  process  could  continue  for  any  level  of  infor¬ 
mation  desired. 

The  boundary  of  the  site  was  hypothetically  drilled  with  25,000  ft  spacing 
between  boring  locations.  The  spacing  was  chosen  somewhat  arbitrarily, 
although  this  spacing  could  have  been  hypothetically  based  on  encountering 
sand  body  width(s)  of  concern,  such  as  that  of  the  Brazos  River.  These  are 
shown  in  Figure  41.  Figure  41  also  shows  which  locations  encountered  the 
assumed  80  ft  thick  sand  body.  Even  though  the  sand  would  have  been 
encountered  early  in  the  hypothetical  drilling,  the  boundary  was  completely 
drilled,  as  would  be  the  case  in  most  environmental  applications. 

The  hypothetical  sand  body  was  then  followed  through  the  site  from  the 
southern  most  boring  location  which  intersected  the  sand  body.  Again,  envi¬ 
ronmentally,  this  would  be  an  exit  point  for  any  pollutant,  so  the  down  gradi¬ 
ent  and  thus  potential  exit  point  would  need  to  be  defined  first. 

Since  the  boundary  had  been  hypothetically  drilled,  several  data  points 
existed.  These  were  used  to  produce  contours  of  the  standard  deviation  of  the 
thickness  estimates  obtained  from  kriging  with  the  Geo-EAS  program,  in  the 
same  manor  as  described  previously  in  the  Case  Study  section.  These  stan¬ 
dard  deviations  are  shown  in  Figure  42. 

From  the  curves  in  Figure  5,  the  interval  for  data  point,  or  hypothetical 
boring,  spacing  of  approximately  17,500  ft  was  obtained.  This  is  the  spacing 
for  a  50  percent  probability  of  dual  penetration  of  a  24,000  ft  wide  sand  body 
in  order  to  define  the  sand  body. 

This  distance,  of  24,000  ft  from  the  boring  which  encountered  the  hypo¬ 
thetical  sand,  was  in  a  area  where  the  sand  thickness  standard  deviation  was 
less  than  the  standard  deviation  of  the  whole  data  set.  One  boring  was  hypo¬ 
thetically  placed  where  the  highest  standard  deviation  occurred  at  that  dis¬ 
tance,  as  shown  in  Figure  43. 

Often  during  the  hypothetical  exploration,  the  boring  spacing  would  not  be 
enough  to  place  an  additional  boring  out  of  the  area  where  the  sand  thickness 
standard  deviations  were  less  titan  the  standard  deviation  of  the  whole  data  set. 
When  this  happened,  usually  only  one  additional  boring  was  placed  where  the 
higher  standard  deviation  would  occur.  Only  in  a  couple  of  instances  were 
two  borings  chosen  at  the  same  time  due  to  a  large  enough  area  (regarding  the 
boring  spacing)  with  the  same  standard  deviation. 

This  procedure  was  followed  during  the  hypothetical  exploration  through 
tiie  she,  defining  the  sand.  The  fifteenth  boring  that  was  added  after  the 
boundary  borings,  produced  estimates  of  the  sand  thicknesses  whose  standard 
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Figure  41 .  Brazos  River  boundary  for  bonng  locations 
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Figure  42.  Standard  deviation,  in  feet,  for  the  Brazos  River  boundary  data 
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deviations  were  less  than  that  of  the  whole  data  set.  This  was  determined  by 
the  contoured  thickness  estimates,  compared  to  die  contoured  standard  devia¬ 
tion  of  those  thicknesses.  This  was  done  for  the  area  in  which  the  estimated 
sand  body  thicknesses  occurred. 

The  proximity  of  the  first  of  the  definition  borings  to  the  south  and  east 
site  boundary  borings  coincidentally  guided  the  first  boring  locations  in  the 
direction  of  the  sand,  and  provided  low  enough  standard  deviations  so  that  no 
borings  wore  needed  outside  the  sand  body.  As  the  borings  moved  further 
into  die  site,  the  procedures  resulted  in  the  highest  sand  thickness  standard 
deviation  and  occasionally  the  standard  deviation  which  was  greater  than  the 
standard  deviation  of  die  whole  data  set  to  occur  in  the  direction  of  the  sand 
body.  When  the  procedures  took  the  location  of  a  boring  out  of  the  sand 
body,  die  next  location  would  be  chosen  in  a  direction  that  resulted  in  inter¬ 
secting  the  sand  body  again,  always  with  the  17,500  ft  spacing  and  in  the 
direction  of  die  highest  standard  deviation. 

On  two  occasions,  the  results  of  a  boring  location  placement  resulted  in 
sand  thickness  estimates  whose  standard  deviations  were  less  than  the  standard 
deviation  of  the  whole  data  set  and  wore  uniform  surrounding  it.  In  these 
instances  the  next  boring  location  was  chosen  in  the  direction  of  the  northern 
boundary  boring  which  encountered  the  sand.  Had  there  not  been  another 
boundary  boring  that  had  encountered  die  sand  body,  three  borings  with 
equal,  17,500  ft,  spacing  from  one  another  would  have  had  to  have  been 
drilled  to  continue  the  exploration. 

The  definition  of  the  sand  body  was  considered  complete  once  the  standard 
deviations  of  the  sand  thicknesses  were  less  than  die  standard  deviation  of  the 
whole  data  set.  This  occurred  after  the  fifteenth  boring.  However,  four  more 
boring  locations  were  chosen  at  locations  with  relatively  higher  standard 
deviations  in  order  to  achieve  relatively  uniform  standard  deviation  for  die 
area  containing  the  sand  body.  All  nineteen  boring  locations  are  shown  in 
Figure  44. 

Appendix  K  contains  the  contoured  standard  deviations  and  each  data  point 
location  as  it  was  selected,  in  die  order  in  which  die  hypothetical  exploration 
proceeded.  The  standard  deviation  for  each  of  the  data  sets  is  contained  in 
Table  8.  Each  data  set  contains  all  die  previous  borings  plus  the  one(s)  added 
at  each  step,  as  described  in  the  Rock  Mountain  Arsenal  test  case. 


Application  of  the  Grid  Method  for  Exploration 

The  hypothetical  sand  body  was  also  defined  by  grid  drilling  for  compari¬ 
son  purposes.  The  spacing  was  kept  the  same  as  the  original  boundary  bor¬ 
ings  (25,000  ft).  The  locations  of  these  borings  are  shown  in  Figure  45.  A 
total  of  thirty-four  borings,  in  addition  to  the  boundary  borings,  would  have 
been  drilled  in  the  site  containing  the  sand  body  using  the  typical  grid  method 
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Figure  44.  Boring  locations  added  to  define  the  hypothetical  sand  body  by 
the  rationale  for  exploration  method,  included  are  the  boundary 
borings 
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Table  8 

Data  Sat  Population  Standard  Deviations  for  the 

Braxos  River 

Boring  Papulation 

Varianoa 

Standard  Deviation 

Boundary  Boringa 

376.00 

19.3 

Sat  1 

528.93 

23.0 

Sat  2 

664.36 

25.8 

Sat  3 

783.67 

28.0 

Sat  4 

888.89 

29.8 

Sat  5 

981.74 

31.3 

Sat  6 

1.043.50 

32.3 

Sat  7 

1.024.00 

32.0 

Sat  8 

1,096.50 

33.1 

Sat  9 

1.161.00 

34.1 

Sat  10 

1,218.40 

34.9 

Sat  11 

1,200.00 

34.6 

Sat  12 

1,251.60 

35.4 

Sat  13 

1,297.50 

36.0 

Sat  14 

1,338.50 

36.6 

SatlS 

1,306.10 

36.1 

Final  Sat 

1,328.70 

36.4 

for  exploration.  This  was  twenty-five  more  borings  than  required  by  die 
rational  for  exploration. 

Had  geology  been  taken  into  account  in  this  typical  grid  method  during  the 
exploration,  and  limited  die  number  of  borings  to  those  just  outside  of  the 
sand  body,  as  shown  in  Figure  46,  twenty-five  borings,  in  addition  to  the 
boundary  borings,  would  have  been  necessary  to  define  the  sand  body.  This 
is  still  ten  more  than  required  by  the  rational  for  exploration,  to  define  the 
sand  body  reliably. 

This  data  set  was  kriged  to  obtain  estimated  thicknesses  and  standard  devia¬ 
tions  of  those  estimates  for  comparison  with  those  obtained  by  following  the 
rationale  for  exploration.  This  would  not  normally  be  done  in  a  gridding  type 
of  exploration  program. 


Chapter  5  Supplementary  Cm*  Study 


75 


Comparison  of  Methods 


Several  parameters  from  the  gridding  method  for  exploration  were  com¬ 
pared  to  their  counterparts  obtained  from  the  rationale  for  exploration  method. 
These  included  the  appearance  of  the  estimated  thicknesses  produced  by 
kriging,  the  number  and  spacing  of  boring  locations  needed  by  each  to  define 
the  sand  body,  and  the  confidence  statements  for  selected  locations  within  the 
sand  body. 

The  rationale  for  exploration  method  predicted  the  sand  body  with  uniform 
sand  thickness  standard  deviations  from  nineteen  borings  in  addition  to  die 
boundary  drilling.  The  gridding  method  for  exploration  required  twenty-five 
borings  to  define  die  sand  body.  The  gridding  method  also  provided  results  in 
which  die  sand  thickness  standard  deviations  were  less  than  the  standard 
deviation  of  die  whole  data  set,  but  kriging  is  not  normally  done  for  this 
method.  The  rationale  for  exploration  boring  spacing  was  17,500  ft  while  the 
grid  method  boring  spacing  was  larger,  at  25,000  ft.  Even  more  boring  loca¬ 
tions  would  have  been  necessary  for  the  grid  method  to  define  the  sand  body 
with  a  spacing  of  17,500  ft. 

The  sand  thickness  estimates  of  both  methods  are  very  similar,  as  shown  in 
Figures  47  and  48.  The  standard  deviation  of  the  estimated  thicknesses  of  die 
grid  method  are  slightly  lower  than  those  from  the  rationale  for  exploration. 
The  sufeodard  deviations  for  each  method  are  shown  in  Figures  49  and  50. 

The  slighdy  lower  standard  deviations  for  the  grid  method  leads  to  slightly 
better  confidence  statements.  However,  the  difference  in  directly  comparable 
locations  (die  same  thicknesses)  is  only  a  few  feet.  The  locations  for  these 
confidence  statements  are  shown  in  Figures  51  and  52.  The  confidence  state¬ 
ments  are  contained  in  Table  9. 

Estimated  thicknesses  were  also  obtained  from  data  of  the  entire  site  on  a 
10,000  ft  grid  spacing,  which  was  used  to  produce  a  variogram  for  die  site. 
The  locations  of  these  data  points  are  shown  in  Figure  53.  Discussion  relating 
to  the  variogram  will  be  contained  in  die  Geostatistics  section.  The  standard 
deviation  for  this  grid  was  uniform,  thus,  no  contours  were  possible.  How¬ 
ever,  the  estimated  thickness  from  this  data  set  closely  approximated  the  hypo¬ 
thetical  sand  body’s  geometry.  These  estimated  thicknesses  are  shown  in 
Figure  54. 

The  sand  body’s  appearance  from  the  estimated  thicknesses  from  both 
methods  of  exploration  are  noticeably  wider  than  that  of  this  large  data  set. 
This  is  due  to  the  differences  in  data  point  spacing,  as  well  as  the  number  of 
data  points  used.  Obtaining  such  a  large  number  of  data  points  obviously 
adds  to  die  definition  and  confidence,  but  is  not  reasonable,  for  reasons 
described  in  the  Introduction  section. 

This  larf  *  data  set  was  useful  in  looking  at  the  confidence  statement  loca¬ 
tions.  Confidence  statements  from  locations  corresponding  to  die  edge  of  the 
sand  body,  which  is  shown  in  Figure  55,  are  shown  for  both  methods  in 
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Figure  48.  Estimated  sand  thickness,  in  feet,  from  a  typical  grid  exploration 
method 
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Figure  50.  Standard  deviation,  in  feet,  for  estimated  sand  thickness  from  a 
typical  grid  method  of  exploration 
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Figure  52.  Locations  of  confidence  statements  in  Table  9  for  a  typical  grid 
exploration  method 
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Table  9 

Confidence  Statements  for  Brazos  River 
Figures  51  and  52 

Locations  Shown  in 

Location 

Orid  Method 

Rational*  for  Exploration  Method 

1 

60  ±  <16 

40  ±  <18 

2 

10  ±  16 

25  ±  32 

3 

IS  ±  18 

<10  ±  >36 

4 

70  ±  <16 

50  ±  12 

I  5 

10  t  16 

45  t  27  1 

6 

10  ±  22 

10  ±  28  1 

7 

80  ±  <16 

55  ±  18  H 

8 

20  ±  16 

10  t  18  I 

9 

10  ±  20 

20  ±  30  | 

10 

70  ±  <16 

60  ±  19 _ I 

Table  10.  Figures  56  and  57  show  die  locations  where  the  confidence  state¬ 
ments  were  obtained  with  respect  to  each  method’s  estimated  sand  body  thick¬ 
nesses.  This  comparison  shows  that  in  the  area  where  the  close  approximation 
shows  10  ft  thicknesses,  both  cases  show  at  least  20  ft  thicknesses  and  greater. 
Assuming  that  the  error  is  negative  in  each  confidence  statement  gives  esti¬ 
mates  that  are  relatively  closer  to  that  of  the  close  approximation,  with  some 
exceptions.  The  exceptions  were  the  rationale  for  exploration  confidence 
thicknesses  which  were  noticeably  higher.  This  results  from  the  comparison 
locations  picked.  By  viewing  die  estimated  thicknesses  for  each  method, 
changing  the  locations  for  confidence  statements  could  cause  the  grid  method 
to  have  noticeably  high  estimates  of  thickness  and  the  rationale  for  exploration 
estimates  of  thickness  to  all  be  relatively  close. 

In  both  methods,  the  thicknesses  are  noticeably  over-estimated  near  the 
close  approximation  sand  boundary.  However,  both  methods  indicate  that  at 
those  boundary  locations  the  sand  body  was  present.  In  other  words,  the 
thickness  estimate  is  larger  than  the  possible  error.  Some  of  the  locations 
from  Figures  51  and  52  show  a  potentially  negative  error  that  if  subtracted 
from  die  estimated  thickness  would  result  in  a  negative  estimate.  In  these 
cases  it  cannot  be  said  for  sure  that  the  sand  body  exists  at  that  location. 

From  this  data  it  appears  that  if  the  negative  error  is  taken  and  subtracted 
from  the  estimated  thickness,  the  sand  body  boundary  is  approximated  by 
resulting  numbers  approaching  0,  but  not  becoming  negative. 

As  stated  earlier,  better  confidence  can  be  established  by  increasing  data 
point  locations.  However,  the  objective  is  to  define  the  sand  body  with  the 
least  number  of  data  points.  The  traditional  grid  method  yielded  similar 
results,  in  estimated  sand  thickness  and  standard  deviations  of  those  estimates, 
to  those  produced  by  the  rationale  for  exploration.  With  the  rationale  for 
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Figure  53.  Data  point  locations  for  the  Brazos  River  10,000  foot  grid 
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Figure  54.  Estimated  sand  thickness,  in  feet,  from  the  Brazos  River 
10,000  foot  grid 
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Figure  55.  Locations  chosen  for  the  confidence  statements  in  Table  10 
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Table  10 

Confidence  Statements  (in  ft)  for  Brazos  River  Locations  Shown  in 
Figures  55,  56,  and  57 

|  Explanation 

Location  | 

i 

2 

3 

4 

6 

e 

I  Grid  Method 

1 95  paroant 
Confidence 

30  ±  <16 

20  ±  <16 

25  ±  16 

20  ±  <16 

30  ±  16 

35  ±  <16 

Error  Subtracted 

<14 

<4 

9 

<4 

14 

<19  § 
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exploration  requiring  fewer  data  point  locations,  then  the  definition  of  the  sand 
body  with  the  least  number  of  data  points  was  accomplished  by  the  rationale 
for  exploration. 

Also  important  is  that  the  sand  body  could  be  followed  through  the  site  as 
in  the  Case  Study  section.  This  is  particularly  useful  for  some  applications, 
such  as  oil  and  gas  exploration.  Table  1 1  contains  a  summary  of  the  compari¬ 
son  of  the  definition  of  the  sand  body  by  the  traditional  grid  method  and  by 
foe  method  described  by  die  rationale  for  exploration. 
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Figure  56.  Locations  of  confidence  statements  in  Table  10  for  the  rationale 
for  exploration  method 
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Table  11 

Summary  of  Comparison  for  Soloctad  Parameters  of  Different 
Exploration  Methods 


Grid  Method  I  Rational*  (or  Exploration  Mothod 


Odd  epaojng  26,000  ft 

Number  of  dofinition  boring* _ 25 _ 

Thick  no**  aatimato*  comparod  Higher 

to  thoaa  of  10,000  ft  grid 


17,500  ft  (boa ad  on  «*nd  body  width) 

IS  (roiiobiy) _ 

Higher  (attghtly  higher  than  thoaa  of 
tha  grid  method) 
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6  Geostatistics 


A  statistical  method  was  chosen  to  produce  estimates  of  sand  thickness  and 
subsequently  errors  for  those  thicknesses.  The  method  chosen,  which  has 
already  been  named  in  earlier  sections,  was  kriging.  Specifically,  ordinary 
block  kriging  was  used. 

Kriging  is  a  form  of  weighted  local  averaging.  Kriging  is  considered 
optimal  for  geologic  data  sets  by  numerous  authors  (Davis  1986,  David  1977, 
Clark  1979,  etc.).  This  is  because  die  method  provides  estimates  of  values  at 
unrecorded  places  without  bias  and  with  minimum  and  known  variance,  pro¬ 
vided  there  is  a  model  for  die  semi-variogram.  Kriging  produces  those  esti¬ 
mates  with  a  lower  number  of  observations  than  that  of  conventional  methods. 
Clark  (1979)  gave  the  points  of  major  importance  found  in  numerous  publica¬ 
tions  as:  (1)  Given  the  basic  assumptions,  no  trend,  and  a  model  for  die 
semi-variogram,  kriging  always  produces  the  best  linear  unbiased  estimator. 

(2)  If  the  proper  models  are  used  for  die  semi-variogram,  and  die  system  is 
set  up  correcdy,  there  is  always  a  unique  solution  to  die  kriging  system. 

(3)  If  you  try  to  estimate  the  value  at  a  location  which  has  been  sampled,  the 
kriging  system  will  return  the  sample  value  as  the  estimator,  and  a  kriging 
variance  of  zero.  In  other  words,  you  already  know  that  value.  This  is 
usually  referred  to  as  an  exact  interpolator.  (4)  If  you  have  regular  sampling, 
and  hence  the  same  sampling/block  setup  at  many  different  positions  within 
die  deposit,  it  is  not  necessary  to  recalculate  the  kriging  system  each  time. 
Figure  58  shows  results  of  Di  et  al.’s  (1989)  comparison  of  a  conventional 
method  to  die  kriging  method  and  that  the  kriging  method  produced  lower 
standard  errors.  Conventional  methods  include  standard  errors  of  die  mean, 
student’s  t  test,  least  squares  analyses,  analysis  of  variance,  etc.  The  above 
descriptions  served  as  die  basis  for  selecting  kriging  for  this  application. 

The  semi-variogram,  also  termed  simply  the  variogram,  is  a  curve  of  the 
variation  of  a  numeric  variable,  thickness  in  this  application,  versus  distance 
between  pairs  of  known,  or  control,  points.  Estimates  of  the  semi-variogram 
are  used  to  determine  the  weights  applied  to  the  data  when  computing  die 
averages  and  are  presented  in  the  kriging  equations. 

Since  kriging  is  an  advanced  technique  involving  intense  computing,  it 
must  be  done  with  a  computer.  For  this  application,  as  mentioned  in  a  previ¬ 
ous  section,  a  public  domain  program  was  used.  This  program  which  is 
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STANDARD  ERROR 


Figure  58.  Di  and  others'  comparison  of  conventional  method  to  the  kriging 
method  (after  Di  et  al.  1 989) 
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called  Geo-EAS  (Geochemical  Enviionmental  Assessment  Software),  was 
developed  by  the  U.S.  Environmental  Protection  Agency.  This  program  con¬ 
tains  sub-routines  which  allow  for  die  generation  of  the  variogram,  kriging  die 
variable  values  with  die  generated  variogram,  and  contouring  die  output  val¬ 
ues.  The  kriging  also  produces  standard  deviation  for  the  kriged  values, 
which  can  be  contoured. 

Because  the  Geo-EAS  program  is  public  domain  and  contains  all  die 
desired  routines  in  one  package,  it  was  used  for  the  calculations  of  die  kriging 
method  and  subsequent  contouring.  Also,  if  used  in  field  application,  speed 
of  obtaining  the  needed  locations  is  crucial,  therefore  a  relatively  easy  to  use 
program  with  all  die  needed  products  is  desirable. 


Obtaining  the  Variogram 

There  is  disagreement  in  the  geostatistical  literature  in  using  semi- 
variogram  or  variogram  as  die  correct  term.  For  this  application,  variogram 
will  be  used,  for  simplicity. 

The  variogram  as  defined  by  Englund  and  Sparks  (1988)  is  a  plot  of  the 
variance  of  paired  sample  measurements  as  a  function  of  the  distance  between 
samples.  Variograms  provide  a  means  of  quantifying  die  commonly  observed 
relationship  that  samples  close  together  will  tend  to  have  more  similar  values 
than  samples  far  apart. 

The  variogram  is  necessary  for  kriging,  and  is  a  critical  part  of  this,  as 
well  as  any  geostatistical,  study.  The  variogram  is  the  interpretation  of  die 
spatial  correlation  structure  of  die  sample  data  set.  It  controls  die  way  that 
kriging  weights  are  assigned  to  samples  during  interpolation,  and  consequently 
controls  die  quality  of  the  results. 

Englund  and  Sparks  (1988)  point  out  that  all  interpolation  and  contouring 
methods  make  the  assumption  that  some  type  of  spatial  correlation  is  present, 
that  is,  they  assume  that  a  measurement  at  any  point  represents  nearby  loca¬ 
tions  better  than  locations  farther  away.  Variogram  analysis  attempts  to 
quantify  this  relationship.  In  other  words,  how  well  can  a  measurement  be 
expected  to  represent  another  location  a  specific  distance  away?  Experimental 
variograms  plot  die  average  difference  of  pairs  of  measurements,  as  one  half 
the  squared  difference  (variance),  against  die  distances  separating  the  pairs.  If 
all  possible  sample  locations  were  measured,  a  true  variogram  could  be  com¬ 
puted  for  the  site  where  die  variance  of  all  pairs  of  measurements  would 
satisfy  each  combination  of  distance  and  direction.  Since,  this  is  not  usually 
possible,  limited  data  is  used  to  compute  variances  and  then  plot  a  graph  of 
the  variances  verms  distance.  Then,  a  curve  is  fitted  to  die  graph.  This 
model  is  assumed  to  be  an  approximation  of  die  true  variogram. 

Several  types  of  variograms  are  possible.  The  Geo-EAS  program  allowed 
choices,  so  that  each  type  could  be  compared.  These  were  the  ordinary 
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variogram,  relative  vartogram,  "mad’ogram",  and  the  non-ergodic  or  Covario- 
gram.  Engiund  and  Sparks  (1988)  describes  each  as  follows.  The  relative 
variogram  is  analogous  to  die  relative  standard  deviation  often  used  to  mea¬ 
sure  analytical  variability.  When  modeled  and  used  for  kriging  the  relative 
kriging  standard  deviations  must  be  multiplied  by  the  estimated  values  to  be 
comparable  with  kriging  standard  deviations  produced  with  ordinary  vario¬ 
gram  models.  The  "mad’ogram"  plots  the  mean  absolute  differences,  but  is 
not  recommended  for  kriging.  The  non-ergodic  or  covariogram  is  based  on 
estimates  of  covariance  rather  than  variance.  The  covariograms  have  the  same 
units  as  ordinary  variograms  and  may  be  modeled  and  used  for  kriging  in  the 
same  way.  The  covariance  values,  rather  than  variogram  values,  are  actually 
used  in  the  Geo-EAS  kriging  matrix  equations  for  greater  computational 
efficiency. 

Once  the  type  of  variogram  is  selected,  there  are  several  mathematical 
models  which  may  define  the  graph.  Again,  Geo-EAS  allowed  choices,  so 
that  each  model  could  be  compared.  These  are  the  spherical,  exponential, 
linear,  and  Gaussian.  The  spherical  model  of  the  variogram  is  observed 
frequently  in  experimental  data  (Engiund  and  Sparks  1988). 

To  fit  any  model  to  the  variogram,  an  estimate  of  the  Y-intercept,  termed 
the  nugget,  is  needed.  The  difference  between  the  nugget  and  the  maximum 
Y  value,  termed  the  sill,  is  also  needed.  Finally,  the  distance  at  which  the 
model  reaches  die  maximum  value,  termed  the  range,  is  needed.  Although 
some  form  of  least  squares  criteria  could  be  used,  the  Geo-EAS  program 
selection  is  subjective,  picking  the  model  by  trial  and  error  which  gives  the 
best  fit. 

The  Geo-EAS  program  also  allows  for  using  data  in  a  specified  direction, 
or  from  all  directions,  to  specify  pair  orientation  criteria  for  the  variogram 
computations. 

As  mentioned  in  the  Case  Study  section,  there  was  a  problem  in  obtaining 
die  variogram  because  of  die  limited  amount  of  data  used  in  those  data  sets. 
Although  the  covariogram  with  a  spherical  model  in  a  specific  direction 
appeared  to  be  die  best,  it  was  questionable.  In  particular,  the  range  was 
poorly  defined.  However,  when  kriging,  the  Geo-EAS  program  allowed  for  a 
minimum  and  maximum  range  so  that  an  area  could  be  bracketed  for  the 
range.  Examples  are  the  variograms  produced  from  the  Rocky  Mountain 
Arsenal  data  set  1,  shown  in  Figure  59  and  from  the  Brazos  River  boundary 
data  set,  shown  in  Figure  60.  The  other  variograms  for  the  remainder  of  the 
data  sets  are  contained  in  Appendix  F. 

To  justify  die  type  and  model  of  variogram  used  for  this  application,  and  if 
it  should  be  directional,  larger  data  sets  for  die  same  variable,  thickness,  were 
used  to  generate  more  graphs.  One  data  set  contains  all  die  data  from  die 
Rocky  Mountain  Arsenal  that  encountered  the  target  sand,  one  contains  all  the 
data  for  a  10,000  ft  grid  for  die  Brazos  River  and  another  contains  data  from 
a  1000  ft  grid  for  the  Salt  River  in  Phoenix,  Arizona.  The  Salt  River  was 
added  to  die  sites  already  used  to  reduce  the  chance  of  a  coincidental 
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Figure  60.  Variogram  for  Brazos  River  boundary  data  set  (each  plot  may  represent  more  than  one  data  point 


agreement  between  the  type  and  model  of  variogram,  and  for  additional 
variety  of  sand  body  size  and  climatic  conditions.  The  sand  width  and  thick¬ 
nesses  for  the  Salt  River  were  obtained  in  the  same  manor  as  that  of  the 
Brazos  River,  described  in  the  Supplemental  Case  Study  section.  Satellite 
imagery  and  the  Phoenix,  Arizona  7-1/2  min  quadrangle  were  used  in  this 
case.  The  Salt  River  site  is  bounded  by  the  Arizona  state  grid  coordinates  of 
435,000  ft  and  775,000  ft  East  and  875,000  ft  and  884,000  ft  North.  The 
general  location  of  the  Salt  River  site  is  shown  in  Figure  61.  Figures  26, 

53,  and  62  show  the  data  point  locations  for  each  of  these.  Tabulation  of  the 
data  is  contained  in  Appendix  E. 

These  data  sites  allowed  for  a  variety  in  orientation  through  the  site,  size 
(width  and  thickness),  climates,  and  ages  of  sand  bodies.  In  each  case,  a 
directional  covariogram  with  a  spherical  model  was  the  best  fit  for  the  data. 
Therefore,  die  type  and  model  of  the  variogram,  which  results  from  the  data 
sets  in  the  case  studies,  are  justifiable.  Figures  63,  64,  and  65  show  these 
variograms.  Appendix  F  contains  other  variograms  used  in  this  project. 


Kriging 


As  previously  described,  kriging  is  a  weighted-moving-average  interpo¬ 
lation  method  where  the  set  of  weights  assigned  to  samples  minimizes  the 
estimation  variance,  which  is  computed  as  a  function  of  the  variogram  model 
and  locations  of  the  samples  relative  to  each  other,  and  to  the  point  or  block 
being  estimated  (Englund  and  Sparks  1988). 

As  the  above  definition  reveals,  kriging  estimates  can  be  for  an  area, 
termed  block,  or  for  a  point.  Point  kriging  usually  provides  estimates  similar 
to  block  kriging.  However,  if  a  point  being  estimated  coincides  with  a  sam¬ 
pled  location,  the  estimate  is  set  equal  to  the  sample  value.  This  is  not 
appropriate  for  contour  mapping,  which  implicidy  requires  a  spatial  estimator 
(Englund  and  Sparks  1988). 

The  kriging  estimates  can  be  produced  with  either  ordinary  or  simple 
kriging.  Ordinary  kriging  estimates  the  point  or  block  values  with  a  weighted 
average  of  die  sample  values  within  a  local  search  neighborhood,  centered  on 
die  point  or  block.  Simple  kriging  also  assigns  a  weight  to  the  population 
mean,  but  makes  a  strong  assumption  that  the  mean  value  is  constant  over  the 
site.  It  also  requires  that  the  available  data  be  adequate  to  provide  a  good 
estimate  of  die  mean  (Englund  and  Sparks  1988). 

The  kriging  portion  of  die  Geo-EAS  program  allows  for  a  selection  from 
these  choices.  From  the  above  statements,  the  ordinary  and  block  kriging 
were  chosen  for  this  application.  These  also  happen  to  be  the  default  settings 
for  the  program.  Default  settings  in  the  Geo-EAS  program  were  used  in  most 
cases  for  sake  of  simplicity  and  speed,  but  mainly  because  they  were  the  best 
choices  for  this  application,  as  evidenced  above. 
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Figure  62.  Floodplain  boundaries  and  data  point  locations  for  the  Salt  River 
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Figure  63.  Variogram  for  all  the  data  at  the  Rocky  Mountain  Arsenal  (each  plot  may  represent  more  than  one  data  point) 
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Figure  65.  Veriogram  for  the  Salt  River  data  (each  plot  may  represent  more  than  one  data  point) 


The  kriging  portion  of  the  Geo-EAS  program  produces  estimates  of  the 
variable  and  the  standard  deviation  of  those  estimates  in  a  grid  spacing  which 
is  contoured.  Thickness  was  die  variable  for  this  application.  The  thickness 
was  contoured  in  order  to  compare  the  results  of  the  computer  program  to  that 
produced  by  other  methods,  as  at  Rocky  Mountain  Arsenal,  or  that  which  was 
known,  as  that  for  the  Brazos  and  Salt  Rivers.  The  Brazos  and  Salt  River 
sites’  actual  boundaries  are  shown  in  Figures  40  and  62.  The  boundaries 
shown  by  die  thickness  estimates  can  be  seen  in  Figures  54  and  66.  The  com¬ 
parison  for  die  sites  with  known  widths  shows  good  reproduction  of  the  sand 
bodies  by  the  estimated  thicknesses.  The  standard  deviations  of  die  thickness 
estimates  were  used  in  part  to  chose  locations  for  additional  information  and 
in  part  to  establish  confidence  of  the  estimated  thickness.  These  have  already 
been  discussed  previously  in  the  Case  Study  and  Supplemental  Case  Study 
sections. 

The  level  of  discussion  presented  in  this  section  provides  an  overview  of 
die  geostatistics  used  in  this  application.  A  tutorial  is  provided  in  Appendix  A 
for  a  more  in  depth  understanding.  Procedures  can  also  be  found  in  more 
depth  from  several  source  texts  on  the  subject,  such  as  in  Davis  (1986),  and 
by  Davis’s  suggestion,  Clark  (1979).  Discussion  of  die. techniques  used  and 
die  operation  of  die  Geo-EAS  program  is  contained  in  Englund  and  Sparks 
(1988). 
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Figure  66.  Estimated  sand  thickness  from  the  Salt  River  1000  foot  grid 


7  Predictive  Model 


As  a  result  of  die  effort  described  in  the  preceding  sections,  a  model  was 
produced  which  predicts  die  location  of  data  points  needed  in  defining  a  dis¬ 
continuous  sand  body.  After  having  survived  trials,  as  described  in  the  Case 
Study,  Supplemental  Case  Study  and  Geostatistics  sections,  the  rationale  for 
exploration  has  become  die  predictive  model.  The  model  is  shown  in  Fig¬ 
ure  67  and  is  described  as  follows. 


The  Model 


The  model  as  shown  in  Figure  67  begins  initially,  as  any  site  investigation 
should,  with  a  literature  survey.  The  question  is  then  asked,  "Did  the  litera¬ 
ture  survey  provide  any  information  on  which  to  base  an  exploration  pro¬ 
gram?”.  If  not,  a  minimum  of  three  stratigraphic  borings  must  be  drilled  to 
obtain  die  minimum  information  needed,  the  direction  of  dip  of  the  bedrock  of 
interest. 

With  die  information,  either  from  the  literature  survey  or  from  boring 
information,  priorities  on  which  site  boundaries  should  be  drilled  first  are 
established.  The  priorities  are  based  on  die  most  likely  orientation  of  any 
possible  sand  body. 

Once  die  boundary  priorities  are  established,  a  decision  must  be  made  as  to 
whether  surface  geophysics  can  be  used  to  explore  for  any  sand  body  along 
the  boundaries.  If  surface  geophysics  is  a  possibility,  the  method  most  appli¬ 
cable  must  be  selected.  If  surface  geophysics  cannot  be  used,  a  drilling  pro¬ 
gram  must  be  initiated.  The  spacing  for  placement  of  die  borings  will  depend 
on  die  application  of  die  model. 

During  die  exploration  program  and/or  upon  completion  of  it,  die  question, 
"Was  a  sand  body  encountered?"  must  be  asked.  If  no  sand  was  encountered, 
the  question  then  posed  is  whether  die  exploration  program  is  complete  or  not. 
If  the  exploration  program  is  complete,  then  the  model  is  ended  with  no  sand 
body  having  been  found.  If  die  exploration  program  is  not  complete,  die 
model  is  reentered  by  continued  exploration  inward  from  the  boundary.  This 
cycle  or  loop  would  continue  until  no  sand  was  found  and  the  exploration 
program  was  completed,  or  until  a  sand  body  was  encountered. 
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Figure  67.  (Concluded) 


At  any  time,  if  a  sand  body  is  encountered,  die  question,  "Can  the  environ¬ 
ment  of  deposition  be  determined?*  must  be  asked.  This  requires  that  a  pos¬ 
sible  sand  body  found  by  surface  geophysical  methods  be  drilled,  to  obtain  the 
information  necessary  to  make  die  environment  of  deposition  determination. 

If  the  environment  of  deposition  can  be  determined,  it  must  be  determined  to 
be  that  of  a  meandering  stream.  If  it  is  not,  the  model  ends,  since  a  require¬ 
ment  for  this  model  is  that  the  sand  body  is  deposited  by  a  meandering 
stream.  If  die  environment  of  deposition  is  fluvial  and  meandering,  or  if  it  is 
assumed  to  be  fluvial  and  meandering  (the  environment  of  deposition  not 
determinable),  die  sand  body  width  is  estimated. 

The  next  step  in  die  model,  following  estimation  of  die  sand  body  width,  is 
dependent  upon  die  purpose  for  which  die  model  is  being  used.  If  the  model 
is  being  used  to  define  a  sand  body,  as  would  be  typical  for  an  environmental 
application,  a  data  or  well  spacing  that  has  a  SO  percent  probability  of  dual 
penetration  of  die  sand  body  is  chosen.  If  the  model  is  being  used  to  reen- 
counter  the  sand  body,  as  would  be  typical  for  an  oil  and  gas  application,  a 
data  or  well  spacing  that  has  the  highest  possible  probability  of  dual  penetra¬ 
tion  is  selected  within  the  spacing  restrictions  for  which  die  model  is  being 
used.  In  either  case,  die  data  or  well  spacing  is  picked  from  the  curves  in 
Figures. 

Using  die  data  or  well  spacing  distance  as  die  radius  from  well(s)  penetrat¬ 
ing  die  sand  body,  circle(s)  or  arc(s)  are  drawn. 

Using  die  Geo-EAS  program,  die  standard  deviations  for  estimated  sand 
thicknesaes  produced  by  kriging  is  contoured  from  the  current  data  set. 

Additional  data  point(s)  or  well(s)  are  placed  in  locations)  on  the  circle  or 
arc  where  die  standard  deviation  shows  that  die  thickness  estimates  standard 
deviation  is  greater  than  that  of  die  standard  deviation  for  the  whole  set  of 
data,  or,  where  die  thickness  estimates  standard  deviation  is  the  highest.  If 
more  than  one  data  location  is  placed  along  die  circle  or  arc,  they  are  sepa¬ 
rated  by  die  data  or  well  spacing  obtained  for  the  radius  of  the  circle  or  arc. 

The  next  step  in  the  model  is  answering  die  question,  "Has  die  data  set,  to 
date,  created  a  reasonable  estimate  of  thickness  for  the  area  of  interest  (i.e. 
die  sand  body)?”.  If  not,  die  model  must  be  reentered  with  the  question  "Was 
a  sand  body  encountered?”  (with  die  most  recent  addition  to  die  data  set).  It 
should  be  noted  that  when  reentered  at  this  point,  each  subsequent  portion  of 
die  model  is  redone.  This  includes  calculating  new  sand  body  widths  esti¬ 
mates,  selecting  new  data  or  well  spacings,  etc.  If  and/or  when  die  area  of 
interest  does  show  sand  thickness  standard  deviations  less  than  die  standard 
deviation  of  die  whole  data  set,  then  a  confidence  statement  is  made. 

The  confidence  statement  is  the  thickness  estimate  +  or  -  the  standard 
deviation  for  die  kriged  area.  If  die  standard  deviation  is  greater  than  the 
thickness,  then  it  cannot  be  said  that  the  sand  is  present.  In  this  case,  addi¬ 
tional  data  or  boring  locations  must  be  added,  in  locations  of  die  highest 
standard  deviations  and  then  reenter  die  model  by  making  a  new  confidence 
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statement.  Once  the  confidence  statement  is  reasonable,  (i.e.  the  standard 
deviation  is  less  than  the  estimated  thickness)  the  model  is  completed. 

If  a  sand  body  is  encountered,  meeting  the  conditions  required  to  take  the 
model  to  completion  will  normally  require  cycling  through  the  model  several 
times,  reentering  through  the  "Was  a  sand  body  encountered?"  section. 
Through  this  process,  die  exploration  program  will  follow  the  sand  body 
through  the  site. 


8  Conclusion 


The  research  described  in  this  report  was  directed  at  developing  a  method 
to  minimize  die  data  needed  to  define  a  discontinuous  sand  body.  Since  most 
populated  areas  are  located  adjacent  to  streams,  the  fluvially  deposited  sand 
body  was  die  specific  type  targeted  by  this  research.  The  purpose  of  the 
research  was  to  minimize  the  hazards  and  costs  associated  with  die  exploration 
of  such  sand  bodies,  particularly  at  contaminated  sites.  Minimizing  the  num¬ 
ber  of  borings  necessary  to  define  a  fluvial  sand  body  was  accomplished  by 
bringing  geology  into  the  exploration  by  predicting  the  location  for  needed 
data,  based  on  determining  die  geometry  of  a  sand  body,  once  it  was  encoun¬ 
tered,  and  by  using  any  data  already  available. 

The  geometry  of  die  sand  body  was  determined  by  establishing  the  envi¬ 
ronment  of  deposition  from  stratigraphic  data  such  as  lithology,  geophysical 
logs,  grain  size  analysis,  and  the  amount  of  quartz  present.  Once  die  environ¬ 
ment  of  deposition  was  determined  to  be  fluvial  and  meandering,  die  thickness 
of  the  sand  was  used  to  estimate  a  width  for  die  sand  body,  using  Leeder’s 
(1973)  and  Lorenz  and  others’  (1985)  equations.  This  procedure  required 
several  assumptions  and  contained  errors  (based  on  geological  variations),  but 
produced  a  satisfactory  estimate  of  die  sand  body  width.  From  the  sand  body 
width,  a  data  point  spacing  was  obtained  from  die  probability  of  penetrating 
die  sand  a  second  time. 

By  using  die  data  from  locations  already  in  existence,  such  as  that  created 
by  boundary  drilling  exploration,  die  thickness  was  kriged  using  die  Geo-EAS 
program.  This  produced  a  grid  of  estimates  of  die  sand  body  thicknesses,  and 
perhaps  more  importantly,  die  error  of  those  estimates  (as  die  standard  devia¬ 
tion).  By  comparing  die  standard  deviation  of  die  estimates  with  the  standard 
deviation  for  the  entire  data  set,  new  data  point  locations  could  be  chosen  in 
areas  where  die  sand  thickness  standard  deviations  were  greater  than  the  stan¬ 
dard  deviation  of  die  whole  data  set,  or  in  areas  with  the  most  error,  at  die 
data  point  spacing  needed. 

The  procedure  described  above  was  followed  until  die  sand  body  was 
defined.  By  using  this  method,  the  sand  body  was  defined  with  significantly 
fewer  data  point  locations  than  required  by  typical  grid  methods  of 
exploration. 
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This  rationale  for  exploration  was  applied  to  the  Rocky  Mountain  Arsenal, 
where  the  sand  body  was  defined  to  the  point  where  the  standard  deviation  of 
foe  sand  thicknesses  were  less  than  that  of  the  whole  data  set.  The  exploration 
already  conducted  at  die  Arsenal  had  not  defined  the  sand  body  to  a  point 
where  die  standard  deviation  of  die  sand  thicknesses  were  less  than  that  of  the 
whole  data  set,  thus,  a  planned  comparison  for  accuracy  could  not  be  made. 
This  resulted  in  a  supplemental  site  being  used  for  comparison  purposes. 

Since  die  Arsenal  site  supported  the  stratigraphic  determination  portion  of 
die  rationale  for  exploration,  the  supplemental  site  was  chosen  as  a  modern 
floodplain  whose  width  could  be  visually  established.  Using  this  known 
width,  a  thickness  was  obtained  and  used  as  the  thickness  of  the  sand  body  at 
all  locations.  The  site  chosen  was  that  of  the  Brazos  River.  Hypothetical 
exploration  programs  wore  then  conducted  to  define  the  sand  body,  using  the 
rationale  for  exploration  and  a  typical  grid  exploration  method.  The  results  of 
each  method  of  exploration  were  then  compared.  The  comparison  showed 
that  the  rationale  for  exploration  had  defined  the  hypothetical  sand  body  with 
significantly  fewer  data  points,  but  with  similar  accuracy  of  that  of  the  grid 
method. 

Variograms  used  to  krige  during  the  hypothetical  explorations  at  the 
Arsenal  site  and  for  die  Brazos  River  site  had  to  be  justified,  because,  nor¬ 
mally  large  data  sets  are  required  to  produce  die  variogram.  All  the  data 
defining  the  sand  body  at  the  Arsenal,  and  close  spaced  grid  data  from  the 
Brazos  River  and  die  Salt  River  were  used  to  create  variograms.  These  vario¬ 
grams  from  the  large  data  sets  indicated  that  variograms  used  in  the  limited 
data  explorations  were  justified. 

Having  survived  the  test  case  applications,  die  rationale  for  exploration 
became  the  predictive  model  for  defining  a  fluvial  sand  body. 

Thus,  based  on  the  research  described,  the  predictive  model  which  was 
developed  can  be  used  to  select  data  point  locations  in  defining  a  fluvial  sand 
body  by  following  die  sand  body  into  and  through  die  site.  This  is  accom¬ 
plished  by  bringing  geologically  based  statistical  methods  into  the  exploration 
program.  By  doing  this,  signiflcandy  fewer  data  points  are  needed  to  define 
the  sand  body  than  needed  by  the  typical  grid  method,  which  is  commonly 
used.  The  lowo-  number  of  data  points  reduces  die  hazard  and  cost  associated 
with  die  exploration  program. 
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Based  on  the  findings  of  this  research  die  following  recommendations  are 
made: 


a.  The  conclusion  that  die  predictive  model  defined  a  discontinuous  sand 
body  with  a  minimum  number  of  data  point  locations  is  based  on 
limited  testing.  Additional  sites  should  be  tested.  These  need  to  be  of 
various  size  fluvial  sand  bodies,  from  different  locations,  in  different 
climates,  of  different  ages,  for  varying  amounts  of  initial  data,  and  for 
different  applications. 

b.  Although  die  findings  of  this  research  jows  that  die  predictive  model 
reduces  die  number  of  data  point  locstioas  needed  to  define  a  discontin¬ 
uous  sand  body,  existing  or  visible  da*  was  used.  The  model  should 
be  used  during  actual  exploration  programs  for  determining  "in  the 
field"  applicability. 

c.  The  confidence  statements  which  were  made  for  die  defined  sand  bodies 
in  this  research  shows  large  errors,  particularly  at  die  smaller  thick¬ 
nesses.  These  confidence  statements  were  made  based  on  die  assump¬ 
tion  that  the  distribution  of  die  data  is  normal,  even  though  a  normal 
distribution  can  give  a  physically  impossible  sand  body  thickness  (nega¬ 
tive  about  a  zero  or  small  thickness).  This  is  because  die  functions 
which  may  better  describe  the  variable’s  distribution,  such  as  log¬ 
normal  or  beta,  need  a  relatively  large  amount  of  data  to  even  establish 
which  function  would  be  best.  More  data  than  this  application  usually 
will  usually  provide.  Additionally,  little  work  has  been  published  to 
advance  the  state  of  die  art  for  these  functions.  Further  investigation 
needs  to  address  the  use  of  non-normal  distribution  functions. 

d.  Investigations  should  be  conducted  to  determine  if  die  predictive  model 
developed  in  this  research  can  be  adapted  to  stacked  fluvial  meander 
belt  sand  bodies,  in  a  layered  sequence. 

e.  Other  environments  of  deposition  for  discontinuous  sand  bodies  need  to 
be  researched  in  order  to  determine  if  predictive  models  can  be  devel¬ 
oped  to  reduce  the  number  of  data  point  locations  needed  to  define 
hem. 
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Appendix  A 
Kriging  Tutorial 


Overview 

Geostatistical  methods  are  useful  for  site  assessment  where  data  are 
collected  on  a  spacial  network  of  sampling  locations.  Kriging  is  a  weighted 
moving  average  method  used  to  interpolate  values  from  a  sample  data  set  onto 
a  grid  of  points  for  contouring.  The  kriging  weights  are  computed  from  a 
variogram,  which  measures  die  degree  of  correlation  among  sample  values  in 
die  area  as  a  function  of  die  distance  and  direction  between  samples  (Englund 
and  Sparks  1988).  This  tutorial  will  provide  a  glossary  of  geostatistical  terms 
adapted  from  Englund  and  Sparks,  1988;  a  discussion  of  kriging,  with  die 
kriging  equation;  a  discussion  of  variograms,  with  explanations  of  die 
different  models;  and  a  simple  example  showing  calculations. 


Glossary 

Block  Kriging  -  Estimating  die  value  of  a  block,  centered  on  a  specific  grid 
node,  from  a  set  of  nearby  sample  values  using  kriging. 

Covariance  -  A  statistical  measure  of  die  correlation  between  two  variables. 
Covariance  is  usually  treated  as  die  simple  inverse  of  die  variogram, 
computed  as  die  overall  sample  variance  minus  die  variogram  value. 

Exponential  Model  -  A  function  frequently  used  when  fitting  mathematical 
models  to  experimental  variograms. 

Gaussian  Model  -  A  function  frequently  used  when  fitting  mathematical 
models  to  experimental  variograms. 

Kriging  Standard  Deviation  -  The  standard  error  of  estimation  computed 
for  a  kriged  estimate.  Kriging  is  the  weighted  linear  estimate  with  die  particu¬ 
lar  set  of  weights  which  minimizes  the  computed  estimation  variance  (standard 
error  squared).  The  relationship  of  die  kriging  standard  deviation  to  the  actual 
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error  of  estimation  is  very  dependent  on  the  variogram  model  used  and  die 
validity  of  the  underlying  assumptions. 


Linear  Model  -  A  function  frequently  used  when  fitting  mathematical 
models  to  experimental  variograms. 

Madogram  -  A  plot  of  mean  absolute  differ  Mice  of  paired  sample  measure¬ 
ments  as  a  function  of  distance  and  direction.  Madograms  are  not  true  vario¬ 
grams,  and  generally  should  not  be  used  in  kriging.  The  kriging  standard 
deviations  will  be  meaningless. 

Nugget  Model  -  A  constant  variance  model  most  often  used  in  combination 
with  one  or  more  other  functions  when  fitting  mathematical  models  to 
experimental  variograms. 

Ordinary  Kriging  -  A  variety  of  kriging  which  assumes  that  local  means 
are  not  necessarily  closed  related  to  the  population  mean,  and  which  therefore 
uses  only  die  samples  in  die  local  neighborhood  for  the  estimate.  Ordinary 
kriging  is  die  most  commonly  used  method  for  environmental  situations. 

Point  Kriging  -  Estimating  die  value  of  a  point  from  a  set  of  nearby  sample 
values  using  kriging.  The  kriged  estimate  for  a  point  will  usually  be  quite 
similar  to  die  kriged  estimate  for  a  relatively  small  block  centered  on  die 
point,  but  the  computed  kriging  standard  deviation  will  be  higher.  When  a 
kriged  point  happens  to  coincide  with  a  sampled  location,  die  kriged  estimate 
will  equal  die  sample  value. 

Range  -  The  distance  at  which  a  variogram  model  readies  its  maximum 
value,  or  sill. 

Semi-Variogram  -  There  is  disagreement  in  die  geostatistical  literature  as  to 
whether  "semi-variogram"  or  "variogram"  should  be  used,  but  they  have  die 
same  meaning. 

Sill  -  The  upper  limit  of  any  variogram  model. 

Simple  Kriging  -  A  variety  of  kriging  which  assumes  that  local  means  are 
relatively  constant  and  equal  to  die  population  mean,  which  is  known.  The 
papulation  mean  is  used  as  a  factor  in  each  local  estimate,  along  with  die 
samples  in  die  local  neighborhood. 

Spherical  Model  -  A  function  frequently  used  when  fitting  mathematical 
models  to  experimental  variograms. 

Variogram  -  A  plot  of  die  variance  (one-half  die  mean  squared  difference) 
of  paired  sample  measurements  as  a  function  of  die  distance  (and  optionally  of 
the  direction)  between  samples.  Typically,  all  possible  sample  pairs  are 
examined  and  grouped  into  classes  (called  lags)  of  approximately  equal  dis¬ 
tance  and  direction.  Variograms  provide  a  means  of  quantifying  the 
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commonly  observed  relationship  that  samples  close  together  will  tend  to  have 
more  values  than  samples  far  apart. 


Variograms 

Englund  and  Sparks  (1988)  state  that  the  computation,  interpretation,  and 
modeling  of  variograms  is  die  "heart"  of  a  geostatistical  study.  The  vario- 
gram  model  is  the  interpretation  of  die  spatial  correlation  structure  of  the  sam¬ 
ple  data  set.  It  controls  the  way  that  kriging  weights  are  assigned  to  samples 
during  interpolation,  and  consequendy  controls  the  quality  of  the  results. 

All  interpolation  and  contouring  methods  make  die  assumption  that  some 
type  of  spatial  correlation  is  present.  They  assume  that  a  measurement  at  any 
point  represents  nearby  locations  better  than  locations  farther  away.  Vario- 
gram  analysis  attempts  to  quantify  the  relationship  of  how  well  a  measurement 
can  be  expected  to  represent  another  location  a  specific  distance  away.  Vario¬ 
grams  plot  the  average  difference  (actually,  one-half  die  squared  difference,  or 
variance)  of  pairs  of  measurements  against  die  distances  separating  the  pairs. 

If  measurements  were  possible  at  all  sample  locations,  a  "true"  variogram 
could  be  computed  for  a  site  showing  die  variance  of  all  pairs  of  measure¬ 
ments  which  satisfy  each  combination  of  distance  and  direction.  In  practice, 
with  limited  data,  variances  are  computed  for  groups  of  pairs  of  measurements 
in  class  intervals  of  similar  distance  and  direction.  Then  a  graph  is  plotted  of 
die  variances  versus  distance.  Then  a  model  curve  is  fitted  to  die  graph.  The 
model  is  assumed  to  be  an  approximation  of  the  "true”  variogram. 

Davis  (1986)  describes  the  semivariogram  in  a  similar  manor,  exemplifying 
die  use  of  a  different  term  for  die  variogram.  Davis  (1986)  continues  his 
presentation  of  die  "semivariogram"  by  assuming  that  the  samples  are  point 
measurements  of  a  property.  For  computational  tractability,  die  assumption  is 
also  made  that  the  samples  are  uniformly  spaced  along  straight  lines.  If  the 
spacing  between  samples  along  a  line  is  some  distance  A,  the  semivariance, 

7h,  can  be  estimated  for  distances  that  are  multiples  of  A: 

Jk  -  £,"■*  <*,  - 


In  this  notation,  X;,  is  a  measurement  of  a  regionalized  variable  taken  at 
location  i,  and  Xj+j,  is  another  measurement  taken  h  intervals  away.  We  are 
therefore  finding  die  sum  of  the  squared  differences  between  pairs  of  points 
separated  by  die  distance  Ah.  The  number  of  points  is  n,  so  the  number  of 
comparisons  between  pairs  of  points  is  n  -  h. 

If  the  semivariances  are  calculated  for  different  values  of  h,  die  results  can 
be  plotted  in  the  form  of  a  semivariogram  (i.e.  variogram).  When  die  dis¬ 
tance  between  sample  points  is  ^ero,  the  value  at  each  point  is  being  compared 
with  itself.  Hence,  all  die  differences  are  zero,  and  die  semivariance  for  y0  is 
zero.  If  Ah  is  a  small  distance,  the  points  being  compared  trad  to  be  very 
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similar,  and  the  semivariance  will  be  a  small  value.  As  die  distance  Ah  is 
increased,  the  points  being  compared  are  less  and  less  closely  related  to  each 
other  and  their  differences  become  larger  resulting  in  larger  values  of  7b.  At 
some  distance  the  points  being  compared  are  so  far  apart  that  they  are  not 
related  to  each  other,  and  their  squared  differences  become  equal  in  magnitude 
to  the  variance  around  the  average  value.  The  semivariance  no  longer 
increases  and  die  semivariogram  develops  a  flat  region  called  a  sill.  The 
distance  at  which  the  semivariance  approaches  the  variance  is  referred  to  as 
the  range  (or  span)  of  die  regionalized  variable,  and  defines  a  neighborhood 
within  which  all  locations  are  related  to  one  another. 

For  some  arbitrary  point  in  space,  the  neighborhood  can  be  envisioned  as  a 
symmetrical  interval  about  the  point.  If  the  regionalized  variable  is  stationary, 
or  has  die  same  average  value  everywhere,  any  locations  outside  the  interval 
are  completely  independent  of  the  central  point,  and  cannot  provide  informa¬ 
tion  about  die  value  of  die  regionalized  variable  at  that  location.  Within  die 
neighborhood,  however,  the  regionalized  variable  at  all  observation  points  is 
related  to  die  regionalized  variable  at  die  central  location  and  hence  can  be 
used  to  estimate  its  value.  If  a  number  of  measurements  are  used,  made  at 
locations  within  die  neighborhood  to  estimate  die  value  of  the  regionalized 
variable  at  the  central  location,  die  semivariogram  provides  the  proper  weight¬ 
ings  to  be  assigned  to  each  of  these  measurements. 

The  semivariogram  expresses  the  spatial  behavior  of  die  regionalized  vari¬ 
able  or  its  residual.  A  reasonable  form  for  die  semivariogram  must  be 
assumed  and  used  as  a  first  approximation.  A  semivariogram  tangent  to  the  X 
axis  at  die  origin  is  described  as  parabolic  and  indicates  that  die  regionalized 
variable  is  exceptionally  continuous.  A  variogram  that  is  linear  in  form  indi¬ 
cates  moderate  continuity  of  the  regionalized  variable.  A  truly  random  vari¬ 
able  will  have  no  continuity  and  its  semivariogram  will  ba  a  horizontal  line 
equal  to  die  variance.  In  some  circumstances  die  semivariogram  will  appear 
to  not  go  through  die  origin  but  rather  will  assume  some  nonzero  value.  This 
is  referred  to  as  die  "nugget  effect".  In  theory,  yQ  must  equal  zero.  The 
nugget  effect  arises  because  die  regionalized  variable  is  so  erratic  over  a  very 
short  distance  that  die  semivariogram  goes  from  zer o  to  die  level  of  the  nug¬ 
get  effect  in  a  distance  less  than  the  sampling  interval. 

In  principle,  the  experimental  semivariogram  could  be  used  direedy  to 
provide  values  for  estimatioL  procedures.  However,  the  semivariogram  is 
known  only  at  discrete  points  representing  distances  Ah.  In  practice,  semi- 
variances  may  be  required  for  any  distance,  whether  a  multiple  of  A  or  not. 
For  this  reason,  die  discrete  experimental  semivariogram  must  be  modelled  by 
a  continuous  function  that  can  be  evaluated  for  any  desired  distance.  Fitting 
a  model  equation  to  an  experimental  semivariogram  is  a  trial-and -error  pro¬ 
cess,  usually  done  by  eye.  Ideally,  die  model  chosen  to  represent  the  semi¬ 
variogram  should  begin  at  die  origin,  rise  smoothly  to  some  upper  limit,  then 
continue  at  a  constant  level. 
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The  spherical  model  has  these  properties.  For  a  distance,  h,  less  than  die 
range,  R,  it  is  defined  as: 


Jtl 
2 R3 


for  all  distances  to  die  range,  R,  of  the  semivariogram.  Beyond  the  range, 
die  semivariance,  yh,  equals  the  variance,  yt.  The  spherical  model,  shown  in 
die  following  diagram,  usually  is  described  as  the  ideal  form  of  the  semivario¬ 
gram.  Another  that  is  sometimes  used  is  the  exponential  model: 

Jh  *  y*  (i -«"*'*) 


The  exponential  model,  shown  in  the  following  diagram,  never  quite  reaches 
die  limiting  value  of  die  sill,  but  approaches  it  asymptotically.  Also,  the 
semivariance  of  die  exponential  model  is  lower  than  the  spherical  for  all  val¬ 
ues  of  h  less  than  the  range. 

The  linear  model,  is  simpler  than  either  the  spherical  or  exponential, 
because  it  has  only  one  parameter,  die  slope  a.  The  model  has  the  form: 

7h  m 


and  plots  as  a  straight  line  through  the  origin.  This  model  cannot  have  a  sill, 
as  it  rises  without  limit.  Sometimes,  as  shown  in  die  following  diagram,  the 
linear  model  is  arbitrarily  modified  by  inserting  a  sharp  break  at  the  sill  value. 
The  use  of  such  a  model  has  been  criticized  because  die  kriging  estimation 
procedure  presumes  die  semivariogram  is  a  continuous  smoothly  varying 
function.  If  die  regionalized  variable  has  been  sampled  at  a  sufficient  density, 
relative  to  die  range,  there  will  be  no  significant  differences  between  estimates 
made  assuming  a  linear  model  and  those  obtained  using  a  spherical  or  other 
model. 


Kriging 

According  to  Davis  (1986)  kriging  addresses  a  regionalized  variable,  which 
is  a  naturally  occurring  property  that  has  characteristics  intermediate  between 
a  truly  random  variable  and  one  that  is  completely  deterministic.  Many  geo¬ 
logical  surfaces,  both  real  and  conceptual,  can  be  regarded  as  regionalized 
variables.  These  are  continuous  from  place  to  place  and  therefore  must  be 
spatially  correlated  over  short  distances.  However,  points  on  an  irregular 
surface  that  are  widely  separated  tend  to  be  statistically  independent.  The 
degree  of  spatial  continuity  of  a  regionalized  variable  1  be  expressed  by  a 
semivariogram.  If  measurements  have  been  made  at  scattered  sampling  points 
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and  die  form  of  die  semivariogram  is  known,  it  is  possible  to  estimate  the 
value  of  die  surface  at  any  unsampled  location.  This  estimation  procedure  is 
called  kriging,  named  after  K.G.  Krige,  a  South  African  mining  engineer  and 
pioneer  in  die  application  of  statistical  techniques  to  mine  evaluation. 

Kriging  can  be  used  to  make  contour  maps,  but  unlike  conventional  con¬ 
touring  algorithms,  it  has  certain  statistically  optimal  properties.  Perhaps 
most  importantly,  the  method  provides  measures  of  the  error  or  uncertainty  of 
die  contoured  surface.  Kriging  uses  the  information  from  the  semivariogram 
to  find  an  optimal  set  of  weights  that  are  used  in  the  estimation  of  die  surface 
at  unsampled  locations.  Since  the  semivariogram  is  a  function  of  distance,  the 
weights  change  according  to  the  geographic  arrangement  of  the  samples. 

Point,  or  Punctual  according  to  Davis  (1986),  is  the  simplest  form  of 
kriging,  in  which  the  observations  consist  of  measurements  taken  at  dimen¬ 
sionless  points,  and  the  estimates  are  made  at  other  locations  that  are 
dimensionless  points.  Punctual  kriging  is  used  in  contour  mapping  where  the 
observations  may  be  from  a  set  of  exploratory  drill  holes.  Constructing  a  map 
requires  that  estimations  of  the  variable  be  made  at  closely  spaced  locations 
over  the  map  area.  Once  made,  contour  lines  can  be  drawn  through  these 
estimates. 

To  simplify  the  operation,  it  can  be  assumed  that  the  variable  being 
mapped  is  statistically  stationary,  or  free  from  drift.  The  value  at  an 
unsampled  location  may  be  estimated  as  a  weighted  average  of  the  known 
observations,  of  weight  W.  The  value,  Y,  at  a  point,  p,  is  based  on  a  small 
set  of  nearby  known  points: 

Y  'p  «  IWfr 


It  is  expected  that  the  estimated  value,  Y’  ,  will  differ  somewhat  from  the  true 
(but  unknown)  value,  Yf  by  an  amount  tnat  is  called  the  estimation  error, 

€p: 


If  die  weights  used  in  the  estimation  equation  sum  to  one,  the  resulting  esti¬ 
mates  are  unbiased.  This  means  that,  over  a  great  many  estimations,  the 
average  error  will  be  zero,  as  overestimates  and  underestimates  will  tend  to 
cancel  one  another.  However,  the  estimates  may  scatter  widely  about  the 
correct  values.  This  scatter  can  be  expressed  as  the  error  variance,  S| , 

E  0”,-V 

- : - 


or  as  its  square  root,  the  standard  error  of  the  estimate,  Se : 
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As  noted,  it  seems  intuitively  reasonable  that  nearby  known  points  should 
be  most  influential  in  estimating  the  value  at  an  unsampled  location,  and  that 
more  distant  control  points  should  be  less  influential.  It  also  seems  reasonable 
to  expect  that  the  weights  used  in  the  estimation  process,  and  the  error  in  the 
estimate,  should  be  related  in  some  way  to  the  semivariogram. 

For  example,  to  estimated  the  value  of  Y  at  a  point  p  from  three  nearby 
points,  using  as  our  estimator  a  weighted  average  of  the  three  known  values: 

Y'p-  WlYl  *  W2Y2  *  w3y3 

The  weights  are  constrained  to  sum  to  one,  so  die  estimate  is  unbiased  if  there 
is  no  trend.  Suppose  that  weight  W,  is  chosen  to  ge  equal  to  1.0.  Then, 
weights  W2  and  W3  must  be  zero  and  die  estimate  at  p  is: 

Y  'p  -  1.0 Yj  ♦  0.072  ♦  O.OK3 


or 


The  estimation  error  is  simply  €  =  Yp  -  Yj,  since  Yj  is  the  estimate  Y’p.  If 
many  other  locations  like  Yp  are  estimated  from  points  arranged  in  a  manner 
spatially  similar  to  Y]t  the  estimation  variance  can  be  calculated  as  die  aver¬ 
age  squared  difference  between  these  pairs  of  points.  For  convenience,  these 
other  estimated  locations  may  be  called  Ypi  and  the  other  estimating  points 
Yu.  Then,: 

S,2  -  ill  Em  <1*  -  jy2 

The  estimation  variance  is  equal  to  twice  the  semivariance  for  a  distance  equal 
to  the  separation  between  points  Ypi  and  YH. 

A  particular  combination  of  weights  have  been  chosen  to  arrive  at  an  esti¬ 
mate  Y’p  and  to  determine  die  estimation  error.  There  are  an  infinite  number 
of  other  possible  combinations  of  weights  that  could  be  chosen,  each  of  which 
will  give  a  different  estimate  and  a  different  estimation  error.  There  is,  how¬ 
ever,  only  one  combination  that  will  give  a  minimum  estimation  error.  It  is 
this  unique  combination  of  weights  that  kriging  attempts  to  find. 

Deriving  die  kriging  equations  requires  calculus  and  will  not  be  considered 
here.  A  simple  discussion  is  contained  in  Clark  (1979)  and  a  complete 
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derivation  is  provided  by  Olea  (1975).  Optimum  values  for  the  weights  can 
be  found  by  solving  a  set  of  simultaneous  equations,  which  includes  values 
from  a  semivariogram  of  foe  variable  being  estimated.  The  weights  are  opti¬ 
mal  in  foe  sense  that  foe  resulting  estimates  are  unbiased  and  have  minimum 
estimation  variance.  No  other  linear  combination  of  foe  observations  can 
yield  estimates  that  have  a  smaller  scatter  around  their  true  values. 

To  make  a  kriged  estimate  of  foe  value  Y*  at  a  point  p  from  three  known 
observations,  Ylf  Y2,  and  Y3,  three  weights,  W,,  W2,  and  W3  must  be  found 
for  foe  kriging  equation.  To  find  these  requires  foe  solution  to  a  system  of 
three  simultaneous  equations: 

W,  y(hu)  *  W2  y(/»12)  +  W3  y(hl3)  =  y(hlp) 

W,  y(hi2)  +  W2  yihn)  +  W3  y(h^)  =  7<V 

W,  y(hl3)  +  W2  yih^)  +  W3  y(h33)  =  y(h3p) 

In  fois  notation,  7(hy)  is  foe  semivariance  over  a  distance  h  corresponding  to 

foe  separation  between  points  i  and  j.  For  example,  7(h13)  is  foe  semivariance 
for  a  distance  equal  to  that  between  known  points  1  and  3;  7(hlp)  is  the  semi¬ 
variance  for  a  distance  equal  to  that  between  known  point  1  ana  foe  location 
p,  where  foe  estimate  is  to  be  made.  The  left-hand  matrix  is  symmetrical 
because  h~  =  h^.  It  has  zeroes  along  foe  main  diagonal  because  hy  represents 
foe  distance  from  a  point  to  itself,  which  is  zero.  Assuming  foe  semivario¬ 
gram  goes  through  foe  origin,  foe  semivariance  for  zero  distance  is  zero. 
Values  of  foe  semivariance  are  taken  from  the  semivariogram,  which  must  be 
known  (or  estimated)  prior  to  kriging. 

A  fourth  equation  is  needed  to  ensure  that  foe  solution  is  unbiased,  by 
constraining  foe  weights  to  sum  to  one: 

W,  +  W2  *  W3  -  1.0 

This  gives  a  set  of  four  equations  but  only  three  unknowns.  Since  there  are 
more  equations  than  unknowns,  an  extra  degree  of  freedom  can  be  used  to 
assure  that  foe  solution  will  have  foe  minimum  possible  estimation  error.  This 
is  done  by  adding  a  slack  variable,  called  a  Lagrange  multiplier,  X,  to  foe 
equation  set.  The  complete  set  of  simultaneous  equations  has  foe  following 
appearance: 

W,  y{hn)  +  W2  y(hl2)  +  W3  y(hl3)  *  X  *  y(hlp) 

W,  7(X12)  +  W2  yih^)  +  W3  y{hx)  ♦  X  =  y%J 

W,  y(hl3)  *  W2  yihx)  *  W3  y{h33)  +  X  =  y(h3p) 

Wl  *  W2  *  W3  +0=1 

Rearranging  in  matrix  form. 
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y(hn)  7012)  t(A,3)  r 

X 

70, P’ 

70 12)  7022)  7023)  * 

• 

W* 

7  (V 

70ij)  7023)  7033)  1 

W3 

70*) 

1110 
m  m 

X 

1 

In  general  terms  the  matrix  equation  must  be  solved: 

[A]  •  [WJ  -  [B] 

for  die  vector  of  unknown  coefficients,  [W].  The  terms  in  matrix  [A]  and 
vector  [B]  are  taken  directly  from  die  semivariogram  or  from  die  mathemati¬ 
cal  function  that  describes  its  form.  Once  die  unknown  weights  have  been 
determined,  the  variable  at  location  p  is  estimated  by: 

Y  'p  «  W,y,  ♦  W2Y2  *  WjY3 

The  estimation  variance  is: 

S?  “  W,  y(hlp)  *  w2  y%J  *  w3  yih^  ♦  X 

The  variance  estimate  is  essentially  die  weighted  sum  of  die  semivariances  for 
die  distances  to  the  points  used  in  die  estimation,  plus  a  contribution  from  the 
X  coefficient  that  is  equivalent  to  a  constant  term.  Kriging  has  two  powerful 
advantages  over  conventional  estimation  procedures  such  as  those  used  for 
contour  mapping.  Kriging  produces  estimates  that,  on  average,  have  the 
smallest  possible  error,  and  also  produces  an  explicit  statement  of  the  magni¬ 
tude  of  this  error. 

If  the  assumption  is  made  that  die  errors  of  estimation  are  normally  dis¬ 
tributed  about  foe  true  value,  foe  standard  error  can  be  used  as  a  confidence 
band  around  foe  estimates.  The  probability  that  foe  true  value  at  point  p  is 
within  one  standard  error  above  or  below  foe  value  estimated  is  68  percent, 
and  foe  probability  is  95  percent  that  foe  true  elevation  lies  within  two  stan¬ 
dard  errors. 

Although  a  normal  distribution  can  give  a  physically  impossible  number  in 
foe  case  of  variables  which  cannot  have  negative  values,  it  is  still  used.  A 
normal  distribution  is  used  because  of  a  lack  of  documentation  for  other  func¬ 
tions,  such  as  log-normal  and  beta,  even  though  they  may  give  a  more  realis¬ 
tic  distribution.  Also,  a  significant  amount  of  data  is  necessary  to  establish 
which  function  would  give  foe  most  realistic  distribution. 
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**»)  7(*i2)  y(hn)  i* 

7(V 

y(H  12)  yQhi )  7(^23)  1 

• 

*2 

S 

7<V 

7<M  7(^33)  1 

7(V 

1  1  i  0 

X 

1 

then: 


7(X,-*,)  7(*2’-*l)  7(*3‘xl) 

X 

7(*o-*i)' 

7(*r*2>  ‘V<x2“x2>  7(*3“jc2>  1 

• 

W2 

7(*o-JC2  ) 

7<*r*3)  7(*2“*3)  7U3-X3)  1 

*3 

7(*0-JC3) 

!  1  !  0 

X 

1 

so: 


7(0)  7(8)  7(6)  r 

7(10)* 

7(8)  7(0)  7(10)  1 

*2 

7(6) 

7(6)  7(10)  7(10)  1 

*3 

7(8) 

11  10 

X 

1 

where: 


7h  =  2h 
so: 


and: 


>(0)  -  0 
7(6)  -  12 
7(8)  -  16 
7(10)  -  20 


'0  16  12  r 

*20* 

16  0  20  1 

*2 

12 

12  20  0  1 

*3 

SB 

16 

1  1  10 

X 

1 

Invert  the  matrix  to  solve  for  Wlf  W2,  W3,  and 

where  the  W’«  are  the  kriging  weights  and  X  is  die  Lagrange  multiplier. 
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K 

'-.05681  .022727  .34090  .227272' 

'20' 

*2 

.022727  -.03409  .011363  .409090 

12 

*3 

.034090  .011363  .04545  .363636 

16 

X 

.227272  .409050  .363636  -10.9090 

■  _ 

1 

w,  -  -.09090 
W2  -  .636363 
W3  -  .454545 
X  -4.363636 


Note  that  the  W’s  sum  to  1  for  the  unbiased  requirement. 
Kq  —  WjATj  ♦  W2K2  +  W3ZT3 


so: 


Ko  -  123.64 


To  obtain  die  kriging  variance: 

£  -  Wl  7(*i P  ♦  W2  7<V  +  *3  7<V  ♦  X 


or: 


■  wi  7<*r*o>  *  w 2  yO^-Xq)  ♦  W3  7(x3-Io)  ♦  X 

which  is: 


-  W,  T(10)  ♦  W2  7(6)  +  W3  7(8)  «■  X 

-  -.09090(20)  ♦  .6363(12)  +  .4545(16)  ♦  4.36 

-  13.09 


and: 


S(  -  3.62 
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If  K  has  a  normal  distribution  then 

JT^  «  Kq±  2Sf  with  95  percent  probability 
Kq  *  2St  -  123.64  ♦  2(3.62)  -  130.88 
Kq  -  25,  -  123.64  -  2(3.62)  -  116.38 

so: 

116.38  £  a  130.88  with  95  percent  confidence 
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Apandix  B  Data  from  ProbatoWty  of  Dual  Ponotrathm  Caleulationa 


B1 


|  Probably 

f  of  Dual  Penetration 

Ml 

uSSSSM 

Sand  body  1  meander  bah)  width,  ft 

60 

100 

HI 

m 

760^ 

1,000 

1.250 

1,500 

1,760 

2j00^ 

2.600 

260 

.129 

.363 

.682 

.788 

.841 

.873 

.894 

.909 

.921 

.937 

500 

.032 

.064 

.163 

.363 

.676 

.682 

.745 

.788 

.818 

.841 

.873 

1,000 

.016 

.032 

.080 

.163 

.253 

.363 

.491 

.576 

.636 

.682 

.745 

1,500 

.011 

.021 

.053 

.107 

.163 

.221 

.285 

.363 

.454 

.523 

.618 

2,000 

.008 

.016 

.040 

.080 

.121 

.163 

.206 

.253 

.303 

.363 

.491 

2,600 

.006 

.013 

.032 

.064 

.096 

.129 

.163 

.197 

.234 

.272 

.363 

3,000 

.005 

.011 

.027 

.053 

.080 

.107 

.135 

.163 

.192 

.221 

.285 

3,500 

.005 

.009 

.023 

.046 

.068 

.092 

.115 

.139 

.163 

.187 

.239 

4,000 

.004 

.008 

.020 

.040 

.060 

mu 

.100 

.121 

.142 

.163 

.206 

4,500 

.004 

.007 

.018 

.035 

.053 

.089 

.107 

.125 

.144 

.182 

5,000 

.003 

.006 

.016 

.032 

.048 

.064 

.Q80 

.096 

.113 

.129 

.163 

6,000 

.003 

.005 

.013 

.027 

.040 

.053 

.067 

.080 

.094 

.107 

.135 

7,000 

B9 

.005 

.011 

.023 

.034 

.046 

.057 

.068 

.080 

.092 

.115 

8,000 

.002 

.004 

.010 

.020 

.030 

.040 

.050 

.060 

.070 

.080 

m 

9,000 

HI 

.004 

Hi 

.018 

.027 

.035 

.044 

.053 

.062 

.071 

.089 

10,000 

.002 

.008 

.016 

.024 

.040 

.048 

.056 

.064 

.080 

15,000 

.001 

El 

.006 

.011 

.016 

ESI 

.027 

.032 

.037 

.043 

.053 

20,000 

.002 

.004 

.008 

.012 

.020 

.024 

.028  | 

.032 

.040 

25,000 

m 

.003 

.010 

ESI 

.016 

.019 

.022 

.025 

.032 

30,000 

.001 

.003 

.005 

.008 

1BH 

.013 

.016 

.019 

.021 

.027 
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Appendix  C  Software  Verification* 


Cl 


SANO-BOOY  WIDTH  Wm  (It.) 


WEa  SPACING  (w) 


w.4000ti. 


Plot  of  sand  body  meander  belt  width  /w  »  „ 

tong penetrated  by  both  wolla,  fj ^^bab"ilv  of  the  sand  body 
,98SI  ”eC“<l  we"  spocings  M  (from  Ms' 


L  SAND-BODY  WIDTH  (ft,) 

^  sandstone  meender  bfilt  wirfth  aa/  t 

STt^r  w*"5- i  5 ^2"; 
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Contour  of  percent  iron  kriged  using  Geo-Eas,  with  Clark's  simulated  iron  ore 
deposit  data 


Contour  of  standard  deviation  of  the  kriging  estimates  of  percent  iron  using 
Geo-Eas,  with  Clark's  simulated  iron  ore  deposit  data 
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Sand  Body  Width  Estimation 

The  measured  thickness  of  19.0  ft  is  multiplied  by  1.1  to  compensate  for 
compaction. 

19.0  ft  X  1.1  -  20.9  ft 


The  20.9  ft  is  then  converted  to  meters  by  multiplying  by  3.048  X 
101  mm. 

20.9  ft  X  (3.048  X  10*1  m/ft)  *  6.37  m 

This  6.37  meters  is  dim  inserted  into  Leeder’s  formula  for  calculating 
channel  width  (W^. 

We  -  6.8ft1 54 

Where  h  is  die  thickness  of  6.37  meters. 

Wc  -  6.8(6.37 m)154 
We  «  117.73  m 


The  117.73  meters  is  converted  to  ft  by  multiplying  by  3.281  ft/m. 
117.73m  X  3.281  film  -  386.27  ft 

Appwtdbc  D  Sample  Calculation* 


This  386.27  ft  is  inserted  into  Lorenz  et  al.’s  equation  for  calculating 
meander  belt  width  (WJ 

WM  -  7.44WJ01 

Wm  -  7.44(386.27  ft)lM 

Wm  -  3050  ft 

The  calculated  meander  belt  width  of  3050  ft  is  the  estimate  of  the  sand 
body  width. 
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Append*  E  Data  Point  Location* 


Data  Points  for  tha  Rational  for  Exploration  at  the  Rocky  Mountain  [ 

Arsenal  | 

Grid  Coordinates 

Sand  body 

in  Feet 

Thickness 

Bast 

North 

in  Feet 

2178931 

175779 

0 

2183891 

175445 

0 

2187216 

175608 

0 

i  2188353 

185171 

0 

2188139 

180921 

0 

2178446 

179361 

0 

2178426 

185575 

0 

2184000 

185000 

19.0 

2181155 

184639 

0 

2183065 

182320 

10.5 

2186014 

184035 

0 

2183045 

180900 

28.6 

2184234 

177874 

10.4 

2181060 

178469 

0 

2186235 

180686 

0 

2181205 

180372 

0 

2181444 

176282 

0 

2186789 

178332 

0 

|  May's  Data  Points  at  tha  Rocky  Mountain  Arsenal  1 

Grid  Coordinates 

Sand  body 

in  Feet 

Thickness 

East 

in  Feet 

2184469 

177203 

46.0 

2183284 

180386 

44.0 

2178561 

184632 

0 

2183045 

180900 

28.6 

2182621 

161912 

21.1 

2183308 

183082 

0 

2183527 

182723 

37.8 

2183023 

183552 

0 

2162572 

183996 

0 

2184218 

185348 

10.0 

2185102 

185785 

0~ 

2183794 

185108 

19.0 

2184285 

181477 

0 
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|  AM  Data  Points  at  tha  Rocky  Mountain  Amend 

Grid  C 

oordinataa 

Sand  body 

Grid  Coordinates 

Sand  body 

in 

Fact 

Thlckaaaa 

in 

Faat 

Thlckaaaa 

last 

North 

la  Faat 

Eaat 

North 

in  Faat 

2164469 

177203 

46.0 

2181155 

184639 

0 

2183264 

160386 

44.0 

2186014 

184035 

0 

2183045 

180900 

28.6 

2181060 

178469 

0 

2182621 

161912 

21.1 

2166235 

180686 

0 

2163927 

162723 

37.8 

2181205 

180372 

0 

2184218 

185348 

1.0 

2186789 

178332 

0 

2183794 

185108 

19.0 

2178561 

184632 

0 

2183841 

179243 

38.9 

2183308 

183082 

0 

2183065 

182329 

10.5 

2183023 

183552 

0 

2184127 

185118 

6.0 

2182572 

183996 

0 

2184378 

184870 

5.5 

2185102 

185785 

0 

2184128 

184868 

16.3 

2184285 

181477 

0 

2183878 

164867 

4.8 

2183928 

184867 

5.8 

2183880 

184617 

9.2 

2183877 

185067 

15.0 

2184130 

184618 

9.1 

2184127 

185018 

14.0 

2184380 

184620 

3.1 

2183894 

182367 

27.4 

2183780 

184666 

6.2 

2183893 

182617 

31.3 

2178931 

175779 

0 

2183891 

182867 

33.0 

2183891 

175445 

0 

2183889 

183117 

25.2 

2187216 

175608 

0 

2183900 

181368 

40.9 

2188353 

185171 

0 

2184150 

184641 

6.0 

2188139 

180921 

0 

2185002 

181127 

3.1 

2178446 

179361 

0 

2184171 

160698 

29.0 

2176426 

185575 

0 

2184234 

177874 

10.4 

2184129 

184218 

10.0 

2181444 

176282 

0 
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Data  Points  lor  the  Rational  for  Exploration  on  the  Brazos  River 


Or id  Coord in* t os 

Sand  body 

Grid  Coordinates 

Sand  body 

la  Foot 

Thickness 

in 

Feet 

Thickness 

loot 

North 

In  Feet 

last 

North 

in  Feet 

3275000 

300000 

60 

3160000 

459000 

60 

3275000 

325000 

0 

3167000 

403000 

0 

3275000 

350000 

0 

3227500 

375000 

0 

3275000 

375000 

0 

3160000 

478000 

0 

3275000 

400000 

0 

3131000 

482500 

80 

3275000 

425000 

0 

3100000 

300000 

0 

3275000 

450000 

0 

3100000 

325000 

0 

3275000 

475000 

0 

3100000 

350000 

0 

3275000 

500000 

0 

3100000 

37SOOO 

0 

3250000 

305500 

80 

3100000 

400000 

0 

3242000 

314500 

60 

3100000 

425000 

0 

3225000 

325000 

80 

3100000 

450000 

0 

3225000 

342000 

80 

3100000 

475000 

0 

3213000 

356500 

80 

3100000 

500000 

0 

3200500 

374000 

0 

3125000 

300000 

0 

3196000 

361500 

80 

3125000 

500000 

80 

3177500 

361500 

0 

3150000 

300000 

0 

3190000 

377000 

80 

3150000 

500000 

0 

3103000 

393000 

80 

3175000 

300000 

0 

3103000 

411500 

80 

3175000 

500000 

0 

3107500 

429000 

0 

3200000 

275000 

0 

3170000 

425000 

80 

3200000 

300000 

0 

3160000 

440000 

80 

3200000 

500000 

0 

3225000 

275000 

0 

3225000 

500000 

0 

3250000 

275000 

0 

3250000 

500000 

0 

3275000 

275000 

0 
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Grid  Method  of  Exploration  on  the  Brazos  River 

f  Grid  Coordinates 

Saad  body 

Grid  Coordinates 

Sand  body 

la 

Fast 

Thickness 

In 

Feet 

Thickness 

East 

North 

la  Fast 

East 

North 

in  Feet 

3100000 

300000 

0 

3175000 

425000 

80 

3100000 

325000 

0 

3175000 

450000 

80 

3100000 

350000 

0 

3175000 

475000 

0 

3100000 

375000 

0 

3175000 

500000 

0 

3100000 

400000 

0 

3200000 

275000 

0 

3100000 

425000 

0 

3200000 

32SOOO 

0 

3100000 

450000 

0 

3200000 

350000 

80 

3100000 

475000 

0 

3200000 

375000 

80 

3100000 

500000 

0 

3200000 

400000 

0 

3125000 

300000 

0 

3200000 

425000 

0 

3125000 

450000 

0 

3200000 

450000 

0 

3125000 

475000 

0 

3200000 

500000 

0 

3125000 

500000 

60 

3225000 

275000 

0 

3150000 

300000 

0 

3225000 

300000 

0 

3150000 

400000 

0 

3225000 

325000 

60 

3150000 

425000 

0 

3225000 

350000 

80 

3150000 

450000 

80 

3225000 

375000 

0 

3150000 

475000 

80 

3225000 

500000 

0 

|  3150000 

500000 

0 

3250000 

275000 

0 

3175000 

300000 

0 

3250000 

300000 

80 

3175000 

350000 

0 

3250000 

325000 

80 

3175000 

375000 

0 

3250000 

350000 

0 

3175000 

400000 

80 

3250000 

500000 

0 

3275000 

275000 

0 

3275000 

300000 

80 

3275000 

325000 

0 

3275000 

350000 

0 

3275000 

375000 

0 

3275000 

400000 

0 

3275000 

425000 

0 

3275000 

450000 

0 

3275000 

475000 

0 

3275000 

500000 

0 
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Data  Points  for  the  Brazos  River 
10.000  ft  Grid 


Grid  Coordinates 
ia  Faat 

Sand  body 
Thickness 

Grid  Coordinates 
in  Faat 

Sand  body 
Thickness 

Bast 

North 

in  Faat 

East 

North 

in  Feet 

3100000 

300000 

0 

3110000 

330000 

0 

3100000 

310000 

0 

3110000 

340000 

0 

3100000 

320000 

0 

3110000 

350000 

0 

3100000 

330000 

0 

3110000 

360000 

0 

3100000 

340000 

0 

3110000 

370000 

0 

3100000 

350000 

0 

3110000 

380000 

0 

3100000 

360000 

0 

3110000 

390000 

0 

3100000 

370000 

0 

3110000 

400000 

0 

3100000 

380000 

0 

3110000 

410000 

0 

3100000 

390000 

0 

3110000 

420000 

0 

3100000 

400000 

0 

3110000 

430000 

0 

3100000 

410000 

0 

3110000 

440000 

0 

3100000 

420000 

0 

3110000 

450000 

0 

3100000 

430000 

0 

3110000 

460000 

0 

3100000 

440000 

0 

3110000 

470000 

0 

3100000 

450000 

0 

3110000 

480000 

0 

3100000 

460000 

0 

3110000 

490000 

0 

3100000 

470000 

0 

3110000 

500000 

0 

3100000 

480000 

0 

3120000 

300000 

0 

3100000 

490000 

0 

3120000 

310000 

0 

3100000 

500000 

0 

3120000 

320000 

0 

3110000 

300000 

0 

3120000 

330000 

0 

3110000 

310000 

0 

3120000 

340000 

0 

3110000 

320000 

0 

3120000 

350000 

0 

3120000 

360000 

0 

3130000 

390000 

0 

3120000 

370000 

0 

3130000 

400000 

0 

3120000 

380000 

0 

3130000 

410000 

0 

3120000 

390000 

0 

3130000 

420000 

0 

3120000 

400000 

0 

3130000 

430000 

0 

3120000 

410000 

0 

3130000 

440000 

0 

3120000 

420000 

0 

3130000 

450000 

0 

3120000 

430000 

0 

3130000 

460000 

0 

3120000 

440000 

0 

3130000 

470000 

0 

3120000 

450000 

0 

3130000 

480000 

80 

3120000 

460000 

0 

3130000 

490000 

80 

3120000 

470000 

0 

3130000 

500000 

80 

3120000 

480000 

0 

3140000 

300000 

0 

3120000 

490000 

80 

3140000 

310000 

0 

3120000 

500000 

80 

3140000 

320000 

0 

3130000 

300000 

0 

3140000 

330000 

0 

3130000 

310000 

0 

3140000 

340000 

0 

3130000 

320000 

0 

3140000 

350000 

0 

3130000 

330000 

0 

3140000 

360000 

0 

3130000 

340000 

0 

3140000 

370000 

0 

3130000 

350000 

0 

3140000 

380000 

0 

3130000 

360000 

0 

3140000 

390000 

0 
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Appendix  E  Data  Point  Locations 


Data  Points  for  the  Brazos  River 
10,000  ft  Grid 


Grid  Coord inataa 
in  Foot 

Sand  body 
Thickness 

Grid  Coordinates 
in  Feet 

Sand  body 
Thickness 

Boat 

North 

in  Feat 

East 

North 

in  Feet 

3130000 

370000 

0 

3140000 

400000 

0 

3130000 

380000 

0 

3140000 

410000 

0 

3140000 

420000 

0 

3150000 

450000 

80 

3140000 

430000 

0 

3150000 

460000 

80 

3140000 

440000 

0 

3150000 

470000 

80 

3140000 

450000 

0 

3150000 

480000 

0 

3140000 

460000 

80 

3150000 

490000 

0 

3140000 

470000 

80 

3150000 

500000 

0 

3140000 

480000 

80 

3160000 

300000 

0 

3140000 

490000 

80 

3160000 

310000 

0 

3140000 

500000 

80 

3160000 

320000 

0 

3150000 

300000 

0 

3160000 

330000 

0 

3150000 

310000 

0 

3160000 

340000 

0 

3150000 

320000 

0 

3160000 

350000 

0 

3150000 

330000 

0 

3160000 

360000 

0 

3150000 

340000 

0 

3160000 

370000 

0 

3150000 

350000 

0 

3160000 

380000 

0 

3150000 

360000 

0 

3160000 

390000 

0 

3150000 

370000 

0 

3160000 

400000 

0 

3150000 

380000 

0 

3160000 

410000 

0 

3150000 

390000 

0 

3160000 

420000 

0 

3150000 

400000 

0 

3160000 

430000 

0 

3150000 

410000 

0 

3160000 

440000 

80 

3150000 

420000 

0 

3160000 

450000 

80 

3150000 

430000 

0 

3160000 

460000 

80 

31500G0 

440000 

0 

3160000 

470000 

80 

3160000 

480000 

0 

3180000 

300000 

0 

3160000 

490000 

0 

3180000 

310000 

0 

3160000 

500000 

0 

3180000 

320000 

0 

3170000 

300000 

0 

3180000 

330000 

0 

3170000 

310000 

0 

3180000 

340000 

0 

3170000 

320000 

0 

3180000 

350000 

0 

3170000 

330000 

0 

3180000 

360000 

0 

3170000 

340000 

0 

3180000 

370000 

0 

3170000 

350000 

0 

3180000 

380000 

80 

3170000 

360000 

0 

3180000 

390000 

80 

3170000 

370000 

0 

3180000 

400000 

80 

3170000 

380000 

0 

3180000 

410000 

80 

3170000 

390000 

0 

3180000 

420000 

80 

3170000 

400000* 

0 

3180000 

430000 

80 

3170000 

410000 

80 

3180000 

440000 

80 

3170000 

420000 

80 

3180000 

450000 

0 

3170000 

430000 

80 

3180000 

460000 

0 

3170000 

440000 

80 

3180000 

470000 

0 

3170000 

450000 

80 

3180000 

480000 

0 

3170000 

460000 

0 

3180000 

490000 

0 

3170000 

470000 

0 

3180000 

500000 

0 
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Data  Points  for  the  Brazos  River 
10.000  ft  Grid 


1  Grid  Coord inataa 

1  in  Foot 

Sand  body 
Thickneaa 

Grid  Coordinates 
in  Feet 

Sand  body 
Thickneaa 

last 

North 

in  Faat 

Seat 

North 

in  Feet 

3170000 

460000 

0 

3190000 

300000 

0 

3170000 

490000 

0 

3190000 

310000 

0 

3170000 

500000 

0 

3190000 

320000 

0 

3190000 

330000 

0 

3200000 

340000 

80 

3190000 

340000 

0 

3200000 

350000 

80 

3190000 

350000 

0 

3200000 

360000 

80 

3190000 

360000 

80 

3200000 

370000 

80 

3190000 

370000 

80 

3200000 

380000 

0 

3190000 

380000 

SO 

3200000 

390000 

0 

3190000 

390000 

80 

3200000 

400000 

0 

3190000 

400000 

80 

3200000 

410000 

0 

3190000 

410000 

80 

3200000 

420000 

0 

3190000 

420000 

80 

3200000 

430000 

0 

3190000 

430000 

0 

3200000 

440000 

0 

3190000 

440000 

0 

3200000 

450000 

0 

3190000 

450000 

0 

3200000 

460000 

0 

3190000 

460000 

0 

3200000 

470000 

0 

3190000 

470000 

0 

3200000 

480000 

0 

3190000 

480000 

0 

3200000 

490000 

0 

3190000 

490000 

0 

3200000 

500000 

0 

3190000 

500000 

0 

3210000 

280000 

0 

3200000 

280000 

0 

3210000 

290000 

0 

3200000 

290000 

0 

3210000 

300000 

0 

3200000 

300000 

0 

3210000 

310000 

0 

3200000 

310000 

0 

3210000 

320000 

0 

3200000 

320000 

0 

3210000 

330000 

80 

3200000 

330000 

0 

3210000 

340000 

80 

3210000 

350000 

80 

3220000 

360000 

80 

3210000 

360000 

80 

3220000 

370000 

0 

3210000 

370000 

0 

3220000 

380000 

0 

3210000 

380000 

0 

3220000 

390000 

0 

3210000 

390000 

0 

3220000 

400000 

0 

3210000 

400000 

0 

3220000 

410000 

0 

3210000 

410000 

0 

3220000 

420000 

0 

3210000 

420000 

0 

3220000 

430000 

0 

3210000 

430000 

0 

3220000 

440000 

0 

3210000 

440000 

0 

3220000 

450000 

0 

3210000 

450000 

0 

3220000 

460000 

0 

3210000 

460000 

0 

3220000 

470000 

0 

3210000 

470000 

0 

3220000 

480000 

0 

3210000 

480000 

0 

3220000 

490000 

0 

3210000 

490000 

0 

3220000 

500000 

0 

3210000 

500000 

0 

3230000 

280000 

0 

3220000 

280000 

0 

3230000 

290000 

0 

3220000 

290000 

0 

3230000 

300000 

0 

3220000 

300000 

0 

3230000 

310000 

0 

3220000 

310000 

0 

3230000 

320000 

80 
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Appendix  E  Data  Point  Locations 


Data  Points  for  the  Brazos  Rivar 
10.000  ft  Grid 


Grid  Coordinates 
in  Faat 

East  North 

Sand  body 
Thickness 
in  Feet 

Grid  Coordinates 
in  Faet 

East  North 

Sand  body 
Thickness 
in  Feet 

3220000 

32000Q 

80 

3230000 

330000 

80 

3220000 

330000 

80 

3230000 

340000 

80 

3220000 

340000 

80 

3230000 

350000 

0 

3220000 

350000 

80 

3230000 

360000 

0 

3230000 

370000 

0 

3240000 

380000 

0 

3230000 

380000 

0 

3240000 

390000 

0 

3230000 

390000 

0 

3240000 

400000 

0 

3230000 

400000 

0 

3240000 

410000 

0 

3230000 

410000 

0 

3240000 

420000 

0 

3230000 

420000 

0 

3240000 

430000 

0 

3230000 

430000 

0 

3240900 

440000 

0 

3230000 

440000 

0 

3240000 

450000 

0 

3230000 

450000 

0 

3240000 

460000 

0 

I  3230000 

460000 

0 

3240000 

470000 

0 

1  3230000 

470000 

0 

3240000 

480000 

0 

1  3230000 

480000 

0 

3240000 

490000 

0 

3230000 

490000 

0 

3240000 

500000 

0 

3230000 

500000 

0 

3250000 

280000 

0 

3240000 

280000 

0 

3250000 

290000 

0 

3240000 

290000 

0 

3250000 

300000 

80 

3240000 

300000 

80 

3250000 

310000 

80 

3240000 

310000 

80 

3250000 

320000 

80 

3240000 

320000 

80 

3250000 

330000 

80 

3240000 

330000 

80 

3250000 

340000 

0 

3240000 

340000 

80 

3250000 

3.50000 

0 

3240000 

350000 

0 

3250000 

360000 

0 

3240000 

360000 

0 

3250000 

370000 

0 

3240000 

370000 

0 

3250000 

360000 

0 

|  3250000 

390000 

0 

3260000 

400000 

0 

3250000 

400000 

0 

3260000 

410000 

0 

3250000 

410000 

0 

3260000 

420000 

0 

3250000 

420000 

0 

3260000 

430000 

0 

3250000 

430000 

0 

3260000 

440000 

0 

3250000 

440000 

0 

3260000 

450000 

0 

3250000 

450000 

0 

3260000 

460000 

0 

3250000 

460000 

0 

3260000 

470000 

0 

3250000 

470000 

0 

3260000 

480000 

0 

3250000 

480000 

0 

3260000 

490000 

0 

3250000 

490000 

0 

3260000 

500000 

0 

3250000 

500000 

0 

3270000 

280000 

0 

3260000 

280000 

0 

3270000 

290000 

80 

3260000 

290000 

0 

3270000 

300000 

80 

3260000 

300000 

80 

3270000 

310000 

80 

3260000 

310000 

80 

3270000 

320000 

0 

3260000 

320000 

0 

3270000 

330000 

0 

3260000 

330000 

0 

3270000 

340000 

0 

3260000 

340000 

0 

3270000 

350000 

0 
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Data  Points  for  the  Brazos  River 

10.000  ft  Grid 

Grid  Coord inat as 

In  Fast 

Sand  body 
Thickness 
in  Feat 

Grid  Coordinate a 
in  Feet 

Sand  body 
Thickness 
in  Feet 

Bast 

North 

East 

North 

3260000 

350000 

0 

3270000 

360000 

0 

3260000 

360000 

0 

3270000 

370000 

0 

3260000 

370000 

0 

3270000 

380000 

0 

3260000 

380000 

0 

3270000 

390000 

0 

3260000 

390000 

0 

3270000 

400000 

0 

3270000 

410000 

0 

3280000 

420000 

0 

3270000 

420000 

0 

3280000 

430000 

0 

3270000 

430000 

0 

3280000 

440000 

0 

3270000 

440000 

0 

3280000 

450000 

0 

3270000 

450000 

0 

3280000 

460000 

0 

3270000 

460000 

0 

3280000 

470000 

0 

3270000 

470000 

0 

3280000 

480000 

0 

3270000 

480000 

0 

3280000 

490000 

0 

3270000 

490000 

0 

3280000 

500000 

0 

3270000 

500000 

0 

3280000 

280000 

0 

3280000 

290000 

80 

3280000 

300000 

80 

3280000 

310000 

0 

3280000 

320000 

0 

3280000 

330000 

0 

3280000 

340000 

0 

3280000 

350000 

0 

3280000 

360000 

0 

3280000 

370000 

0 

3280000 

380000 

0 

3280000 

390000 

0 

3280000 

400000 

0 

3280000 

410000 

0 
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Data  Points  for  the  Salt  River 
1,000  ft  Grid 


Grid  Coordinates 
in  Feet 

Kaet  North 

Sand  body 

Thickness 
in  Feet 

Grid  Coordinates 
in  Feat 

Bast  North 

Sand  body 
Thickness 
in  Feet 

435000 

875000 

17.5 

437000 

877000 

17.5 

435000 

876000 

17.5 

437000 

878000 

17.5 

435000 

877000 

17.5 

437000 

879000 

0 

435000 

878000 

17.5 

437000 

880000 

0 

435000 

879000 

0 

437000 

881000 

0 

435000 

880000 

0 

437000 

882000 

0 

435000 

881000 

0 

437000 

883000 

0 

435000 

882000 

0 

437000 

884000 

0 

435000 

883000 

0 

437000 

885000 

0 

435000 

884000 

0 

438000 

875000 

0 

435000 

885000 

0 

438000 

876000 

17.5 

436000 

875000 

0 

438000 

877000 

17.5 

436000 

876000 

17.5 

438000 

878000 

17.5 

436000 

877000 

17.5 

438000 

879000 

0 

436000 

878000 

17.5 

438000 

880000 

0 

436000 

879000 

0 

438000 

881000 

0 

436000 

880000 

0 

438000 

882000 

0 

436000 

881000 

0 

438000 

883000 

0 

436000 

882000 

0 

438000 

884000 

0 

436000 

883000 

0 

438000 

885000 

0 

436000 

884000 

0 

439000 

875000 

0 

436000 

885000 

0 

439000 

876000 

0 

437000 

875000 

0 

439000 

877000 

17.5 

437000 

876000 

17.5 

439000 

878000 

17.5 

439000 

879000 

0 

441000 

881000 

0 

439000 

880000 

0 

441000 

882000 

0 

439000 

881000 

0 

441000 

883000 

0 

439000 

882000 

0 

441000 

884000 

0 

439000 

883000 

0 

441000 

885000 

0 

439000 

884000 

0 

442000 

875000 

0 

439000 

885000 

0 

442000 

876000 

0 

440000 

875000 

0 

442000 

877000 

17.5 

440000 

876000 

0 

442000 

878000 

17.5 

440000 

877000 

17.5 

442000 

879000 

0 

440000 

878000 

17.5 

442000 

880000 

0 

440000 

879000 

0 

442000 

881000 

0 

440000 

880000 

0 

442000 

882000 

0 

440000 

881000 

0 

442000 

883000 

0 

440000 

882000 

0 

442000 

884000 

0 

440000 

883000 

0 

442000 

885000 

0 

440000 

884000 

0 

443000 

875000 

0 

440000 

885000 

0 

443000 

876000 

0 

44i000 

875000 

0 

443000 

877000 

17.5 

441000 

876000 

0 

443000 

878000 

17.5 

441000 

877000 

17.5 

443000 

879000 

0 

441000 

878000 

17.5 

443000 

880000 

0 

441000 

879000 

0 

443000 

881000 

0 
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Data  Points  for  the  Salt  River 
1,000  ft  Grid 


Grid  Coordinatea 

Sand  body 

Grid  Coordinates 

Sand  body 

in 

Faat 

Thicknaas 

in 

Faat 

Thickness 

East 

North 

in  Faat 

East 

North 

in  Feet 

441000 

880000 

0 

443000 

882000 

0 

443000 

883000 

0 

445000 

884000 

0 

443000 

884000 

0 

445000 

885000 

0 

443000 

885000 

0 

446000 

875000 

0 

444000 

875000 

0 

446000 

876000 

0 

444000 

876000 

0 

446000 

877000 

17.5 

444000 

877000 

17.5 

446000 

878000 

17.5 

444000 

878000 

17.5 

446000 

‘879000 

17.5 

45 1000 

877000 

0 

446000 

880000 

0 

444000 

879000 

0 

446000 

881000 

0 

444000 

880000 

0 

446000 

882000 

0 

444000 

881000 

0 

446000 

883000 

0 

444000 

882000 

0 

446000 

884000 

0 

444000 

883000 

0 

446000 

685000 

0 

444000 

884000 

0 

447000 

875000 

0 

444000 

885000 

0 

447000 

876000 

0 

445000 

875000 

0 

447000 

877000 

17.5 

445000 

876000 

0 

447000 

878000 

17.5 

445000 

877000 

17.5 

447000 

879000 

17.5 

445000 

878000 

17.5 

447000 

880000 

0 

445000 

879000 

0 

447000 

881000 

0 

445000 

880000 

0 

447000 

862000 

0 

445000 

881000 

0 

447000 

883000 

0 

445000 

882000 

0 

447000 

884000 

0 

445000 

883000 

0 

447000 

885000 

0 

448000 

875000 

0 

450000 

877000 

0 

448000 

876000 

0 

450000 

878000 

0 

448000 

877000 

0 

450000 

879000 

17.5 

448000 

878000 

17.5 

450000 

880000 

17.5 

448000 

879000 

17.5 

450000 

881000 

17.5 

448000 

880000 

0 

450000 

882000 

0 

448000 

881000 

0 

450000 

883000 

0 

448000 

882000 

0 

450000 

884000 

0 

448000 

883000 

0 

450000 

885000 

0 

448000 

884000 

0 

451000 

875000 

0 

448000 

885000 

0 

451000 

876000 

0 

449000 

875000 

0 

453000 

878000 

0 

449000 

876000 

0 

453000 

879000 

0 

449000 

877000 

0 

453000 

860000 

17.5 

449000 

878000 

17.5 

453000 

861000 

17.5 

449000 

879000 

17.5 

453000 

882000 

0 

449000 

800000 

17.5 

4S3000 

883000 

0 

449000 

810000 

0 

453000 

884000 

0 

449000 

882000 

0 

453000 

885000 

0 

449000 

883000 

0 

454000 

875000 

0 

449000 

884000 

0 

454000 

876000 

0 

449000 

885000 

0 

454000 

877000 

0 
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1  Data  Points  for  the  Salt  River 
|  1,000  ft  Grid 

Grid  Coordinate* 
la  Foot 

Beat  North 

Sand  body 
Thicknasa 
in  Feat 

Grid  Coordinate* 
in  Feat 

Eaat  North 

Sand  body 
Thicknaaa 
in  Feat 

450000 

875000 

0 

454000 

878000 

0 

450000 

876000 

0 

454000 

879000 

0 

|  454000 

880000 

17.5 

456000 

882000 

0 

454000 

881000 

17.5 

456000 

883000 

0 

454000 

882000 

17.5 

456000 

884000 

0 

454000 

883000 

0 

456000 

685000 

0 

454000 

884000 

0 

457000 

875000 

0 

454000 

885000 

0 

457000 

876000 

0 

455000 

875000 

0 

457000 

877000 

0 

455000 

876000 

0 

457000 

876000 

0 

455000 

877000 

0 

457000 

879000 

0 

455000 

878000 

0 

457000 

880000 

17.5 
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Data  Points  for  the  Salt  River 
1,000  ft  Grid 
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6/4/80  16.5  15.0  16.5  12/12/17  Sand,  brown  flneto  coarse, (SW) 
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bedded,  wet,  very  dark  gray 
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to  57.3-very  hard  sandstone) 
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Boring  No.  1148  Project t  RMA _ Date*  11/24/80  Location! 

Drill  Rig:  Mobile  B-53  Inspector:  R.Hunt _ Operator!  B. Harried 
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Appendix  H  Crocs  Sectional  and  Tabular  Stratigraphic  Information 


Cross-sectional  information  for  boring  746  (from  Braughton  et  al.  1979) 


Appendix  H  Cross  Sectional  end  Tabular  Stratigraphic  Information 


5260 


£  5240 
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5200 


Tabular  Data  for  all  Wells  in  the  Rocky  Mountain  Arsenal  Study  Site  (after 
Ebasco  at  al.  1989) 


Borins  # 

Section 

Zona 

Sandstone 

and 

Vail  # 

or 

Unit 

Top 

Elevation,  ft 

Base 

Elevation,  ft 

Thickness,  ft 

01005 

AS 

5201.6 

5174.6 

27.0 

01005 

AU 

5211.6 

5204.6 

7.0 

722 

01008 

AL 

5181.2 

5176.2 

5.0 

722 

01008 

AM 

5209.9 

5190.7 

19.2 

746 

01015 

AU 

5216.5 

5206.1 

10.4 

748 

01017 

AM 

5182.6 

5177.5 

5.1 

748 

01017 

AU 

5193.7 

5183.1 

10.6 

1143 

01022 

AM 

5154.5 

5147.5 

7.0 

1143 

01022 

AU 

5171.5 

5162.9 

8.6 

1143 

01023 

1U 

5107.0 

5095.0 

12.0 

1155 

01025 

AU 

5173.9 

5166.9 

7.0 

1155 

01026 

AL 

5152.4 

5146.9 

5.5 

1154 

01028 

AS 

5197.2 

5190.2 

7.0 

1154 

01028 

AU 

5206.1 

5202.2 

3.9 

1154 

01029 

AL 

5156.2 

5142.7 

13.5 

1154 

01029 

AM 

5184.7 

5181.7 

3.0 

1162 

01031 

AU 

5208.1 

5206.2 

1.9 

1162 

01032 

AM 

5184.0 

5177.6 

6.4 

1157 

01034 

AM 

5174.0 

5168.6 

5.4 

1157 

01034 

AU 

5193.5 

5192.2 

1.3 

1157 

01035 

AL 

5162.0 

5156.5 

5.5 

1236 

01036 

AU 

5202.9 

5201.6 

1.3 

1236 

01037 

AL 

5161.9 

5160.6 

1.3 

1236 

01037 

AM 

5172.6 

5165.8 

6.8 

1237 

01039 

AU 

5192.9 

5191.4 

1.5 

1237 

01040 

AL 

5165.0 

5157.1 

7.9 

1237 

01040 

AM 

5173.4 

5170.3 

3.1 

1238 

01042 

AL 

5171.2 

5168.8 

2.4 

1238 

01042 

AU 

5202.2 

5201.0 

1.2 

1238 

01043 

1 

5112.0  . 

5106.8 

5.2 

01046 

2 

0.0 

0.0 

49.0  Eat. 

1240 

01047 

1 

5093 . 8 

5074.0 

19.8 

1240 

01047 

1U 

5114.9 

5107.6 

7.3 

1240 

01047 

AL 

5157.8 

5156.3 

1.5 

1240 

01047 

AM 

5185.6 

5182.6 

3.0 

1240 

01047 

AU 

5192.3 

5189.3 

3.0 

1240 

01048 

2 

5070.0 

5045.3 

24.7 

1241 

01050 

AS 

5203.4 

5157.4 

46.0 

(AP29) 

01067 

AS 

5199.0 

5160.1 

38.9 

(AP29) 

01067 

AU 

5218.0 

5208.9 

9.1 

01068 

VC 

5283.9 

5238.6 

45.3 

01071 

1 

5092.6 

5078.7 

13.9 

01071 

1U 

5129.1 

5118.6 

10.5 

01071 

2 

5075.1 

5070.6 

4.5 

01071 

AL 

5174.6 

5172.1 

2.5 

Appendix  H  Cram  Sectional  and  Tabular  Stratigraphic  Information 


H5 


Tabular  Data  for  aN  Walt  in  tha  Rocky  Mountain  Arsenal  Study  Site  (after 
Ebasco  st  al.  1989) 


Section 

Zona 

Sandatona 

and 

or 

Top 

Baae 

Borina  # 

Veil  # 

Unit 

Elevation,  ft 

Elevation,  ft 

Thickneaa,  ft 

01071 

AM 

5186.1 

5185.1 

1.0 

01071 

AO 

5200.5 

5195.1 

5.4 

705 

02004 

AS 

5208.1 

5162.8 

45.3 

1122 

02009 

1 

5105.7 

5103.7 

2.0 

1122 

02010 

2 

5086.8 

5077.9 

8.9 

1122 

02010 

3 

5072.9 

5044.9 

28.0 

1124 

02012 

10 

5114.6 

5109.6 

5.0 

1124 

02013 

2 

5063.6 

5048.9 

14.7 

1123 

02015 

10 

5149.2 

5134.2 

15.0 

1123 

02016 

2 

5091.7 

5075.2 

16.5 

1128 

02018 

AO 

5221.4 

5208.7 

12.7 

1128 

02019 

10 

5165.0 

5159.0 

6.0 

1128 

02019 

AL 

5187.5 

5169.4 

8.1 

1148 

02021 

AM 

5182.0 

5167.4 

14.6 

1148 

02022 

10 

513S.3 

5125.5 

12.8 

1148 

02022 

AL 

5154.0 

5143.0 

11.0 

1153 

02024 

AL 

5178.3 

5177.2 

1.1 

1153 

02024 

AM 

5191.1 

5186.2 

4.9 

1158 

02027 

AL 

5153.4 

5142.9 

10.5 

1158 

02027 

AM 

5160.2 

5156.6 

3.6 

1158 

02028 

10 

5117.4 

5103.7 

13.7 

1161 

02030 

AL 

5177.9 

SL76.0 

1.9 

1161 

02030 

AM 

5196.4 

5195.0 

1.4 

1161 

02030 

AO 

5219.6 

5208.3 

11.3 

1161 

02031 

10 

5135.5 

5128.6 

6.9 

1242 

02032 

AO 

5190.1 

5181.9 

8.2 

1242 

02033 

10 

5127.2 

5102.6 

24.6 

1243 

02035 

AL 

5180.0 

5177.8 

2.2 

1243 

02035 

AM 

5197.5 

5191.9 

5.6 

1243 

02035 

AO 

5207.0 

5202.5 

4.5 

1244 

02038 

AM 

5204.1 

5190.6 

13.5 

1244 

02039 

10 

5154.0 

5147.0 

7.0 

1244 

02039 

AL 

5175.9 

5159.7 

16.2 

1246 

02041 

AL 

5179.2 

5167.7 

11.5 

1246 

02041 

AM 

5200.0 

5197.0 

3.0 

1246' 

02042 

10 

5164.0 

5143.5 

20.5 

1247 

02043 

AO 

5216.1 

5206.2 

9.9 

1247 

02044 

10 

5149.2 

5134.2 

15.0 

1247 

02044 

AL 

5185.2 

5176.6 

8.6 

1248 

02045 

AL 

5195.1 

5184.0 

11.1 

1248 

02045 

AM 

5206.6 

5194.1 

12.5 

1248 

02045 

AO 

5227.1 

5217.6 

9.5 

1248 

02046 

10 

5154.6 

5128.9 

25.7 

1249 

02047 

AS 

5218.7 

5174.7 

44.0 

1249 

02048 

10 

5138.7 

5136.0 

2.7 
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Appendix  H  CroM  Sectional  and  Tabular  Stratigraphic  Information 


Tabular  Data  for  al  Wells  in  the  Rocky  Mountain  Arsenal  Study  Site  (after 
Ebasco  at  al.  1989) 


Boring  # 

Saction 

Zona 

Sandatona 

and 

Vail  * 

or 

Unit 

Top 

Elavation.  ft 

Baaa 

Elevation,  ft 

Thickness,  ft 

6 

35001 

10 

5165.0 

5153.0 

12.0 

139 

35004 

AL 

5178.0 

5172.5 

5.5 

17 

35005 

10 

5168.0 

5148.0 

20.0 

15 

35006 

ID 

5159.0 

5133.0 

26.0 

650 

35009 

10 

5175.0 

5156.0 

19.0 

650 

35009 

4L 

5197.0 

5181.0 

16.0 

649 

35010 

1 

5148.0 

5142.0 

6.0 

1454 

35012 

1 

5145.0 

5129.3 

15.7 

1454 

35012 

10 

5163.0 

5156.0 

7.0 

702 

35015 

40 

5213.4 

5212.0 

1.4 

723 

35016 

10 

5175.0 

5156.0 

19.0 

723 

35017 

1 

5128.0 

5122.0 

6.0 

725 

35018 

10 

5172.0 

5136.0 

36.0 

725 

35019 

2 

5127.0 

5115.0 

12.0 

726 

35021 

10 

5163.4 

5143.8 

19.6 

726 

35021 

4L 

5192.5 

5181.9 

10.6 

730 

35024 

AS 

5215.8 

5178.0 

37.8 

732 

35027 

AL 

5173.6 

5166.9 

6.7 

732 

35027 

40 

5211.1 

5204.6 

6.5 

732 

35028 

10 

5147.0 

5142.0 

5.0 

757 

35030 

AS 

5210.4 

5199.9 

10.5 

816 

35032 

1 

5143.0 

5121.0 

22.0 

816 

35033 

10 

5161.0 

5143.0 

18.0 

816 

35033 

2 

5097.0 

5091.0 

6.0 

817 

35035 

10 

5188.0 

5159.0 

29.0 

817 

35036 

1 

5143.0 

5123.0 

20.0 

818 

35038 

1 

5148.5 

5138.5 

10.0 

818 

35039 

2 

5128.5 

5090.0 

38.5 

819 

35041 

1 

5137.0 

5124.0 

13.0 

819 

35041 

2 

5110.0 

5089.0 

21.0 

822 

35042 

2 

5105.0 

5084.0 

21.0 

771 

35045 

10 

5169.0 

5157.0 

12.0 

823 

35046 

10 

5168.4 

5152.4 

16.0 

823 

35046 

AL 

5173.4 

5172.4 

1.0 

651 

35051 

10 

5173.0 

5154.0 

19.0 

1127 

35053 

AH 

5200.3 

5188.3 

12.0 

1127 

35053 

40 

5212.3 

5205.3 

7.0 

1127 

35054 

AL 

5195.9 

5177.6 

18.3 

1141 

35055 

AL 

5184.6 

5182.2 

2.4 

1141 

35055 

40 

5212.9 

5205.8 

6.1 

1141 

35055 

B 

5250.4 

5231.4 

19.0 

1141 

35056 

10 

5151.0 

5127.0 

24.0 

1145 

35059 

10 

5162.0 

5148.0 

14.0 

1145 

35060 

2 

5128.0 

5121.0 

7.0 

1147 

35062 

AL 

5179.3 

5166.5 

12.8 

Appendbt  H  Cron  S actional  and  Tabular  Stratigraphic  Information 
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Tabular  Data  for  an  Wetis  in  the  Rocky  Mountain  Arsenal  Study  Site  (after 
Ebaaco  at  al.  1989) 


Section 

Zona 

Sandatone 

and 

or 

Top 

Baaa 

Borina  # 

Vail  # 

Pnlt 

Elevation,  ft 

Elevation,  ft 

Thickneaa.  ft 

11*7 

35063 

10 

5152.0 

5131.0 

21.0 

116* 

35066 

AL 

5191.0 

517*. 3 

16.7 

Hi* 

35067 

10 

5169.0 

5153.5 

15.5 

lie* 

35066 

1 

5136.0 

5115.0 

21.0 

lie* 

35066 

2 

5115.0 

5097.0 

16.0 

lie* 

35066 

3 

5093.0 

5077.0 

16.0 

lies 

35070 

10 

5156.3 

5153.2 

3.1 

1250 

35071 

10 

5135.7 

511*. 2 

21.5 

1250 

35071 

AS 

5209.6 

5181.0 

28.6 

1250 

35072 

1 

5102.3 

5093.0 

9.3 

1251 

35073 

AS 

5209.0 

5181.9 

27.1 

1251 

3507* 

AL 

5175.9 

5170.9 

5.0 

35076 

1 

5125.0 

5120.8 

4.2 

35076 

10 

5170.0 

5156.5 

13.5 

35078 

2 

5108.0 

5100.2 

7.8 

35081 

1 

5136.7 

5133.7 

3.0 

35061 

10 

5170.7 

5161.4 

9.3 

35081 

2 

5122.7 

5101.0 

21.7 

35082 

1 

5112.0 

5106.0 

6.0 

35082 

10 

51*7.0 

5136.0 

11.0 

35062 

2 

5097.7 

5091.0 

6.7 

35062 

3 

5077.1 

5044.0 

33.1 

35082 

AL 

518*. 0 

5182.0 

2.0 

35082 

AM 

5208.0 

5200.0 

8.0 

35082 

AO 

5226.0 

5224.0 

2.0 

35088 

1 

5119.0 

5108.3 

10.7 

35088 

10 

5166.5 

5147.3 

19.2 

35089 

2 

5091.5 

5063.0 

8.5 

35089 

3 

5077.5 

5046.5 

31.0 

CPI  13 

36002 

AL 

5200.5 

5197.0 

3.5 

CP11* 

36003 

AS 

5217.8 

5211.8 

6.0 

CP115 

3600* 

AS 

5221.6 

5218.7 

2.9 

CP111 

36007 

AS 

5220.3 

5214.8 

5.5 

CP110 

36008 

AS 

5219.2 

5202.9 

16.3 

CP109 

36009 

AS 

521*. 9 

5210.1 

4.8 

CP105* 

36010 

AS 

5210.5 

5201.3 

9.2 

CP  106 

36011 

AS 

5210.9 

5201.8 

9.1 

CP107 

36012 

AS 

5213.6 

5210.5 

3.1 

*0 

36020 

AM 

5222.9 

5206.9 

16.0 

COll  6 

3602* 

AL 

5205.5 

5198.7 

6.8 

RP105 

36025 

AS 

5210.6 

5204.4 

6.2 

RP106 

36026 

AS 

5213.6 

5203.6 

10.0 

RP10* 

36027 

AM 

5216.7 

5206.7 

10.0 

RP109 

36029 

AS 

5216.6 

5210.8 

5.8 

API  13 

36033 

AS 

5222.0 

5207.0 

15.0 
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Appendix  H  Crete  Sectional  and  Tebi  'ar  Stratigraphic  Information 


Tabular  Data  for  all  Walls  in  the  Rocky  Mountain  Arsenal  Study  Site  (after 
Ebasco  at  al.  1989) 


Boring  # 

Sactlon 

Zona 

Sandatona 

and 

Hall  « 

or 

Onit 

Top 

Elevation.  ft 

Baaa 

Elavatlon,  ft 

Thickness,  ft 

RP114 

36034 

AS 

5223.0 

5209.0 

14.0 

C0101 

36036 

AS 

5218.9 

5191.5 

27.4 

C010S 

36037 

AS 

5216.8 

5185.5 

31.3 

C0109 

36038 

AS 

5214.1 

5181.1 

33.0 

C0113 

36039 

AS 

5209.1 

5183.9 

25.2 

C0116 

36043 

AM 

5196.3 

5190.0 

6.3 

C0201 

36044 

AS 

5221.2 

5180.3 

40.9 

707 

36061 

AL 

5182.3 

5182.1 

0.2 

707 

36061 

AM 

5199.9 

5191.1 

8.8 

707 

36061 

AU 

5209.8 

5209.6 

0.2 

708 

36062 

AL 

5174.8 

5154.9 

19.6 

710 

36063 

AL 

5176.8 

5159.0 

17.8 

711 

36066 

10 

5146.7 

5141.3 

5.4 

711 

36066 

AL 

5169.9 

5156.3 

13.6 

711 

36066 

AU 

5216.8 

5206.4 

10.4 

36071 

AM 

5202.5 

5193.0 

9.5 

36072 

AL 

5184.8 

5170.7 

14.1 

727 

36076 

AU 

5223.6 

5205.4 

18.2 

734 

36078 

AS 

5214.0 

5208.0 

6.0  i 

734 

36079 

1 

5142.0 

5132.0 

10.0 

734 

36079 

1U 

5163.0 

5158.0 

5.0 

36081 

1U 

5145.0 

5140.0 

5.0 

36081 

AL 

5181.7 

5166.9 

14.8 

36104 

AM 

5196.5 

5176.9 

22.6 

718 

36105 

AL 

5169.9 

5162.2 

7.7 

718 

36105 

AM 

5186.8 

5169.9 

16.9 

718 

36105 

AU 

5209.2 

5201.7 

7.5 

781 

36110 

AS 

5196.6 

5193.5 

3.1 

1149 

36113 

AL 

5168.0 

5167.5 

0.5 

1149 

36113 

AM 

5201.0 

5198.0 

3.0 

1149 

36113 

AS 

5207.3 

5206.3 

1.0 

1149 

36114 

1 

5146.0 

5126.0 

20.0 

1149 

36114 

2 

5126.0 

5100.0 

26.0 

36116 

AU 

5257.8 

5244.8 

13.0 

36117 

AM 

5224.3 

5209.8 

14.5 

1160 

36118 

AU 

5209.0 

5201.0 

8.0 

1160 

36119 

AM 

5176.9 

5158.6 

18.3 

1188 

36121 

AM 

5180.6 

5174.8 

5.8 

1188 

36122 

AM 

5158.6 

5151.9 

6.7 

LM01 

36147 

1U 

5162.6 

5161.4 

1.2 

LM01 

36147 

AL 

5212.3 

5204.9 

7.4 

LM01 

36147 

AM 

5219.0 

5216.5 

2.5 

LM01 

36147 

AU 

5224.2. 

5222.7. 

1.5 

1228 

36148 

2 

5110.0. 

5090.0, 

20.0 

1228 

36148 

3 

5090.0 

5074.0 

16.0 

Appsndfe  H  Cross  Sactional  and  Tabular  Stratigraphic  Information 


H9 


Tabular  Data  for  afl  WeiU  in  tha  Rocky  Mountain  Arsenal  Study  Site  (after 
Bbaseoatal.  1989) 


rim  # 

Section 

Zone 

Sandstone 

and 

Well  « 

or 

Unit 

Top 

Elevation,  ft 

Baae 

Elevation,  ft 

Thickness,  ft 

1228 

36149 

10 

5175.0. 

5155.% 

20.0 

1228 

36150 

1 

5144.0 

5110.0 

34.0 

1228 

36150 

AS 

5223.6 

5204.6 

19.0 

1234 

36154 

1U 

5126.7 

5116.3 

10.4 

1234 

36155 

AL 

5160.1 

5156.0 

4.1 

1234 

36155 

B 

5243.3 

5231.3 

12.0 

1235 

36156 

10 

5125.0 

5117.5 

7.5 

1235 

36156 

AL 

5155.0 

5153.0 

2.0 

1235 

36156 

AM 

5199.5 

5171.2 

28.3 

36169 

AM 

5169.0 

5165.0 

4.0  , 

36170 

1 

5114.0 

5095.0 

19.0 

36170 

10 

5137.0 

5134.0 

3.0 

36170 

2 

5095.0 

5073.0 

22.0 

36170 

AL 

5158.0 

5153.0 

5.0 

36179 

1 

5141.0 

5118.0 

23.0 

36179 

10 

5163.0 

5152.0 

11.0 

36179 

2 

5118.0 

5090.0 

28.0 

36182 

AS 

5222.0 

5174.0 

48.0 

36183 

AL 

5157.0 

5143.0 

14.0 

36183 

AM 

5164.0 

5162.0 

2.0 

1242 

02032 

AL 

5170.1 

5164.1 

6.0 

Approximate 

Motet  Zone  or  Unit  AS  la  the  "target*  sand. 
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Geophysical  Logs 


Appmdbc  I  Gaophyaioal  Loga 
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•11  Legging  Co. t  Western  Well  Logging,  Inc. 
lientt  U.S.  Amy  Corp  of  Engineers 
lele  #t  1123  (SP-8)  Date:  18  July,  1980 

[•ooatiom 

epth  Logged!  147'  Depth  Drilled! 


ire  Line  Operator!  P.  O' Brian 
leelogist/Witness  t 


Fluid  Level! 


riuid  Typei  Flui 

lit  8ise! 

[aside  Casing  Diameter! 


Unit/ instrument  #i  L-4 
Engineer!  Maj  Zebell 
Time  8inoe 
Circulation! 

Cased  Interval! 

Casing  Thickness! 


’robe  #!  646 

sags!  500  (full) 
ime  constant!  1 
ihart  Scale  (CP8/in)  * 

ging  speed  (Ft/min) i  10 


Resistivity  Scale  (Ohm/ in) t  40 
Spontaneous  Potential  (mv/in) : 
Logging  speed  (Ft/min)!  10 
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Appendix  I  Geophysical  Logs 


Boring  1123  Continued 


Appandte  I  G*ophy*ioal  Log* 
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I* 


h 


•11  Logging  Co. :  Western  Well  Logging,  Inc. 
lieat:  Army  Corp  of  Engineers 

ole  #:  1124  (SP-13)  Date:  18  July,  1980 

oca t ion: 

eptta  Logged:  194'  Depth  Drilled: 

ire  Lias  Operator:  P.  O' Brian  Dnit/lastruaeat  #:  L-4 

eologist /Witness:  Eagiaeer:  Maj  Zebell 

Tins  Since 

luid  Type:  Fluid  Level:  Circulation: 

it  Sise:  Cased  interval: 

aside  Casing  Diameter:  Casing  Thickness: 

axUMX  ELECTRIC 

robe  #:  64G  Resistivity  Scale  (Ohm/in) :  40 

ange:  500  (full)  Spontar. ious  Potential  (mv/in): 

ima  Constant:  l  Logging  speed  (Ft/min) :  10 

hart  scale  (CPS/ in) : 
ogging  Speed  (Ft/min) :  10 
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Appendix  I  Geophysical  Lope 


Boring  1124  Continued 


Appendix  I  Geophysical  Lope 
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•11  Logging  Co.:  Colorado  Well  Logging,  Inc. 
ilient:  U.S.  Corps  of  Engineers 

[ole  #t  1143  (SP-16)  Date:  Nov.  21,  1980 

locations  Rocky  Mountain  Arsenal  ( Fremont  Co.)  Colorado 
•pth  Logged!  165'  Depth  Drilled!  167' 

ire  Line  Operators  Luby  unit/ instrument  #:  93 

leologist/Witnesss  Lawson  Smith  Engineer: 

Tine  Since 

luid  Type:  Bentonite  Fluid  Level:  72*  Circulation:  2  Hr 
it  8ise:  4.5"  Cased  Interval:  0-69' 

nside  Casing  Diameter:  Casing  Thickness: 


robe  #:  47 

ange:  50  (full) 
ime  Constant:  1 
hart  scale  (CPS/ in) :  10 

ging  Speed  (Pt/min) :  15 


ELEgTRIg 

Resistivity  Scale  (Ohm/ in) :  10 

Spontaneous  Potential  (mv/in) :  5 
Logging  Speed  (Pt/min) :  25 
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Appendix  I  Geophysical  Lops 


•11  Legging  Co. :  Digilog,  Inc. 

lient:  Corps  of  Engineers 

els  it  1153  (SP-9) 

ooetion:  Rocky  Mountain  Arsenal, 

•pth  Logged:  109 • 

ire  Line  operator:  R.  Bouffard 

eologist/Vitness :  Lt  Col  Zebell 

luid  Type:  Fluid  Level: 

it  Rise:  5  5/8N 

nside  Casing  Disaster:  6M 


Date:  02/14/81 
sec.  2,  T2S,  R67W 
Depth  Drilled:  109' 
Unit/Znstruaent  it  D-2 
Engineer: 

Tine  Since 
circulation: 

Cased  interval:  0-31' 
Casing  Thickness: 


robe  it  1555  Re 

ango:  10/. 5"  s 

iae  Constant:  2  Lo 

hart  Scale  (CF8/in) :  10 

ing  Speed  (Ft/nin) :  10 


Resistivity  Scale  (Oha/in) :  20 

Spontaneous  Potential  (av/in) : 
Logging  Speed  (7t/nin) :  10 


Appcndbt  I  Gtophyiiod  Logs 
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Appendix  I  Geophysical  Lope 


•11  Logging  Co. »  Digilog,  Inc. 
liont:  Army  Corp  of  Engineers 
ole  #:  1160  (SP-2) 

ooation:  R.K.A. 

•pth  Logged!  92' 

ire  Line  Operator t  c.  Jones 

eologist/Witness  t  Richard  Hunt 

laid  Type t  Fluid  Level: 

it  Rise: 

nside  Casing  Diane ter: 


Date:  June  16,  1981 

Depth  Drilled:  92.8' 
Unit/ Instrument  #:  D-l 
engineer: 

Tine  Since 
Circulation: 

Cased  Interval: 

Casing  Thickness: 


robe  #:  1555  R 

ange:  20/. 5" 
ine  constant:  2  L 

hart  Scale  (CP8/in):  20 

ing  Speed  (Pt/min) :  20 


Resistivity  Scale  (Ohn/in) :  20 

spontaneous  Potential  (av/in) :  1 
Logging  Speed  (Pt/nin) :  20 


Appsndbt  I  Geophysical  Logs 
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•11  Logging  Co. t  Digilog,  Inc. 
lients  Corps  of  Engineers 

ole  #S  1165  (N-6)  DO to l  6/23/81 

oeation:  SNA,  Adans  Co.,  Colorado 

epth  Loggodi  112 •  Dopth  Drillod:  117* 

iro  Lino  Operators  C.  Davis  Unit/Xnatrumont  #s  D- 

oologist /Witness s  Richard  Hunt  Bnginoors 

Tins  Sines 

luid  Typos  Tlttid  Levels  circulations 

it  Sisos  5  5/8"  Cassd  Interval:  0-39* 

aside  Casing  Disasters  Casing  Thickness: 


robe  #:  1482 

angos  10/.5" 
iao  Constants  2 
hart  Scale  (CPS/ in) s  10 
ging  Speed  (Ft/ain) s  15 


Resistivity  Seale  (Oha/in) s  15 
Spontaneous  Potential  (av/in) : 
Logging  Speed  (Ft/ain) s  15 
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Appendix  I  Geophysical  Logs 


•11  Loggia?  Co. s  Digilog,  Inc. 
Hants  Army  Corps  Engineers 
ole  #t  1188  (E-l) 

OOetiOBS  R.N.A. 

•pth  Logged i  81' 

ire  Line  Operators  c.  Jones 

eologist/Witnesss  Richard  Hunt 

luid  Types  Fluid  Levels 

it  Sines 

aside  Casing  Diameters 


Dates  June  16,  1981 

Depth  Drilleds  83' 
Unit/ Instrument  # s 
Engineers 

Time  Since 
Circulations 
Cased  intervals 
Casing  Thickness s 


robe  #s  1555  R 

anges  20/. 5"  8 

ime  Constants  2  L 

hart  Seale  (CF8/in) s  20 
ging  Speed  (Ft/min) s  20 


Resistivity  Scale  (Ohm/in) s  20 
Spontaneous  Potential  (mv/in) s  10 
Logging  Speed  (Ft/min) s 


Appvndbt  I  Gaophyaical  Logs 
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tell  Logging  Co. t  Digilog,  Inc. 
ilient:  Army  Corps  of  Engineers 
[ole  #!  1228  (AP-l) 

looetiont  Sec.  36,  Colorado 
pth  Logged!  168* 
ire  Line  Operator!  Nike  Hughes 
teologist /Witness!  Richard  Hunt 

luid  Type!  Fluid  Level! 

it  Sise: 

nside  Casing  Diameter! 


Date!  30  April,  1982 

Depth  Drilled! 
Unit/Xnstrument  #x  I 
Engineer! 

Time  Since 
Circulation! 
Cased  Interval: 
Casing  Thickness: 


robe  #:  1555  R 

i  10/. 5" 
ime  constant:  2  L 

hart  scale  (CPS/ in) :  10 

ing  speed  (Ft /min) :  20 


Resistivity  Scale  (Ohm/ in) :  10 

Spontaneous  Potential  (mv/in) x  l 
Logging  Speed  (Ft /min) x  20 
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Appendix  I  Geophysical  Lops 


Boring  1228  Continued 


Appvndbc  I  Geophysical  Logs 
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Appendix  I  Geophysical  Logs 


Boring  1239  Continued 


Appandbc  I  Geophysical  Logs 


117 


•11  Logging  Co.:  Digilog,  Inc. 
lient:  Corps  of  Engineers 
Ole  #1  1247  (AP-21) 

ooetion: 

ptn  Logged :  132* 

ire  Line  Operator:  Hohaus 
eologist /Witness :  Richard  Hunt 


Date:  5-27-82 


luid  Type:  Flui 

it  8ise: 

aside  Casing  Diameter: 


Fluid  Level: 


robe  #:  1430  R 

ange:  10/. 5" 

iae  Constant:  4  L 

hart  8eale  (CP8/in) :  10 

ogging  Speed  (Ft/ain) :  20 


Depth  Drilled:  132' 
us  Unit/ Instrument  #:  D-2 

rd  Hunt  Engineer: 

Time  Since 

Level :  Circulation : 

Cased  Interval: 

Casing  Thickness: 

Resistivity  Scale  (Ob  /in):  20 

Spontaneous  Potential  (mv/in) :  2 
Logging  Speed  (Ft/ain) :  20 
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Appendix  I  Geophysical  Lops 


•11  Logging  Co.*  Digilog,  Inc. 
client:  U.S.  Any  Engineers 
[Ole  #*  1251  (AP-25) 

ooetion*  Adana  Co.,  Colorado 
epth  Logged*  141' 
ire  Line  Operator*  C.  Jones 
eologist/witnesss  Richard  Hunt 

luid  Type*  Fluid  Level* 

it  Sise: 

nside  Casing  Diameter t 


Date*  5-6-82 

Depth  Drilled*  143' 
Unit/ instrument  #* 
Engineer: 

Time  Since 
Circulation* 
Cased  Interval* 
Casing  Thickness* 


robe  #:  1489  R 

*  20/. 5" 

ime  Constant*  4  L 

hart  scale  (CP8/in) *  20 

ogging  Speed  (Ft /min) *  20 


Resistivity  Scale  (Ohn/in) :  40 

Spontaneous  Potential  (mv/in) *  1 
Logging  Speed  (Ft/min) t  20 
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Appendix  I  Geophysical  Logs 


Boring  1251  Continued 


Appendix  I  Geophysical  Logs 
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Appendix  J 
Gradation  Curves 


Appendix  J  Gradation  Curvaa 
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GRADATION  CURVES 


Appendix  J  Gradation  Curve* 
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Appendix  J  Gradation  Curves 
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Standard  Deviation  Contours 
with  Data  Point  Locations 


Appondbt  K  Standard  Deviation  Contours  with  Data  Point  Locations 


Appendix  K  Standard  Deviation  Contours  with  Data  Point  Locations 
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Appsndbc  K  Standard  Deviation  Contours  with  Data  Point  Locations 
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